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ABSTRACT 

The robustness of the stability of multivariable linear time-invariemt 
feedback control systems with respect to model uncertainty is considered 
using frequency dcwiain criteria. Available and new robustness tests are 
unified xinder a common framework based on the nature cind structure of model 
errors. These results are derived using a multiveiriable version of Nyquist's 
stability theorem in which the minimum singular value of the return dif- 
ference transfer matrix is shown to be the multivariable generalization 
of the distance to the critical point on a single-input, single-output 
(SISO) Nyquist diagram. Using the return difference transfer matrix a 
very general robustness theorem is presented from which all of the robust- 
ness tests dealing with specific model errors may be derived. These latter 
robustness tests regarding the stability of the feedback system under model 
variations may be divided into two categories: (a) those that use only 
the magnitude of the model error and (b) those that use some aspect of the 
model error structure , in addition to its magnitude. The robustness tests 
that explicitly utilize model error structure are cible to guarantee feed- 
back system stability in the face of model errors of larger magnitude 
than those robustness that do not utilize model error structure and thus 
represent an improvement of these latter robustness tests. 

The robustness of Linear-Quadratic-Gaussian (LQG) control systems are 
analyzed via this robustness theory and multiloop stability margins are 
presented; in particular, a new type of margin, a crossfeed margin, is 
introduced. Other frequency dcsnain analysis and design techniques are 
also briefly discussed and their relation to the present robustness 
analysis is examined. 
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1. INTRODOCTIC»l 

The importance of obtaining robustly stable feedback control systems 
has long been recognized by designers. Indeed, a principal reason for 
using feedback rather than open-loop control is the presence of model 
vmcertainties. Any model is at best an approximation of reality, and the 
relatively low order, linear, time-inv^u:iant models most often used for 
controller synthesis are bound to be rather crude approximations. 

More specifically, a given system model can usually be characterized 
as follows. There is a certain range of inputs typically bounded in 
an^litude and in a certain frequency reuige for which the model is a 
reasonable engineering approximation to the system. Outside of this 
range, due to neglected nonlinearities and dynamic effects, the model 
and system may behave in grossly different ways. Unfortunately, this 
range of permissible inputs is rarely spelled out explicitly along with 
the model, but is rather implicit in the technology that the model came 
from - there is no "truth in modelling" law in systems theory. 

The term robustness as used in this thesis will refer to the 
extent to which a model of a open-loop system may be changed from the 
nominal design model without desteJailizing the overall closed-loop feed- 
back system designed to control the outputs of the open-loop system. 

We stress that in this definition, we implicitly assume that the dynettaic 
compensator is fixed, that is, it does not change if, for whatever reason, 
one suspects that the actual open-loop dynamics are different from those 
used in the model. Real time changes in the compensator structure (gains 
or onher changes) lead to adaptive control systems, a topic that will not 
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be addressed in this thesis. Thus, the tern robustness refers the 
preservation of closed-loop system stability in the face of model un- 
certainty not accounted for in the compensator design. 

Robustness issues are not new in control system design. In classi- 
cal single-input, single-output (SISO) servomechanism designs, robust- 
ness specifications were often specified in terms of gcdn margin and 
phase margin requirements. However, for multiple-input, multiple- 
output (MIMO) control systems, similar robustness measures are not 
strai^t forward, and their interpretation must be done with care. Thus, 
the major theme of this thesis is to address robustness issues for MIMO 
control designs. 

The robustness problem can be logically divided into three distinct 
questions: 

(a) given a model of a feedback control system how close to 
instability is it? 

(b) given the class of model errors for which the control 
system is stable, does this class include the model 
errors that can be reasonably expected for this 
particular system? 

(c) how can a robust feedback syst^ be designed? 

Question (a) is am analysis problem that can be solved exactly by an 
appropriate mathematical formulation. This problem will be addressed 
extensively in this thesis amd is by far the easiest of the three 
questions to answer. Question (b) cannot be answered without a proper 
understanding of the physics of the physical system to be controlled 
and the assvmptions that were made in constructing a model to be used 
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in controller design. Even with a good understemding of SKtdelling de- 
ficiencies it is difficult to characterize this knowledge in a fona that 
is mathonatically easy to deal with froa the analysis point of view. 
Question (c) combines aspects of both questions (a) and (b) in that a 
designer must be able to tell if there exists a controller that would 
be ^dble to tolerate the class of modelling errors he believes is reason- 
able for a given open- loop system design model. 

However, the robustness properties of a feedback system cemnot be 
optimized without regard to the deterministic and noise performeuice 
requirements for the control system. For <^n-loop stable systons, this 
is clearly demonstrated since the roost robust control system is the open- 
loop system with no feedback. Of course, for this open- loop stable systmn 
the transient response to a step input comnand or the response to dis- 
turbances may not meet the performance specifications. This underscores 
the fact that there is a f\indamental tradeoff between robustness, 
determinstic performance and stochastic performance (performance with 
respect to stochastic disturbance and/or sensor noise inputs) . Speci- 
fication of any ono of these system ch<iracteristics may place constraints 
on the achievable performance or margin of s:tability for the other two 
system characteristics. For example, with linear-quadratic-gaussian 
(lOG) regulators one may obtain acceptable deterministic responses to 
comnand inputs cmd have an adequate margin of stability but the adequate 
robustness properties may be obtained at the expense of an increased 
response to process noise driving the open- loop plant if the deterministic 


performance must be maintained. 
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In signal-lnputr single-output (5IS0) control system design these 
issues are well understood. The classical frequency domain techniques 
for SISO design naturally handle the robustness characterization^. 

These techniques »ploy varioxis graphical means (e.g. , Bode« Hyquist, 
inverse-Nyquist > Nichols diagrams) of di^laying the system sndel in 
terms of its frequency response. From these plots, it is very very easy 
to determine (by inspection) the minimum change in model frequency response 
that leads to instability. Fr(»t the same plots the system's transient 
response and response to various inputs ceui also be estimated. Thus, 
the classical control system designer can observe the fundeucental trewie- 
offs that must be made from these plots. 

This is in contrast to the multiple-input, multiple-output (MHO) 
case where these tradeoffs eure often obscured. Many design techniques 
for MIMO systems such as pole placement completely neglect the robustness 
issue in placing poles to obtain a good tramsient response. Other state 
space methods atten^t to overcome this problem by using state- space 
models whose parameters may vary and t)^n assuring that for a range of 
parameter values the closed- loop feedbaclc system will be stable. How- 
ever, these parameterized state-space models caimot characterize mcxlelling 

errors arising from neglected dynamics and, therefore, c»Bit an important 

2 

class of variations in the nominal design model for stability analysis. 

In short, many state space methods do not naturally lead to techniques 
that adequately account for modelling error. 


See the fundamental work of Bode [6] , and any good classical textbook, 
but especially [9] . 

2 

If the dimension of parameterized state-space model is allowed to increase 
then neglected dynamics could be accounted for. 
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The presently available frequency dooain NIMO design techniques 
[1,2,4,5,56] also have the probl^ that they <k> not ensure stability for 
a sufficiently Icurge class of modelling errors. They basically treat a 
MIMO system as a series of single-loop design problons that are essentially 
decoupled. They give good stability margins in a coordinate system that 
makes the design problem simple but not in the coordinate system of the 
input and output of the physical plcmt, the coordinate systan in which 
it is important to have robustness emd good stability margins . For this 
reason, these methods may not detect small modelling errors that could 
potentially destabilize the closed-loop feedback systaa. The meastires 
of the robustness of a MIMO feedback control system presented in this 
thesis do not suffer the above deficiency; they will always detect the 
near instability of a feedback control system. However, in many cases 
these robustness tests are conservative and therefore a significant 
section of this thesis is devoted to eliminating this conservatism. 

These results are derived in the frequency dom^d.n using a multivcuriable 
version of Nyquist's criterion, singular values and the singular value 
deccanposition familar from mmerical linear algebra [44] . The approach 
taken in this thesis is similar in nature to that of Doyle in [14] and to 
that of Safonov in [18] . 

1.1 Thesis Contributions 

The main contributions of this thesis are: 

(1) a simplified derivation of available and new robustness 
results for linear time-invariant systems. 

(2) the unification of these robustness results under a coimion 
framework based on a classification of various types of 
modelling errors 
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(3) the reducti(»i of conservatism of rolmistness results using 
only information about the magnitude of modelling error 
by including information about the struct\ire of the 

modelling error. ^'2 

(4) the interpretation of robustness properties of LQG control 
systons via the framework based on model error type. 

The results of this thesis sxmroarize and extend the state of the 
art on the robustness of multivariable control syst^. However, the 
practical aK>lication of these results is far from trivial and requires 
sound engineering judgment edaout the nattire of modelling errors based 
on the physics of the controlled systan . However, it is hc^ed that 
practical experience with physical systems may provide further insight 
as to how to successfully apply these new results since engineering 
knowledge about modelling errors is not always easily interpreted in the 
mathematical framework required by these results.^ 


1.2 Sunnary of Thesis 

In chapter 2 some matrix theory results that are useful in later 
chapters are collected to enable a clearer discussion of control related 


"The original motivation for e]q>loring this problem of conservatism that 
lead to the development of these results were discussions with Mr. Jenes 
Lewis, a former classmate, working on MIMD control systems for auto- 
motive engines [42]. The application of then current robustness results 
proved conservative for an engine control system similar to one that had 
worked satisfactorily for years on production automobiles. This in 
turn lead to the question of how the robustness of a control system 
may be assessed when the sufficient conditions for stability are violated. 
The nature of the solution of this problem was first suggested by Dr. 

David Castanon. 

2 

These results were further developed due to discussions with Dr. Sherman 
Chan, who raised meuiy thought provoking questions with regard to their 
practical application. 

3 

For an application of some of these results to control of multitezminal 
DC/AC power systems the reader is referred to [46], 
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robustness issues In those later chapters. The basic problsn solved in 
Chester 2 is that of finding the nearest singuleu: complex matrix to a 
given nonsingular con^lex matrix under constraints on the class of 
singular matrices considered. This is done using singular values and 
the singular value decomposition of a matrix which are scxne of the 
fundamental mathematical tools explained in Chapter 2. 

The structure of the matrix E of smallest norm that makes A-t-E 
singular, where A is a given nonsingulau: matrix, is given in the 
solution to Problem A. The main new matrix theory result is given in 
Problem B. Problem B poses the problem of finding the matrix of smallest 
norm that makes AtE singular but where E is constrained to be unlike in 
structure to the E matri.{ of Problem A. Problem C extends a special 
case of Problem B to include more complicated structural constraints on 
the matrix E. 

Chapter 3 formulates the fundamental robustness theorem (Theorem 
3.2) using a multivariable version of Nyquist's criterion from which 
all robustness tests for linear systems in this thesis may be derived. 
These robustness tests (Theorens 3.3 to 3.6 and 3.9) are formulated in 
terms of the size or magnitude of different types of modelling errors. 

They are first explained for SISO systems to demonstrate that the 
MIM3 case simply generalizes the idea that if magnitude of the change in 
the Nyquist diagram of the nominal system, induced by modelling error, 
is less tban the distance of the Nyquist diagreun to the critical (-1,0) 
point, then the closed-loop system will r^ain stable. These tests, 
employing various model error criteria are then used to formulate multi- 
loop gain, phase emd crossfeed stability margins. Corollaries 3.3 and 3.4 


give bounds on the amount of cross co\q>llng betiraen feedback channels 
that the fee^Mck system will tolerate; that is, they ^ecify crossfeed 
margins. The various robustness tests employing different model error 
criteria are then related to the well-known small gain and passivity 
theorems [12] (Theorems 3.7 and 3.8). Extensions to simple nonlinear 
systems are also given. The chapter concludes with a discussion of the 
relative merits of the various results. 

Chapter 4 begins with a discussion of how to distinguish between 
model errors that increeise the mar>jin of stability of the feedback system 
and those that decrease the margin of stability of the feedback system. 
This is first explained in the SISO case and generalized to the MIMO case 
by using the matrix theory results of chapter 2. The basic results in- 
volve defining the smallest error that destabilizes the feedback system. 
If this type of model error can be ruled out on physical grounds, the 
results describe the next smallest destabilizing model error and its 
minimum magnitude (Theorems 4.1 and 4.2). This extends the robustness 
tests of Chapter 3 enad>iing them to consider modelling errors of larg^ir 
magnitude (that violate the original tests) by eliminating only those 
model errors, of smaller magnitude, that would destabilize the feedback 
system, on the grounds that they are not physically realistic or 
plausible types of modelling errors. The interpretation of the smallest 
destabilizing modelling errors is discussed via block diagrams and 
the singular value decomposition of the return difference transfer 
matrix. An exeunple is also used to illustrate the nature of these 


results. 


-9- 


In Chapter 5 the robustness properties of LQG control systems are 
considered. Using the results of Chapter 3 and the multivariable Kalmw 
ineqxiallty the robustness properties of LQ state feedback regulators are 
derived (Theorem 5.2). The multiloop gain, phase and crossfeed margins 
for LQ regulators (Corollaries 5.1 and 5.2) and some variations (Corolleuries 
5.3 and 5.4) of the LQ regulator with better margins hold in a coordinate 
system specified by the control weighting (R) matrix. It is shown that 
if R is not selected properly the gain, phase aind crossfeed margins 
may become arbitrarily small. Using the results for the LQ state feed- 
back regulator the stability margins for LQG regulators are explored 
(Theorem 5.4). In general there are no guaranteed stability margins for 
LQG control systems unless the Kalman filter embedded in the controller 
possesses the correct dynamic model of the perturbed system amd then the 
stability margins for LQ regulators hold. This is not a practical 
assumption and robustness recovery procedures for asymptotically recover- 
ing the LQ guaranteed stability margins at either the input or output of 
the open- loop system are discussed. Next, the possibility of recovering 
stability margins at both input and output is discussed and related to 
the problem of obtaining a characterization of the expected model error. 

Chapter 6, very briefly discusses current frequency domain techniques 
for MIMO design and robustness analysis (characteristic loci, inverse 
Nyquist array and principal gain and phase methodologies) . These are 
placed in perspective with respect to the approach of this thesis. 

Chapter 7 sunmarizes the key results, gives some conclusions, and 


outlines some future research directions. 
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1.3 Notation 

The following conventions will be adopted in this thesis. All 
matrices will be denoted by ci^ital letters^, all scalars by lower 
case letters and all vectors fay underlined lower case letters. Outside 
of the chapter in vfhich they occur, all equation nvBnbers, theorem and 
corollary numbers and figure niSBbers will be prefaced with the chapter 
mmiber followed by a period wd the number occuring within the diapter. 
Thus, for example, equation (32) of Chapter 3 will be referred to as 
(32) within Chapter 3 and as (3.32) outside of Chapter 3. 


One exception to this convention is the matrix functions f(*) and h(.) 
which take matrix arguments and are themselves matrices. These functions 
are found in Chapters 3 and 4 respectively. 
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Notation 


s.t. 

subject to 

a: (A) 

trace of the matrix A 

A» 

complex conjugate transpose of the matrix A 

A* 

complex conjugate of the matrix A 

A*^ 

transpose of the matrix of A 

det A 

determinant of A 

<• ,•> 

innerproduct 

a^(a) 

i— eigenvalue of A 

a^(A) 

i— singular value of A ■ x|(A**A) 

M-ilp 

th 

p— order ni rm 

I 

identity matrix 

j 


ll-ll. 

Euclidean (or Frobenius) matrix norm 

A-l 

inverse of the matrix A 

IR 

the real numbers 

a 

the complex niimbers 

<c 

the space of nxm matrices with elements in ^ 

a e A 

a is an element of the set A 

n 

."“i 

i“l 

the product (a,a_ . . .a ) 
12 n 

|x! 

magnitude of the scalar x 

A>B 

A-B is positive definite 

A>B 

A-B is nonnegative definite 


(A,B,C) realization of the linear system specified by the 
time domain description 
X ■ Ax + Bu 
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d«fin«d as 
G(s) loop-transfsr sstrix 

8 (a) psrturbsd loop-transfer matrix 

L(s) multiplicative perturbation transfer matrix 

*^r.T open-loop characteristic polynomial 

w^.^(8) closed-loop characteristic polynomial 

<^rvT perturbed cpen-loop characteristic polynomial 

(s) perturbed closed- loop characteristic polyncmUal 

hl[Q,f(8),C) number of clockwise encirclosents of the point 0 
by the locus of f(s) as s traverses the closed 
contour C in the complex plane in a clockwise sense. 

Dj^ Nyquist contour of Fig. 3.10 

segment of for which Rets] ^ 0. 

SISO single-input , single-output 

MIMO multiple- input, multiple-output 

ORHP (CRHP) open-(closed) right-half-plane 

OLHP (CLHP) open- (closed) left-half-plane 
LQ linear-quadratic 

LQG linear-quadratic-Gaussian 

KF Kalsian filter 
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2 . MATRIX THEORY 


2. 1 Introduction 

The purpose of this chapter is to introduce in^rtont tools fron 
matrix theory and prove some results which form tite backbone of Um robust- 
ness theory of later chapters. The ^>ecific problem considered in this 
chapter is the following. Given a nonsingular c<M^lex matrix A, find the 
nearest (in seme sense) singular matrix A which belongs to a certain class 
of singular matrices. Essential use of the singular value decomposition 
for complex matrices is made in the solution of this problem as well as in 
the definition of an appropriate constraint set to which A must belong. 

It is thus necessary to review some preliminary definitions and 
properties of special complex matrices and different vector and matrix 
rK>rms. After this preliminary review and some specialized results for 
2x2 matrices the singular values of a complex matrix are defined and 
related to size of the error matrix E which is simply the difference 
A -A. Next the singular value decomposition (SVD) is presented and the 
expansion of an arbitrary matrix in the orthonormal basis generated 
by the SVD is discussed. In the final section of this chapter the structvure 
jf the error matrix E is studied via the SVD when E is both unconstrained 
and constrained to a certain set of matrices. 

2 . 2 P reliminary Definitions and Properti e r 

The following definitions euid properties are elementary and can be 
found in the many books on linear algebra [44]. 
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2.2.1 Vector and Matrix Homs 

It is useful to have a single number to measure the sixe of a vector 
or matrix. This number is called a norm 2 ind is denoted by { j < 1 1 . 

For vector norms the following relations must hold 

1 15 .I I '0 unless X » 0 (1) 

||ax|| « }a| I ix| I for any scalar a (2) 

1 1* + zl I 1 ! 1*1 1 + I Izi 1 <3> 

Three vector norms that are commonly used are given by 

1 1*1 Ip = (|*iP + 1 * 2 !’^ fP “ 1»2,») 

where x^ are the ccxnponents of x and | |n| is interpreted as max|x. | . The 

i 

norm Il*ll2 is the usual Euclidean length of the vector x. 

Two vectors x and are said to be orthogonal if their innerprodw:t, 
<x,^>, defined by 

<x,jj^> = (5) 

is zero. The Schwartz inequality, 

1 Il*ll2llzll2 

bounding the magnitude of the innerproduct, is important in solving least 
squares and minimxzn norm problems of the type we are dealing with. 

Ttirning to matrix norms, we denote the norm of a matrix A also by 1 1 a| 1 
where the following relations must hold 

1 |a| I >0 unless A * 0 (7) 

!|oa||» |a| I 1 a| I tor any scalar a (8) 
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||ab|| < 1|a|| ||b|| . dO) 

Corre^K>nding to each vector norm there is an associated induced matrix 
norm defined by 



which satisfies the conditions (7) to (10) emd is said to be subordinate 
to the vector nom. For the three vector norms given by (4) the three 
induced subordinate matrix norms are 


l|A||^=max I |a I 

j i 

I |a| '„ = ^ Uijl 

I Ihj 1^ = max (a\) 

i 

H 

vdiere X . (A A) are necessarily real as shown later, 

1 

From the definition of these norms it is apparent that 


( 12 ) 

(13) 

(14) 


I 1^1 Ip 1 I 1^1 Ip 1 1x1 Ip ; p = 1,2,00 (15) 

is satisfied for all jc. Any matrix norm which satisfies this inequality 
is said to be consistent or compatible . Another matrix norm which is used 
frequently that is compatible with the vector norm 11*112 the Euclidean 
norm. The Euclidean norm for a matrix A is defined by 

||A||g = [ E I |a (16) 

i j 

The IIaII^ norm is referred to as the spectral norm . Some useful relation- 
ships involving the spectral and Euclidean norms that cam be developed 


are 
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11*1 Isi ”‘''11*1 Ij <”> 

where A is an n x n matrix. These inequalities follow from the fact that 
A A is positive semidefinite cuid 

max X.(A**A) < | |a| « tr(A**A) < n maxX.(A^A) . (18) 

i 1 - E - ^ 1 

Also, if X is an eigenvalue of A and x is a corresponding eigenvector, 
then for consistent matrix and vector norms 

llAxIl = jXl ||xl! < llAli 11x11 (19) 

|X1<1|a11. (20) 

From this we can obtain 

1 1a1 I 2 = max X^(A*^A) < I 1 a”a1 1^ < 1 1a 1 Ij^l 1 a1 1^ . (21) 

i 

2.2.2 Special Matrices 

There are two types of matrices that will play a special role in 
the ensuing analysis. They are known as hermitian and unitary matrices 
and have special properties that make them useful. 

II 

Definition 1 ; A complex matrix A is hermiticui if A = A . 

Definition 2: A conplex matrix U is unitary if U** «= U 

Property 1 : All of the eigenvalues of a hermitian matrix are real. 

Property 2; All of the eigenvalues of a vinitary matrix have unit 

magnitude. 

Property 3 ; 1 jul * 1 if U is unitary. 
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Property 4; If a matrix A is hermitian then there exists a unitary 

II 

matrix U such that A » U A U where A is a diagonal 
matrix of eigenvalues of A. 

Note that Property 4 means that euiy hermitian matrix has a full linearly 
independent set of eigenvectors <diich are orthogonal to each other. For 
example, the columns of the U matrix in Property 4 are eigenvectors of 
the matrix A. 


II 

Definition 3; A complex matrix S is skew-hermitian if A = -A. 


Property 5 : All of the eigenvalues of a skew-hermitian matrix are 

purely imaginary. 

Property 6; If S is skew-hermitian then jS is hermitian. 

Property 7 ; The diagonal elements of a hermitian (skew-Hermitiem) 
matrix are purely real (imaginary) . 


Rayleigh's 

Principle 


If A is hermitian then 


and 


min ^ = X . (A) 

X,i0 x"x 


H 


(22a) 


= X 


max ~ = A (A) 
H max 
x/0 X X 


(22b) 


H H 

where the ratio x fix/x x is known as Rayleigh's quotient 


which achieves its minimum (maximvan) when x is an eigen- 
vector corresponding to X . (A) (X (A)). Note that 

min maix 

1 IilI I t ~ ^ always be assumed and thus the Rayleigh 

^ II 

quotient becomes simply jt Ax. 


2.2.3 Some Useful Results Involving 11‘llg ^nd |1*||2 

Any complex matrix A can be decomposed into the sum of a hermitiaui 
and skew-hermitian matrix as 
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A - + A. 


SH 


(23) 


where 




SH 


i (A -A«) 


(24) 

(25) 


A result we will use later is the following 




(26) 


vrtiich can be seen directly from the following equations 
! 1 a 1 = tr(A“A) = tr((A^ -Agy) + A^^) ] 


(27) 


tr[Aj^l - tr[Agj^] + 


- * "‘''sh*sh' 


where we have used the fact that tr(AB) = tr(BA). Some very specialized 
formulas for 2x2 matrices that are useful in deriving results in section 
2.4 are the following. If A is hermitian then 


|a1|J = X (a”a) = X_^ (A^) 

' ' ' 2 max n» 2 uc 


= max I X . (A) I ' 
i 


(28) 


or 


[Aj 1^ = max|X^(A) [ , 
i 


(29) 


For any matrix A with eigenvalues X . 


n 

tr (A) = T. X . (A) 
i=l ^ 


n 


n X. (A) 

i-1 " 


(30) 


det(A) 


(31) 
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Thus for a 2x2 matrix A 


det(Xl-A) = (X-X^MX-X^) = X - (X^+X2>X + Xj^X^ 
= X^ - tr(A)X + det A 


( 32 ) 


and the eigenvalues of A are given by 




- det A 


Thus for a hermitian 2x2 matrix we have that 

tr(A) 


^ tr(A) 2 _ 


det A 


(33) 


(34) 




When A is not hermitian, but still 2x2, we sin 5 >ly replace A by A A in 
(34) and obtain 


|a"a| 




tr(A^) ^ . ,,H 

— - - det (A A) 


(35) 


or 






(36) 


2.3 Singular Values and the Singular Value Decomposition [19,38,39,40,44,53,54] 
.’’he singular values of a complex nxm matrix A, denoted (A) , are 
the k largest nonnegative square roots of the eigenvalues of A A where 
k = min(n,m), that is 


o^(A) = Xj^“ (A*^A) 


X — 1,2,..., k 


(37) 


where we assume that a, are ordered such that o. > o. 

1 1 — 1+1 


The maximum 
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amd minimum singular values may alternatively be defined by 

I I 2 

a (A) = nutx » I I a| L (38) 

max I I I I I I I 1 2 


‘ '2 I I -1 I i-l -1 

a . (A) = min = A | if A exists (39) 

||x|i2 

The smallest singular value o . (A) measures how near the matrix A is to 

min 

being singular or rank deficient (a matrix is rank deficient if both its 
rows emd columns are linearly dependent) . To see this consider finding 
a matrix E of minimum spectral norm that makes A+E rank deficient. Since A+E 
must be rank deficient there exists a nonzero vector x such that | |x| 12 *^ 
and (A+E)x = o^thus by (38) and (39) 

‘^min^^’ - Ma^IIs = MexII^ 1 MEII 2 = • <40) 

Therefore , E must have spectral norm of at least a . (A) otherwise ."V+E 

min 

cannot be rank deficient. The property that 

o . (A) > a (E) (41) 

min max 

implies that A+E is nonsingular (assuming square matrices) will be a basic 
inequality used in the formulation of various robustness tests. The 
inequality (41) implies that 

A*^A > (42) 

which is a useful inequality for algebraic manipulation. However (42) 
does not imply (41) except when A^A and E*^ share the same eigenvector 
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for their minimvnn and maximum eigenvalues respectively. 

A convenient way of representing a matrix that exposes its internal 
structure is known as the singular value decomposition (SVD) . For an 
nxm matrix A, the SVD of A is given by 

k 


A = = I 0. (A)u.v” 

. , 1 - 1-1 
1=1 


(43) 


where U and V are unitary matrices with column vectors denoted by 


U = 






(44a) 


V = 




' 


(44b) 


and I contains a diagonal nonnegative definite matrix of singular values 
arranged in descending order as in 



n ^ m 


n < m 


(45) 


and 


= diag[a^, o^, — , ; k = min (m,n) . 


(46) 


H H 

The columns of V and U are unit eigenvectors of A A and AA respectively 
and are known as right and left singular vectors of the matrix A. 

Any unitary matrices, such as the U and V produced by ccanputing the 
SVD of a matrix, can be used to generate an orthonormal basis in which 
to exprer T an arbitrary matrix E. Let U and V be nxn unitary matrices 
with colimns as in (44) and express E as 


E E <u.v . , E> u.v 
i-i j-i 


( 47 ) 


where the innerproduct for matrices is defined by 


A 11 

<A,B> = tr(A B) 


(48) 


for complex matrices A cmd B. Note t)iat with this innerproduct the n 


H 


rank one matrices are orthogonal to each other and have unit spectral 

and Euclidean norms and thus form an orthonormal basis. The matrix 


<u^Vj , E > HjX j simply the projection of the matrix E onto the one- 

H H 

dimensional subspace spanned by u^v^ . If the elements of HiZj formed 

2 H 

into a n length vector x by stacking the n rows of same 

procedure is used to reduce the matrix E to a vector y then <u^y^» is 


H 2 

equal to the usual x ^ innerproduct between these n length vectors. 


H H 

This makes it clear that <u.v., E> u.v. can be rearramged into a vector 

- 1-0 - 1-3 


H 


(x y^)x which is just the projection of y in the direction of the vector x. 


Also, if all the matrices are formed into vectors, they will all be or- 

thogonal to each other and have unit Euclidean length. We will thus think of the 


2 2 
n rank one matrices as representing n orthogonal directions and refer 


H H 

to <u^yV, E>as the projection of E along the direction Vj • This type of 


perspective is useful in studying the structure of the error matrix 
E » A - A. 


2.4 Error Matrix Structure 


In this section we will use the tools developed in earlier sections 


to solve the problem of finding a singular matrix A nearest to 
a given matrix. T)iis can be formulated more precisely as a mathonatical 



optimization problon: 
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Problan A ’• 

min 1 |e| I 
E 

s.t. det(A+E) « 0 (49) 

In this formulation the matrix A is singly A-fE, where we refer to E as 
the error matrix . This is the simplest problem to solve since E is 
unconstrained. In what follows we make the following technical assunption. 


Assumption 1 ; The matrix A is nxn nonsingular and has distinct singular 
values . 

The assumption of nonsingularity of A assiues us of a nontrivial problem 
otherwise E is identically zero when A is singular. The assumption of 
distinct singular values is a technical one which allows us to avoid 
some combinatoric problems associated with multiple solutions. Once 
this section's material has been understood by the reader it is not dif- 
ficult to remove this assumption. 


Solution to Problem A: 

Suppose that A has the SVD given by 

A = urv” 


where 


diagfOj^ , . . . 


n ] 
n 


k+1 


U = 


H2- 


V 




V ] . 

— n 


(50) 


(51) 

(52) 


(53) 
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Now since A4E must be pingular there exists a unit vector « 


(A+E)x « 0 


and thus frc»n (38) and (39) or (11) we have that 


< 11)^1 I 2 - IIEXII 2 < llEllj 


For a minimum E L equal to 0 . (A) it is necessary that 

2 min 

arbitrary 6 that 

j9 

X » e V 


otherwise 


|A*11 >a„,„(A) . 


This can be seen by considering 


|ax| I2 “ I |ueA| 


i2 - 


and defining a unit vector z as 


A 

z * V X 


and thus 


i !*il I2 * I I2 * ^ ^ 

i=l 


unless 


T i6 

z » to, 0, ...,0,l]e , 0 arbitrary 


such that 

(54) 

(55) 

for some 

(56) 

(57) 

(58) 

(59) 

(60) 
(61) 


Therefore, Ax is given by 
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Ax ~ Av ■ r< u 


-Ev e 
~n 


j0 


and hence 


Ev 


- 0 u 
n-n 


By aimilar arguments involving the equation 
H 

X (A+E) - 0 


ore can show that 


H 

u E » -O V 
— n n-T) 


(62) 


(63) 


(64) 


(65) 


From equations (64) and (65) we can characterize the form that all solutions 
to Problem A must have, namely 


E - U 


-0 


where P is (n-1) x (n-1) and 
s 




( 66 ) 


(67) 


but is otherwise arbitrary. 

H H 

Recall from equation (47) the interpretation of < u v , E> u v 

~n — n ~n ~n 

as the projection of E onto the direction u v . From (66) we see that 

all solutions to Problem A have the same projection in the direction 

a v which we shall call the most sensitive direction since this is 
'•n — n 

the direction it is "easiest" to make A singular by changing its elements 
the "least". Note also the additional conditions that for any two 
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solutions to Problsn A say ai^ that 




(68) 


and 


_ it _ ^ 


^ H , 

< u .V , B, 

— j-n 2 


■ 0 , 


j ¥ n 


(69) 


requiring the projections of E^ and to be equal along any direction 
H H 

u^v^ and u^ ^v^ where j ■ 1,2,,. .,n. In fact, the matrix P given by 


u"ev 


( 70 ) 


2 ft 

is just the matrix of projection^ onto each of the n directions 

II 

(slightly abusing the notion of projection to mean <u,v.,E> instead of 

-i-j 

<u.v*?,E> u, V*?) that is, 

_i_:j 


p. . ■ < u.y, , E> 
ID -i-j 


( 71 ) 


Now suppose that we construct a constraint set for E so that E cannot 


have a projection of magnitude in the most sensitive direction u ^v ^ . 


This means that the matrix A+E cannot becoeie singular along the direction 
H 


u v and thus 
~ n" n 

find out just how much larger { |E 
strained optimization probloa: 


Ell^ &mst increase if A4’E is to be singular. To 


must become we formulate the con- 


Problem B; 

min ||e|| 

E 

s.t. det(A-fE) ■ 0 

TT~n ~ n 


(72) 
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Solution to Problem B ; 

The error matrix E given by 


E ■ U 


❖ Y 
Y* -4» 




where arbitary and 


|p II < «io"o , + 0(0 - o ) 

s' ' — \ n n-1 ' n n-1 


IEII2' 


where Y is given by 


Y “ \ITt +0 ,) (0 - o )e^®, 6 arbitrary 

V n-1 n 

and A has the SVD 


A - U 


n -> 




( 73 ) 


( 74 ) 


( 75 ) 


( 76 ) 


The proof of the solution to Problem B is somewhat involved and will be 
broken into several steps. However, a geometric interpretation of a 
simplified version of Problem B will be given at the end of this section 
{this is how I actually first worked the problem). Nevertheless, there 
is a need to understand the .1x2 analytical proof of the probl^ as well 
as the next simple Lemma in order to understand the geometric interpre- 
tation. 


Lenana 1 : If the SVD of A is giver, by 

A = u::v» 


( 77 ) 


k 
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then A-t-E is singular if and only if E-t-P is singular where 

P > u\v (78) 

and fu.'.thermore = IIeII^ • 


Proof ; 

A+E = U(E+P)V** (79) 

and thus 

I det (A+E) I = I det (E+P) | (80) 

since |det u| = |det v| = 1 because unitary matrices have only eigenvalues 
of the form e^^ (Property 2) . Therefore A+E is singular if and only if 
E+P is. To show I |e| I 2 = 11^112 write 

E = UPv” (81) 


then from (10) , (78) and Property 5 we have 



emd from (81) that 


and thus 


(82) 


(83) 


|e| 


2 



(84) 


The significance of Lemma 1 is that we need only consider trfe cas'^ where 
A is the diagonal matrix of singular values E, for once this probl'jm 
is solved for P, all we need do is use equation (81). Therefcie, from this 
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point on we make the following assumption. 

Assumption 2 ; The matrix A is diagonal. 

At this point we describe the steps in the solution to Problem B. 

The first step is modify Problem B in accordance with Assumption 2. The 

next step is to show that for the SK)dified Problem B, if a nonhermiti 2 ui E 

solves the modified Probl^ B, then a hermitian solution E„ also exists. 

n 

We proceed by finding all hermitian solutions to the 2x2 case emd then 
show that in the 2x2 case the solution is unique. The final step is to 
show that the 2x2 case can be extended to the nxn case. 

We now give the modified version of Problem B. 

Modified Problan B (MPB) : 

min I |e| I 
E 

s.t. det(A+E) =0, A diagonal and A>0 

|e I < 4) < o 
nn' — n 

This form of the last constraint occurs since U and V in the SVD of a 
positive definite diagonal A are both sin 4 )ly the identity matrix. 

We proceed to the next step in Lemma 2 . 

Lemma 2 : If a nonhermiticui E solves MPB then there exists a hermitian 

solution E to MPB where 
H 

Ejj = I (E+e“) (85) 


ana 


Proof; Since A-t-E must be singular there exists a unit vector x such that 


(A+E)x - 0 (87) 

and 

x”(A+E**) « 0 (88) 

II 

since A > A euid thus 

x”(A+E)x = 0 (89) 

and 

x“(A+e”)x = 0 . (90) 

Adding (89) amd (90) emd dividing by 2 we obtadn 

x®(A+E„)x = 0 . (91) 

n 

Now since A+E is hermitian, we know from Rayleigh's Principle (55) that 

H 

0 = x”(A+Ejx > X . (A+E„) (92) 

— n mxn n 

where X . (A+E„) is strictly real (Property 1) . If X . (A+E„) < 0 , then 
mm n mm h 

since A > 0 there exists a positive scalar a < 1 such that 

because the eigenvalues of a matrix are continuous functions of their 
elements . However , 

lla£Hll2l («/2) IIIeU^ + llE^il^J- allEllj < IIeU^ (94) 
which meams that | |e| could not be a minimum if 
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e < <b < o 


(95) 


n,n 


which is satisfied since 


* Re(e ) < e < $ < O 

H ' ' nn ' — ' nn ' — ^ n 

n,n 


(96) 


and thus it must be true that 


X . (A+E„) = 0 
mxn H 


(97) 


Therefore, since |e| ^ is minimum, from (94) with a « 1 we conclude that 


IVI2 = IIEM2 


(98) 


The signif icemce of Lemma 2 is that any solution E which is not hermit ian 

II 

has a hermitian part (i.e., E = 1/2 (E+E )) which is also a solution to 

H 

MPB. Thus by finding hermitian solutions to MPB we need only determine 

if E.. + E .where is s)cew hermitian, can be a solution to MPB. 

H SH SH 

Continuing our proof, we now find all hermitian solution to MPB 
for the 2x2 case . In this MPB may be restated as 


MPb» 2x2 Hermitian Case ; 

min 1 |e| 1 
E 

s.t. det(A+E) = 0 

|d| < 4> < Of 

where 



(99) 


( 100 ) 
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E** - E - 


a b 

b* d 


, a,d real 


Solution : 


E 


(j) Y 

r* 


^ere 


\{o~+<j>) (Op-^) , 6 arbitrary 


je 


Proof : Calculating llsjlp via equation (34) we have 

INIL 


ra4d I ^ ^^a+d» 2 . i. i 2 


f(^) - ad + |b| 


or 




( 101 ) 


( 102 ) 


(103) 


(104) 


(105) 


From the singularity constraint on A+E we can determine |b| as 


lb|^ - (Oj^+a) (Op+d) (106) 

and then substitute for |b| in (105) and obtain 

IIEII 2 - 1^1 + (a^+a)( 0 p+d) . (107) 

Note from (106), and the bound on |d| in (99), that 

-o^ < a (108) 

since |b|^ ^ 0. Next we calculate the partial of ||e||„ with respect to 
a and d to locate a miniouxD. This results in 


9a 2 


sgn(a+d) + 


+a) (Oj+d) 


- |sgn(a+d) + 


(¥> - °x 


+a) (a^+d) 


where sgn(’) is the sign function defined by 

1 , x>0 

sgn(x) = undefined, x=0 . (Ill 

-1 , x<0 

Note that the partial in (109) and (110) do not exist at a = -d because 
there is a i’Jnp discontinuity in their values at that point. It also 


happens that 9 




9a is never zero but changes sign at a = -d in a 


way to indicate a minimum at a * -d. To see this consider the ratio z 


given by 






(Oj^+a) (a^+d) 


which can be showTi to have magnitude less them unity by the following 
computations 


02 < 0 ^ 


( 0 -+d)o_ < 0 , (o.+d) 


(114 



-34- 


since by (99) a^-fd is positive emd thus 


or 


ad + aOj ( 02 +d )02 < ad + aO^ + Oj^Co^+d) 


( 115 ) 


+ a 2 + ( 02 +d )02 < ad + a 02 + 0 j^(O 2 +d) ( 116 ) 


(¥) - f?) 

* °2 " ^ ♦ <V> <“2«' 




j(^) * 


< 1 


(117) 

(118) 

(119) 


Now 9 ] |e| 12 / 5 ^ can be written as 


3|Wlj 1 1 

- I l8qn(a«) + il - |l+H-2l (120) 

and thus a-fd and 3 1 | E 112/93 have the seune sign which indicates that a 
global mininvBn occurs at a - -d and inplies that |a| ” |d| £ ^ < O 2 * 

By similar single eurguments, it can be shown that 9 ||e 1 1 2 / 
is strictly positive for all |a| > |d| £ ^ indicating that the c^timuB 
value for d occurs on the boundary d « Thus using (106) the value 
of b may by calculated as 


since 
and thus 


b - V (0j^+^) ( 0 ^-^) e^® 
a -d ■ ^ 

IIEII 2 -^♦^+(0^+ 4i)(02“^> 


and specify E as 


E 


(j) b 
b* -d> 


( 121 ) 

(122) 

(123) 


(124) 
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with b given by (121) . 


The next step is to show that the only solution to the 2x2 case of 
NPB is hemitian. To do this we use (36) to express ||e ||2 

1 I®1 I 2 • 2 I 1 = 1 Ie - Met El^ (125) 

\diinh is only valid in the 2x2 case. Now suppose that E is nonhermltian 
and is deconposed into 


E 


E E 
H 


SH 


(126) 


vdiere E 

H 

matrix . 


is of the form given in (124) and E 

on 

Prom (124) and (125) we note that 


is a nonzero skew hemitian 


and that from (126) and (26) we have 

M'sh'Ie 

Now using (127) and (128), IIeII^ i'’ (125) may be written as 


llEll^ 

llEll' 






'H' 


2 

2 


(127) 


(128) 


(129) 

(130) 


where c is some real scalar , the inequality in (130) follows because 

E is not identically zero. Thus, any nonhermit Ian matrix E with )iermitian 

part E„ that makes A+E singular must have a spectral nom strictly greater 
H 

than ||Ejj|| 2 - 

Continuing our proof, we must extend the 2x2 case to the nxn case. 
This is done by considering two special cases of the case where A is nx2 


36- 


and reformulating HPB to the case in which must be rank defioient. 
When A is of the form 


0 ■ 


n 

0 


A, 

2 

m 

-i_ 

0 


0 


and E is conformably partitioned with A and given by 


(131) 


E - 



(132) 


it is clear that for A+E to rank deficient it is necessary that A^ 
is singular cuid thus E^ must lae of the form given by (124) . Also since 


lUli^ - * llEjII^ (133) 

we can conclude by the same argument used in (125) tiiroi^h (130), that E^ 
is identically zero. This can be generalized to the case %dtere A is 
composed of two orthogonal vectors x and ^ so that the SVD of A is given 


hy 


where 



(134) 


(135) 


Since Z is of the form of A in (131) we have juzt considered, «re need 
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only solve for the that aekes ^ + Eq deficient with mininaB I I^qI I2 
and use the unitary matrices U and V of (134) to calcttlate the B that suUces 
A4-E rank deficient with miniravn | |e| This is of course given by 

E - UEqV** (136) 

We are now finally ready to calculate the nxn solution to MPB. Zn 
the general nxn case when A-fE is singular there exists a vector x such 
that 


(A+E)x - 0 


(137) 


with 

X ■ [Xw X ] (138) 

— i n 

where x is an n-1 dimensional vector, with llx, IL = 1, and x ,, the last 
element in x, is strictly real and nonnegative. Note tii^t x could never be such 

that x^ would be zero since that would require the last column of A+E to 
be identically zero which is inconsistent with the bounds in (99) . By 
defining a special matrix 2 we may rewrite (137) as 


where 


(AZ + E2) 


1 

X 

n 


0 


(139) 



( 140 ) 


Now note that AZ + EZ must be rank deficient and that 
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llEllIjl Il*ll2 INI 2 - llsllj 

and also that AZ has the form (134) and is given hy 


(141) 


AZ 


w 

0 ' 

-f' 



where 


(142) 


W . (CJjXj, VlVl' 


(143) 


Since | |x^| 1^ " ^ real ai^ nonnegative we have that 


Vl 


(144) 


with equality twlding if and only if 


Xj^ - [0,0,. . . ,1] 


(145) 


Using (124 ) , (142) and (144) we conclude that 


I !=z| [2 iV 

with eqiiality holding if and only if (145) holds. As will be shotm later 
11^112 “ 11^^112 bound in (141) is achieved and thus it is necessary 

that I |ez |2 is minimized requiring (145) to hold tdiich isg>lies that 



(147) 


Thus since EZ is known the last two columns of E are determineed. Using 
analogous argvnents it can be shown that the last two rows of E are also 
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conplstely dettnnined and thus E, the solution to MPB, has the fon 



( 148 ) 


where 

Y ■ ^ (o^-^) » 6 arbitrary (149) 

||pg| Ij 1 ^4>^+(Oi+'» (02“<^) (150) 

but otherwise Is arbitrary so that 

MeM^ - + (Oj+^HOj-^) - MezII^ . (151) 

Using L^naa 1 with these results gives the desired solution of Problem B 
given in (73) to (75). 


At this point we will give a geometric interpretation of a special 
case of MPB and possible extensions of constraint set for E. 


2.5 Ge<xnetric Interpretation 

A special case of MPB with ^ > 0 has a nice geometric interpretation 
using vectors. With 4> “ 0 we require that E must have no projection 

II 

ir the direction u v . If we thin)c of the columns of the matrix A, 

-n-^i 

where 


o 


1 


A ■ 


0 1 


'0 

n 


(152) 


0 
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as orthogmal vvctora than in tha 3x3 eaaa tha h matrix could ba diiplayad 
at in Pig. 1. 



Fig. It Colusm Vactora of A Matrix 

In tha 2x2 cata, KPB ainply poses tha problm of making tmo vactors parallel 
with ainimtsB "effort" with tha proviso that the original o«^>onant in the 
x-diraction of tha shortast vector must remain un^uutgad. This is illustrated 
in Fig. 2 where and are the resultant vactors tdtan ai^ 0^ xra changed 

minimally in order to align them. If there was no d bound as in the case 
in Pr^lam A the optimal changa would ba to shrink the vector to raro. 

Note that for MFB tfhe- d is rero that the magnitude of the change (i.e. « 

IIeII^ is sisg>ly tha geometric mean of the two smallest (in this case the 
only) singular values which is computed from (103) with d">0. If we now 
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<^2 


Fig. 2: Column Vectors of S - A+E 

proceed to the 3x3 case ol Fig. 1 MPB poses the problem of making the 
shortest vector parallel to some linear unitary ccxabination of the other 
two vectors witteut changing the original x-cooponent of the vector. 
The term linear unitary combina t ion of vectors is nonstandard. It means 

I 

that in a weighted sum of vectors , the weights themselves form a vector 
of umit length. The solution to this problmn is the same as in Fig. 2 
at^ the fact that there is an additional orthogonal vector makes no dif- 
ference. The solution is depicted in Fig. 3 where only the two shortest 
vectors are changed. The x vector is a unitary linear ccxnbination of the 
Oj and Oj vectors and its tip sweeps out an elipse in the y-z pl^^^e as 
the particular unitary combination changes. The "effort" required to 



Pig. 3: Solution to MPB with 

align and x is the geometric i^an of their lengths. Thus since the 
0^ vector always is shorter than x, it is the best vector in the set of 
vectors generated by unitary linear combinations of 0^ and to align 
with the cjj vector. 

The observations allow us to generalize Problan B by accomodating 
a much larger constraint set for the E matrix. This is suggested naturally 


by supposing that in addition to the constraint 


BsS ^ 1 
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in Problem B we also have the constraint 


< u ■•» E > <u , ,E > » 0 
— n-n-l -n-l-n 


( 154 ) 


which rules out a solution of the form given by (73) to (75) with 0«O. 
Thus we have ruled out what we will call "t)» worst perturbation" (i.e., 
an E on the form of (66)) as well as the "next worst perturbation", that 
is an error matrix E of the form given by (73) to (75) . We use the ad- 
jective worst because in our robustness work the singularity of A+E is 
associated with control system instability and thus the smallest error 
matrix that makes A+E singular is considered as the wor^t possible type 
of perturbation that is possible. The term "next worst" arises because 
in Problan B with 4>=0 we eliminate the worst type of perturbation from 
consideration. Now we could continue this process and eliminate the 
next worst perturbation by imposing (154) and ask what is the "next 
next worst perturbation" and so on. If we do this a nice structure of 
the "successively worst perturbations" emerges and cam be formalized 
in the following optimization problem. 


Problem C; 

min i | e| I 
E 

s.t. det(A+E) = 0 

<u.v*^,E><u.v*!*,E> = 0 for all (i, j) SO 

where A has the usual SVD given by 

H " H 

A = UEV G. (A)u.v" 

1 - 1-1 


(155) 


( 156 ) 
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Solution ; 

Let 



■ min 
(i, j)^ 



( 157 ) 


VV)i 


tu. 


H j0 
v« e-" 


H 




where 


0 


120", if k-H 

arbitrary, if k»<£ 


then 


E - + Eq 

where 

'fiTl ■ iieII, 

and satisfies the projection constraint in (155) and also 


(158) 


(159) 


(160) 


(161) 


<uv”,E->- <u.,^,E-> -0, j - 1,2,..., n (162) 

i “ k or £ 

but Eq is otherwise arbitrary. 

What the solution to Problem C formalizes is the procedure of finding 
the miniminn effort required in aligning any two colunn vectors in the Zi 
matrix or shrinking any of its column vectors and then determining which 
of these is possible given the constraints or how each of the vectors may 
or may not bo changed or perturbed. 
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To m£ike this clearer we will illustrate the solutions to the problem 
of finding the matrices E of minimun spectral norm that make A+E singtilar 
under various constraints on the E matrix. 


Exaagle : 

Let A be given by 


A = 


9 

0 

0 


0 

4 

0 


0 

0 

1 


and consider the various constraints on E. 


(163) 


Unconstrained Case: 



where i I^gll2 — ^ otherwise E^ is arbitrary. 


= 0 Case; 


E = 


'11 


0 ^ 

j0 




2e 


0 


(164) 


( 165 ) 


where 2 and otherwise e^^^ 


and 0 are arbitrary. 


0 Case: 



where 16^2 1 £ 3 and otherwise 622 and 0 are arbitrary 



but otherwise and are arbitrary. 


632! £ 1/2 Case: 
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and and 0 are otherwise arbitrary. 

It is important to point out that we have limited ourselves to 
constraints on B of a very special form and in general arbitrary con- 
straints on the form of B lead t.o a mathematical nonlinear programming 
problem that does not in general have a closed form solution. However, 
it turns out that the special form of constrednts on E will be useful in 
obtaining robustness results of Chapter IV. Next, however, we turn to the 
basic robustness problan formulation of the next chapter. 

2.6 Concluding Remarks 

This chapter has briefly introduced singular values and the singular 
value decomposition of a matrix and shown their use in finding the nearest 
singular matrix A to a given nonsingular matrix A. The main results are 
the solutions to Problems B and C which give the structure of the error 
matrix E = A - A when E is constrained to belong to a certain set. The 
norm of the matrix E is given by the geometric mean of the two smallest 
singular values of the matrix A, when E has no projection in the s\ibspace 
spanned by u v , where u and v are the left and right singular vectors 
associated with o^, the smallest singular value of A. 

These resiilts were collected in this chapter in order not to entangle 
the algebraic aspects of this problem with the robustness issues of feedback 
control systems discussed in later chapters, which utilize these results in 
the frequency domain via Nyquist's stability criterion. 
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3. BOBUSTOESS ANALYSIS FOR LINEAR SYSTEMS 
WITH UNSTRtX:TURB) NOISL ERROR 




3.1 Introduction 

The purpose of this chapter is to give a very simple ei^Ianatlon of 
how to measure the stability robustness of multivariable feedback con- 
trol systems using singular values of certain frequency response matrices. 
The difference between multiple-input-multiple-output (MIMO) , multivariable 
and single-input-single-output (SISO) fee<fi)ack control systens with re- 
spect to the robustness probl^ is illustrated by a worked example and 
scxne of the shortcomings of treating a multi-loop system as a series of 
single-loop syst@ns eure exposed. In this chapter, we assune that the 
only information we possess about the nxtdel uncertainty or model error 
is described by a single frequency dependent number which measures the 
size or magnitude of the model error. Results that u~® only error 
magnitude information 3u:e called unstructured . Those that use more 
than just error magnitude information £ure called structured . The un- 
structured robustness results of this chapter are first presented in the 
SISO case in section 3.2 for additive and multiplicative types of 
modelling errors to clearly illustrate the ideas that are later gener- 
alized in the MIMO case. 

In section 3.3 a multivuiable version of Nyquist's theorem is 
given and the worked example is given to show that although the stability 
of a MIMO system may be determined from the multivariid)le Nyquist diagram 
the stability margins for the MIMO feedback system cannot be determined 
from the multivariable Nyquist diagram. This is in contrast to the 
SISO case where the stability margins can be determined by inspection 
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of the Nyquiat diagram. In fact, in the SISO case, this is the nmin 
reason for the value of the Nyquist diagram} it determines the stability 
of a whole class of systems near the nominal system and not just the 
stability of the mmiinal system. 

In order to efficiently generalize the SISO results of section 3.1, 
a very general robustness theorem is derived in section 3.4 that forms 
the basis of the derivation of all subsequent robustness theormns in 
this chapter. This theorem is based on the idea of deforming the Nyquist 
diagram of the nominal feedback system into a Nyquist diagram of the 
actual system without changing the msnber of encirclements of the critical 
point required for stability in the multiveuriable Nyquist ^eorem. 

In section 3.5, different kinds of modelling errors are defined and 
it is sho%m that if the magnitude of these errors are bounded appropri- 
ately then the feedback syst^ will remain stable despite these modelling 
errors. It is then shown, in section 3.6, hew these errors can be in- 
terpreted from a block diagram of the pertxirbed or actual system that 
incorporates these model errors and a ccwiparison of the different theorems 
guarding against different types of errors is made. 

From bounds on the modelling errors it is shown in section 3.7.1 
how guaremteed mtiltivariedale gain and phase margins may be defined and 
determined. Section 3.7.2 introduces a new type of margin which places 
bounds on the allow^d}le anK>unt of crossfeed from one feedback channel 
to another. This crossfeed m£u;gin is also derived from boiinds on the 
modelling error obtained in the theorems of section 3.5. Using these 


robustness results the example of section 3.3 is reworked in section 3.8, 
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the stability margins are caltmlatedf and the near instability of the 
feedback systsin that was undetected by single-lo^ methods is detected 
by the methods of this chapter. 

In the section 3.9 additional rc^stness theorsM are derived 
and related to separating f\mctions . ^ese additional theorests include 
versions of the well-known small gain and passivity theorems [12]. The 
separating functions are used to show the basic similarity of the various 
robustness criteria to the sstall gain theorem. Section 3.10 gi^s seise 
simple extensions of the theorems for linear systems to the nonlinear 
case. Concluding remarks about the relationships and use of the various 
theorems is given in section 3.11. 

The major new results of this chapter are contained in Theorems 
2, 5 and 6. Theorem 2 is the general robustness theorem from which all 
subsequent theorems aze derived. Theorems 5 and 6 concern robustness 
results for modelling errors not previously considered in the literature. 
Versions of Theorems 3, 4, 7, 8 and 9 have previously ag>peared in the 
literature [12, 13, 14, 19, 47, 48, 49] and are presented herein so as 
to place in perspective the newly obtained results and s\g>port the 
explicit interpretation of the robustness criteria as bounds on the 
allowikble modelling errors. 

3. 2 Robustness emd the SISO Myquist Criterion 

The robustness of a SISO feedback system is determined by the 
distance that its Nyquist diagram avoids the critical (-1, 0) point in 
the complex plane. Suppose that we have the SISO control system of 


Fig. 1 



Fig. 1: SISO System vinder consideration 

where g{s) represents the ncaninal open-loop plant transfer function 
together with any other compensation that has been introduced. Now 
due to modelling errors the actual compensated plant is better repre- 
sented by the transfer function g{s) , a perturbed g(s). Thereforei 
we would like to know if the closed-loop system will remain sted>le when 
g(s) is replaced by g{s). This question is answered by drawing the 
Nyquist diagram of g(s) and determining if the Nyquist diagram of g(s) 
encircles the (-1, 0) point the same nxanber of times as the Nyquist 
diagram of g(s) does, (this assumes g(s) and g(s) have the same number 
of unstable poles). Suppose the Nyquist diagrams of g(s) and g(s) are 
those illustrated in Fig. 2. 

From Fig. 2 one would conclude that the perturbed closed-loop 
system is stable since the number encirclenents of (-1, 0) is unchcuiged. 
If d(u>) denotes the distance to the critical point (-1, 0) «md p{ui) 



Perturbed 


ed 

/ 


p(o^o) 


?(jw) 


Nominal 


Fig. 2: Myquist diagrams of nominal and perturbed systems 

denotes the distance between g(ju) amd g(ju), then it is afferent, from 
Fig. 2, that the closed- loop system will rmnain stable if d(U) > p(u) 
for all u). That is we could draw a graph, as in Fig. 3, derating the 
distamce to the critical point (-1, 0) for all ai and guarantee the 
stability of the perturbed closed-loop systmn if the p(u) c\irve lay 
below the d(oj) curve. 

There are several ways to define d(u)) and p(u) but the most 
natxural seems to be 

( 1 ) 


d(w) ■ |l+g{jw)( 
p{u) • |g(ju)) - g{jo))| 


( 2 ) 
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With this dsfinition of aodsl srror iht block disgrasi of psrturbsd 
cIo8sd*Ioop systsa is shown in Tig, 5 



Fig. 5: Multiplicative model error e(s). 

with this type of model error, the measures of distances, d(u) 
and p (u)) , become 


d(o)) - |l+g’^(ja))| 


(4) 


and 


p(u) - |e(ju))| 


g(jw) - g(ju) 
g( ju>) 


(5) 


Equation (4) is sis^ly obtained by letting g(ju)) ■ -1 in (5) or by using 
the additive srror robustness criterion that 


|g(joj)-g( ju) I < |l+g(jw)| (6) 

then dividing by |g(ju>)l to obtain the multiplicative error robustness 
criterion that 


g(jti)) - q(jh) ) 
g(ju)) 


< ll+g"^(jw)| - d<w) 


P (u)) ■ |e< ju) 1 


( 7 ) 
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and identifying the qitantity d(u) - {I'fg ^(ju>)| •• that which boimda the 
magnitude of the modelling error. 

In the NINO case, the multivariable analogs of t)M criteria of (6) 
and (7) will be developed using singular values as i#ell as robustness 
tests involving other types of modelling errors. The basic test is to 
tg>per bound the magnitude of some type of model error (i.e., the distance 
between the ncsninal and perturbed systems) by a general! zaf ion of the 
distance to the critical point. The )cey problem in the MIMO case is that 
these distances can no longer be measured off of a multiv^^•iable Nyquist 
diagram or a series of single-loop Nyquist diagrams. 

3.2.1 Gain and Phase Margin s 

Classically, in the SISO case, a measure of the nearness of the 
Nyquist diagram to the critical point is given by the gain and phase 
margins. These margins are defined with respect to Fig. 6 



Fig. 6: System for definition of SISO gain and phase margins 

The gain margin, denoted GN, is the largest interval (c^, c^) such that 
if a(s) ■ k, k a real constant, then the system of Fig. 6 is stable for 






nearness of a syst^ to instedsillty but nay not be accurate indicators of 
robustness in pathological cases such as the one shown in Pig. 8. 


I 



/ 

/ 

/ 

/ 

/ 

t 

I 



Fig. 8: Nyquist diagram with GM = and PM = [-180®, 180®] 

In the MIMO case it is also possible to define multiloop gain and 
phase margins which also provide an indication of svston robustness but 
do not rule out the type of situation shown in Fig. 8 appropriately 
generalized to the MIMO setting. This will be done in section 3.7.1 
but first we turn to the development of multivariable generalizations of 
the robustness tests of (6) and (7) . 

3.3 Robustness and the Multivariable Nyquist Theorem 

In this section we discuss a version of tl^ multivariable Nyquist 


theorem [1] and work a simple illustrative example that shows single-loop 
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type of stability euialysis is inadeqiiate when applied to MIK) systms. 
The feedback system to be discussed is depicted in Fi?. 9 where the lo <9 
tremsfer matrix G(s) is assumed to incorporate both the qpen-loop plant 
dynamics emd any ccnipensation employed. 



Fig. 9; Feedback system where G(s) represents the 
open- loop plant plus a cjf.pensator 

In addition, G(s) is assianed to be derived frcmi a r te space realization 
so that ^(s) a G(s) u(s) is given in the time domain by 


X = Ax + Bu 

(8) 

X = Cx 

(9) 


and thus 

G(s) =C(Is-A)"^B . (10) 

The basic issue of concern is to characterize the robustness of the 
feedback system, i.e., the extent to which the elonents of the loop 
transfer function matrix G(s) can vary from their nominal design values 
without compromising the stability of the closed-loop system of Fig. 9. 
The analysis is based on the multivariable Nyquist theorem which is de- 
rived from the following relationship 
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where 


detH+G{s) ] 




CL 


*^0L 


(s) 

Ti) 


( 11 ) 


4>qj^(s) = det(sI-A): open-loop characteristic (12) 

polynominal 

<j>^(s) = det(sI-A+BC) : closed- loop characteristic (13) 

polynominal 

and from the Principle of the Argument of complex variable theory. 

Definition (Number of encirclements) ; Let N(Q,f(s),C) denote the 
number of cloclcwise encirclements of the point by the locus of f(s) 
as s traverses the closed contour C in the ccxnplex plane in clockwise 
sense . 

A simple version of the multivariable Nyquist theorem can now be 
stated in the following form. 


Theoron 1 (Multivariable Nyquist Thee am) ; The syston of Fig. 9 is 
closed-loop stable (in the sense that 4> (s) has no closed-right -half- 

CXi 

plane (CRHP) zeros) if and only if for all R sufficiently large 

W(0, det[I+G(s)], = -P (14) 

or equivalently 

W(-l, -l+det[I+G(s) ] , Dj^) = -P (15) 

where D_ is the contour^ of Fig. 10 which encloses all P CHRP zeros 
R 


The indentations on the imaginary axis are made to include open-loop 
jio-axis poles which will be considered as unstable. 


of 

S(,e C . 


and where N(Q, 'f(s), C) is indeteminate if f(8^) ■ H for some 

Ims 



Fig. 10: Nyquist Contour Dp which encloses all zeros 

of in the CRHP, avoiding zeros on the 

imaginary axis by identations of radius 1/R. 


Notice that no controllability or observability assumptions have 
been made. If [A,B,C,1 is a nonminimal realization [50], pole-zero 
cancellations will occur when G(s) is formed, eliminating uncontrollable 
or unobservable inodes. Nevertheless, it is important to count these 
modes in the Nyquist criterion since infinitesimal changes in the matrices 
A,B, and C may make them controllable and observable even though it is 
not possible to detect the instability of these modes in terms of G(s). 
However, by using the zeros of 0 (s) instead of the poles of the loop 

transfer matrix G(s) , this version of the Nyquist theorem allows one to 


test for the internal stability of the closed-loop system. For other 
multivariable versions of Nyquist theory refer to [1-6] . 
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Remark t When compared with the classical Nyquist theoron for the SISO 
case, the multiv£u:iable Nyquist theorem is much more difficult to use, 
for two reasons. 

First, the dependence of det[I+G(s)] on the compensation is^licit 
in G(s) is complicated, and cannot easily depicted with a Nyquist, Bode or 
related plot. This fact has motivated a considerable amount of research 
on synthesis methods, e.g. , [1] - [6]. These will not be discussed at 
length since the main thrust of this thesis is primarily analysis . 

Second, and this is the key observation, one cannot get a satisfactory 
notion of the robustness of a feedback system directly from the multi- 
varicdsle Nyquist diagram. The following extremely simple example illustrates 
this fact. 


Example 1: 

Consider the linear system ^ specified by 



which is illustrated in Fig. 11. 


(16) 


(17) 


This example is a modified version of one found in 158]. 


Fig. 11: Internal Structure of Example 1 


If the feedback compensation 



is used the closed-loop system is given by 



The eigenvalues of this closed-loop system matrix are obviously -2, - 
indicating a sted>le system. The return difference matrix , ItG(s), is 
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not encircle the (-IfO) point^. 

Sv^pose now that we attonpt to Interpret this aultivariable Nyqulst 
diagram as a SISO Nyquist diagram emd read off the gain and phase margins. 

We find that a SISO system with this Nyquist diagram has an infinite 
upward gain margin, a gain reduction margin of -1/3 and a phase margin 
in the neighborhood of +106*. These margins are usually Indicative of 
a highly robust system. For exanqple, it is typically assuned that a 
^ dB gain margin (i.e. , ®i - [1/2, 2]) and a 30* to 45* phase margin 
is adequate insurance against model uncertainty within a limited band- 
width in which the model is accurate and 20 dB upward gain margin and 
+ 180® phase margin cUsove the frequency range for which the sradel is 
valid. 

In practice, stability margins for multiloop systems are often 
calculated for each feedback loc^ separately by opening one feedback 
loop at a time while keeping the remaining loops closed and determining 
the gain and phase margins for the resulting SISO systms. To make this 
clear, consider Fig. 13 where a(s) has been inserted in one of the feed- 
back channels. 

By determining the allowable values of a(s) for a(s) a real constant 
or of the form e^^ a gain and phase margin for the feedback channel with 
a(s) in it may be determined. Moving t^» a(s) to different channels, 
a gain and phase margin may be usociated with each feedback loc^. 


Hlote that the mere determination of stid>ility is accomplished more siiq>ly 
in the time domain by calculating eigenvalues of A-BC. 









Carrying this procsdure out on our exaa^ls we obtain (Mtting u , 

Cl 

Loop 1 open I loop 2 closed: 

Loop 2 c^en ; loop 1 closadi 
The Nyquist diagram for — ~ is given in Fig. 15. 


Im 



Re 


Fig. 15: Nyquist diagremi of 

Thus we see that in each feedback loop with the other held at its 
nominal value we have the following stability margins 


CM - (-1, ®) 


T 
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Again the system «K>uld sera to be hl^ily robust (using the pre- 
viously mentioned typical margin requireMnts) . In fact, if two dif- 
ferent a^(s) are inserted simultaneously (instead of on% at a time as 
in Fig. 13) in the two feedbac)( channels, the closed-lo^ system will 
remain stable if both a^(s) are such that a^(s) 6 (-1, *) for a. (s) 
real and constant and 9^ 6 (-180®, 180*) for a^(s) - e That is, 

(24) and (25^ hold simultaneously in both feedback loc^s^. 

Note, however, that the Nyguist diagrams of Figs. 12 and 15 do 
not depend on the value of the parameter b^^ thAt as l»oomes 
large the closed-loop system is close to insted>ility in the following 
sense. If the open-loop systra of Fig. 11 is perturbed slightly to 
obtain tlie system of Fig. 16, the closed- loop systra obtained by nega- 
tive Identity feedback (i.e. , u <■ -^) is unstable and has closed- 
loop poles at (l+»^)/2. This situation cannot be detected fay inspection 
of the multivariable Nyquist-diagram or a series of single- loop Nyqviist 
diagrams. It also cannot be detected by characteristic loci plots 
(4,5) which are merely polar plots of the eigenvalues of G(s) for 
8 e Dp tdiich in our case are both given by l/(stl) plotted in Fig. 15. 
Clearly, these eigenvalues do not depend upon the value of 
hence are unable, to detect the near instability problra just described. 
An example is given in reference [43] which shows also that Rosenbrock's 
synthesis procedure [1] based on diagonal dominance has similar de- 
ficiencies. This deficiency cein be interpreted as a failure to account 

^his is not true in general and is one of the deficiencies of the loop- 
at-a time method of determining stability margins; one cannot expect 
model uncertainty to only affect one lq<^ at a time! 


I PERTURBED PLANT 


NOMINAL PLANT 


Fig. 16: P*rturbatlon in Noninal Open-Loop Plant 

that oaJees cloaed-loop syat«n unstable 


for certain types of utodelling error. 

The difficulty we have uncovered can be ei^lained in the fol 
lowing way. A multivariable system will not be robust with reg>ec 
to nodelling errors if its return difference treuisfer function SMtrix 


since then a ssmII chan 


in G(ju)>.) will make !■«<;( 


singular. When this happens. 


det[ (I+G{ 1 ■ 0 and the ntsaber of encirclements of the origin counted 

in the multivariable Nyquist criterion changes. 


In this exMple, a aBall change in itG(ju^) produces a large 
change in dettl-KK ju^) ] showing that the near singularity of a eatrix 
cam^t be 'letected in terns of its determinant. Instead, tests such 
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«8 thoM devAlopad in th« following sections must be esg>loyed which 
utilise the minimum singular value to mMsure tlM near singularity of 
matrices <see equations (2,70) and (2.41)). 

3.4 Fundamental Rt^stness Characterisation 

From the example of the previous section, we can see that the 
prdolem of determining the robustness of a multivariable feedbaclc 
syst^, (i.e., its distance from instability), is uf fundamental im- 
portance. Scmie recent work in this area i« due to Safonov [7,18], who 
generalized an approach of Zames [10, 11]. Safonov's work heavily 
utilizes concepts of functior.al analysis, as is standaurd in the modem 
input-output formulation of stability theory^. However, in tlie finite 
dimensional linear-time-invariant case, a powerful robustness charac- 
terization cem be derived more simply in terms of the multivariable 
Nyquist theorem. 

In order to present the basic robustness theory from which all 
the other robustness results that work with specific model error 
criteria may be derived, we need the following notation. 

Definition ; Let 5(s) denote the perturbed loop transfer function 
Biatrix, which represents the actual system and differs from the nominal 
transfer faction matrix G(s) because of the uncertainty in the open- 
locp plant model. We will assume that S(s) has the state space reali- 
zation (A,B,d) and cpen- and closed-loop characteristic polynomials 

^See, e.g. , [12] or [13]. 


given by 


$ 0 l< 8 ) “ det(al-A) (26) 

and 

$ (s) = det(sl-A«-BC) (27) 

CL 

respectively. Furthermore, we define G(s, 6 ) as a matrix of rational 
transfer functions with real coefficients tdiic) are ccoitinuous functions 
of e for all e such that 0 < 6<^1 and for all s 6 D^, which satisfies the 
following two conditions 

G(s,0) = G(s) (28) 

2Uld 

G(s,l) = G(s) (29) 


With these definitions we may state the fuxtdamental robustness 
theorem of this chapter. This theorem does not adopt a specific model 
error criteria but works directly with a perturbed G(s) frcxn which any 
particular model error may be computed. 


Theorem 2 (Fundamental Robustness Theorem) : The polynomial $ (s) has 

ci» 

no CRHP zeros and hence the pertvirbed feedback systan is stable if the 


follcsrinq conditions hold: 


1. (a) T (s) have the same number of CRHP 

Jli Oii 

zeros 


(b) if = 0, then <|)^, (jU3„) 
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(c) (|)_ (s) has no CRHP zeros 

CL 

2. det[I-H:(s,e)] 0 for all (s,e) 6 Dj^x[0,l] for all R 
sufficiently large. 

Proof : For any 6 in [0,1] and for all R sufficiently large the contour 

Dj^ will enclosed all opcn-right-half-pleme (ORHP) zeros of 
$ (s). By virtue of condition lb and the identation construction of 

D , D will enclose all CRHP zeros of <j>-_ (s) and (s) . Also, for 

K R UL UL 

R sufficiently large, avoids all open-left-half-piane (0LHP1 zeros 
of Theorem 1 (multivariable Nyquist 

theorem) and conditio*i Ic we conclude that 

N(0, det[I+G(s,0)l , Dj^) = -P (30) 

where P is the nunber of CRHP zeros of 4 >qj^(s) and also of by 

condition la. Clearly, det[I+G (s,0) ] is a continuous function of 6 
for all s e Dj^. 

Now supp-ise that as 6 is varied continuously from zero to unity 

that the number of encirclements given by N(0, det(I+G(s,e) ] , D^^) 

changes. Since det (I+G(s,e) ] is continuous in (s,6) in Dj^x[0,l], its 

locus on D forms a closed bounded contour in the complex plcuie for 
R 

any e in [0,1] . The only way to change the number of encirclements 
of the critical point (0,0) is for the locus for some 6 in [0,1] to 
pass through the critical point, that is for son« 6^ in [0,1] 

det[l+G(s,e^)] = 0 (31) 

R" 


for seme s in D, 
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Condition 2 eliminates the possibility that det[l4G(s,6Q)] equals 
zero. 'This contradicts the assus^tion that N(0,det[I-K>(s,6) ] < D^) 
changes as 6 is V2u:ied from zero to unity, and thus it must be true that 
it retains constant at -P for all 6. However, this fact along with (20) 
imply that 

N(0, det[I+5(s)l, Djj) = -P (32) 

• nd thus by condition la and ^Theorem 1 (Nyquist's theorem), 
no CRHP zeros. 

Q.E.D. 

Remark; The basic: idea behind this theoren is that of continuously de- 
forming the Nycjuist diagram for the ncjminal system G(s) into one cor- 
responding to the Nyquist diagram of the perturbed or actual system 
G(s) without changing the nvanber of encirclements of the critical point. 

If this can be done and the manber of encirclmnents of the critical 
point recjuired for G(s) and G(s) are the same, then no CRHP zeros of 
(s) will result frcsn this perturbation. 

CJj 

Imbedding arguments of this type have been previously used, 
implicitly by Rosenbrock [1] and explicitly by Doyle (14), in connection 
with linear systems and in the more general cont 5Xt of nonlinear and multi- 
dimensional systems by DeCarlo, Saeks and Murray [15] - (17), utilizing 
homotopy theory from algebraic topology. 

Ranark ; The significance of Theorem 2 is that various multivariable 
robustness characterizations can be stated in terms of conditions that 
guarantee condition 2 is satisfied. In checking condition 2, it is 
unnecessary to consider all s 6 D if | |G(s,e) | | -♦■ 0 as |s| 

K 2 
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This will be the case in what follows and it is related to the assunption 
that the state-space realizations of G(s) emd G(s) have no direct feed- 
through from input to output so that | |g(s) | A nd | |g(s) | approach 
zero as |s| ->■ 0 °. It is therefore convenient to define the segment 
^ as 

Op = {s|s e Dp and Re(s) £ o} (33) 

which is the only part of the Nyquist contovir Dp on which condition 2 
need be verified. 

3.5 Robustness Theorems and Unstructured Model Error 

In this section, we develop theorems that guarantee the stability 
of the perturbed closed-loop system for different characterizations of 
model uncertainty (i.e., different types of model error). This is done 
via Theorem 2 by using a specific error criterion to construct a tremsfer 
matrix G(s,e) continuous in 6 on Dp x[0,l] that satisfies (28) and (29). 
Then a single test bounding the magnitude of the error is devised which 
guarcintees that condition 2 of Theorem 2 is satisfied. This procedure 
is carried out for four different types of errors. These tests use 
only the magnitude of the modelling error amd do not exploit any other 
characteristics or structure of the model error and hence are based on 
the unstructured part of the model error. These different types of model 
errors will emphasize different aspects of the difference between 
the nominal G(s) and 6{s) and thus under certain circxanstances will 
give essentially different assessments of the robustness or margin of 


stability of the feedback control system. 
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Prdbably the most familiar types of errors are those of absolute 
and relative errors. Absolute errors are additive in nature whereas 
relative errors are multiplicative in nature, (hie can use both types 
of errors to derive robustness theoresis. However, the fasilliar notions 
of gain cuid phase margins are associated only with relative type of 
error since these mcurgins are multiplicative in nature. 

If we let the matrix E(s) generically denote the p^u^ticul^u: 
modelling error under consideration, the absolute error is obviously 
given by 

E(s) » 5(s) - G(s) (34) 

cuid the relative error, in a matrix sense, by 

E(s) = g"^(s) [G( s)-G(s) ] . (35) 

In (35) G ^(s) could post-multiply the absolute error emd serve as ui 
alternative definition of relative error in the matrix sense but all 
subsequent results will still hold with trivial modifications. Using 
these errors we will prove two rotnistness theoroas. However, first 
G(s,6) must be constructed. 

Using (34) and (35) can define G(s,e) by replacing S(s) in 
(34) and (35) by G(s,€) and E(s) by 8E(s) cutd solving for G(s,e). If 
we do this we obtain 

G(s,e) - G(s) + eE(s) (36) 

where E(s) is the absolute error given by (34) or 


G(s,e) * G(s) [I+€E(s)) 


(37) 
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where E(s) is the relative error given by (35). Both (36) and (37) inply 
the same G(s,6) although they eo^loy different types of errors to eurrive 
at G(s,e). In either (36) or (37) G(s,6) is simply given by 

G(s,e) = (l-€)G(s) + e5(s) (38) 

showing that G(s,6) is continuous in 6 for 6 on [0,1] and for all 
s G and that G(s,G) satisfies (28) and (29). 

In deriving stability margins based on theorems i-.'sing different 
error criteria, we will find it useful to define a multiplicative un- 
certainty matrix L(s) to accovint for modelling errors in the open- loop 
plcint. The perturbed or actual system 2(s) in this case is given by 

5(s) = G(s)L(s) (39) 

which implicitly defines L(s) . Notice that for the relative error 
criteria that L(s) is very simply given by 

L(s) - (1+E(s)) (40) 

where E(s) is given by (35). However, as will be shown later (40) is 
not the only description of L(s); there are other types of relative 
errors yet to be discussed in which the relationship between L(s) amd 
the generic E(s) is not so simply given by (31). We will use both L(s) 
as defined implicitly in (30) and a variety of error matrices denoted 
by E(s) in stating the subsequent robustness theorems. 

Two robustness theor<3ns based on the preceedings definitions of 
absolute and relative errors in (34) and (35) respectively are the 


following. 
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Theorem 3 [48,49]: The polynomial $^(s) has no CRHP seros and hence 

“ ■ CL 

the perturbed feedback system is stable if the following conditions 
hold: 

1. condition 1 of Theorem 2 holds 

2. a_. [I+G(s)l > a^.^[E(s)] for all s e JL 

mxn max k 

where E(s) is given by (34), and was defined 
by (33). 

Proof : From (36) we see that I+G(?,e) is given by 

I+G(s,e) = I+G(s) + eE(s) . 

Frcxn the properties of singular values (see (2.41)) we know that 
I+G(s) + eE(s) will be nonsingular if 

0 . [I-Kj(s)] > a reE(s)] - SO [E(s)} (41) 

min max max 

which is clearly guaranteed by condition 2 since € is always between 
zero and unity 2 uid thus condition 2 of Theoron 2 holds. 

Q.E.D. 

Theorem 4 [14,48,49] : The polynomial $ (s) has no CRHP zeros and hence the 

” CL 

perturbed feedback system is stable if the following conditions hold; 

1. condition 1 of Theorem 2 holds 

2. 0 . [I+g“^(s)] > o [E(s)] for all s e (L 

min max a 

where E(s) is given by (35) . 


Proof: From (37) we see that I+G(s,6) is given by 
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l4G(s,e) - I-K3(s) [Z-feE(s) ] 

= G (s) [I+ g"^(s) +eE(s)] . (42) 

Here by writing G ^(s) we assume it exists^ so that XtG(s,e) is singular 
if and only if [I+G ^(s) + 6E(s)] is singular. As in the proof of 
Theorem 3, we know frcan (2.41) that condition 2 guareuitees that 
I+G ^(s) + 6E(s) is nonsingular, hence Theorem 2 is satisifed. 

Q.E.D. 

Theorem 4 was first proved by Doyle [14] using singular values 
and Nyquist’s theorem but under the slightly stronger condition that 
E(s) be stable. An operator version of Theorem 3 is due to Sandell 
[48] who was the first to consider additive perturbations. Laub [49] 
provides further numerical insights to the relationship of Theorems 
3 and 4. 

Before we give some discussion of these theorens and some possible 
corollaries, we will develop seme additional robustness theorems which 
are complementary to Theorems 3 and 4 and are derived on the basis of 
alternate definitions of the error matrix E(s). 

Suppose that instead of measuring the absolute relative errors 


The assumption that G exists guarantees that any pertrubed system 
G can be represented as 5 = G(I+E). However, if G is singular but 
S is in the range space of G then E may be implicitly defined as a 
bounded solution of GE = G-G. In this case Theorem 4 still holds 

if n . (I+G'l) is replaced by its equivalent for all G, [G(I+G) ] , 

r.un r j -I ^ 

which is bounded by condition 1 if G ? 0. If G is not in the range of 
G then G cannot be represented as 2 * G(I+S). 
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between &(s) wd G(s) , we neasiire the absolute and relative errors 

between fi~^(s) and G~^(s) . In the SZSO case, this would correspond 

to measuring the absolute and relative errors bet«reen the noninal and 

perturbed systems on em inverse Hyquist diagram in which the inverse 

-1 ~-l 

loop transfer functions g (s) and g (s) are plotted. (The Inverse 
Nyquist diagram cam also be used to determine stability by counting 
encirclements of the critical points (0,0) and (-1,0) in the complex 
plane. )^ Therefore, it would be natural to define the absolute and 
relative errors between the nominal and perturbed systems as 

E(s) * e“^(s) - g'^(s) (43) 

for the absolute error and 

E(s) » [g"^(s) - g"^(s)]G(s) (44) 

for the relative error. Using (43) and (44) we may define a G(s,6), 
again by replacing G(s) by G(s,6) and E(s) by eE(s) in (43) and (44), 
and then solving for G(s,6). If this is done, we obtain 

G(s,e) = (g"^(s) + eE(s))'" (45) 

where E(s) is given by (43) and 

G(s,e) « G(s) (I+€E(s)]’^ (46) 

where E(s) is given by (44). Both (45) and (46) give the same G(s,0) 
which written in terms of G(s) amd S(s) is 

It is not intended to give a discussion of the inverse Nyquist criterion 
[1] but only mention it to suggest that the use of G~l(s) is as reasonadsle 
as G(s) in a definition of model error. 
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G(s,e) « t(i-e)G’^(8) + eB"^{*)r^ (47) 

where now we see that 6 enters nonlinearly and it is not clear that 
G(s,6) is continuous in e in [0,1] for all sen but is is clear that 
it does satisfy (28) and (29). The type of G(s,e) in (47) could be 
replaced by the one in (38) emd theorems worlced out in terms of the 
errors described by (43) and (44). This approach was taken by Lehtomaki, 
Sandell and Athans [51] and led to more restrictive and complicated 
conditions to check than the approach using (47) . 

Since (45) , (46) and (47) are all equivalent in that they give 
rise to the same G(s,6) we may work with any one of them to prove 
assertions eibout the continuity of G(s,6) required by Theorem 2. If 
2 ^(s) and G ^(s) exist, so that E(s) in (44) is well-defined, then we 
can see that for G(s,6) to be continuous in 6 for (s,e) 6 x[0,l] all 

we need to guarantee is that [I4^(s)] is nonsingular . Notice that in 
this case L(s) is simply 

L(s) = [I+E(s)]’^ (48) 

jmd that [I+€E(s)] is nonsingular for all 6 in [0,1] if L(s) defined by 
(39) has no strictly negative eigenvalues. This is true since 

if L(s) has no zero or negative eigenvalues, neither does I-tE(s) and 
thus E(s) cemnot have eigenvalues in the interval (-“’, -1] so that 
SE(s) never has eigenvalues of -1. Therefore with these restrictions 
G(s,e) is continuous in 6 on D x[0,l]. We also see from (46) that 
if E(s) is bounded (i.e., 2 ^(s) and 2 ^(s) exist) and L(s) has no 
rero or negative eigenvalues that | |g(s,€) j 0 as |s| for any 
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e in (0,1] . This allows us to check for the nmsingularity of X-fG(B,€) 
only on Qj^x(0,l]in Theorem 2. We may now state the theorems analogous 
to Theorems 3 emd 4. 

Theoran 5; The polynomial (s) has no CRHP zeros and hence the 



perturbed feedback system is stable if the following conditions holdt 

1. condition 1 of Theorem 2 holds 

2. L(s) of (39) has no zero or strictly negative real 
eigenvalues for any s 6 (2^^ 

3. 0 [i-k; ^(s)J > o „tE(s)] for all s e 

mn inaX R 

where E(s) is given by (43) 

Proof ; From (45) we have that 

I-K3(s,e) - I + [g"^(s) + eE(s)]’^ 

- (I+g‘^(s) + 6E(8)] (g"^( 8) + eE(s)]“^ 

- (I+g"^(s) + eE(8)]G(8,e) (49) 

and since G(s,6) is nonsingular^ , I-fG(s,<3) is nonsingular if and only 
if [I-fG ^(s) * eE(s)] is nonsinguiar which is trxie by condition 3. 
Condition 2 merely ensures that we have a G(s,6) continuous in 0 to 
work with as required to apply Theorem 2. Thus, Theorm 2 holds 
and 2 ( s) has no CKHP zeros . 

CL O.E.D. 


In this proof essential use of the fact that G(s) and S(s) are both 
invertible on is made. This is different from the case of the 
footnote of Theorem 4. 
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Th« next th«oran work* with th* relative error between fi ^(*) and 
G ^(s) and play* a fxuidaaental role in e*tabli*hing the propertie* of 
LQ (linear-quadratic) state feedback requlator* which will be discussed 
in Chapter 5. 

Theorem 6s The polyn<»nial (s) has no CRHP zeros and hence the 

CL 

perturbed feedback system is stadale if the following conditions hold; 

1. condition 1 of Theorem 2 holds 

2. L(s) of (39) has no zero or strictly negative real eigenvalues 

3. o . tl+G(s)] > o (E(s)J for all s 6 (L 

nxn BiAX R 

where E(s) is given by (44) 

Remark: If condition 3 is satisfied and 0 . [I+G(s)] < 1 then it can 

min — 

be easily shown via (46) that condition 2 is autcsnatically satisfied. 
Proof ; Prom (46) we have that 

1+G(s,e) * I + G(s) [I+€E(s) - (I+G(s) + 6E (s) ][I+€E (s) ] (50) 

and condition 2 not only ensures G(s,e) is continuous^ on Dj^x(0,l] but 
also that I+€E(s) is nonsingular on the same set. Thus I-fG(s,6) is 


In this proof no essential use of the fact that 5 and G exist is 
made. If E is implicitly defined by G ■ 5{T+E) rather t an (34), then 
Theorem 6 still holds. However, if d is not in the range space of G, 
and vice versa, it is not possible to represent 5 as 5 = G(I+E)”1. 
Thus, even if G“^ and d“l do not exist it may be that G(s,6) is 
continuous on Dj^x[0,l] by using the implicit definition of E if d 
can be so represented. 
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nontingular if and mly if tX4C(s} i- 9B(s)] is nmisingular idiich is 
guarantssd by condition 3. Hanes, again Tbaorssi 2 is satisfiad and 
tharafore ClUff zaros. 

wL 

Q.B.D. 

Theoren 6 is an inpzovad version of a thaoram fowkl in (51). 

Obsarvation ; Tha ctmdition that L(s) hava no strictly raal and nagativa 
aigenvalues or be singular can be intarpreted in tarns of a phasa 
reversal of certain signals between the n<minal and perturbed systasui 
or as the introduction of transnission zaros by the siodelling error. 

To make this precise, suppose that for sMie that L(ju)^)x • \x for 
some covplex nonzero vector x and some real X £ 0. Then there exists a 
vector u(t) of input sinusoids of various phasing and at frequency 
which when aj^lied to the ncninal syst^ produces an output ^(t) and 
produces an output X ^(t) when applied to the perturbed systan. This 
is depicted in Fig. 17. 






*Xy(t) 


Fig. 17: Relationship between ncninal and perturbed 

system for special input ^(t) %dtan L(ju) ) 
has eigenvalue X. 
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ThuB when X it negative the phase difference between the sinueoidal outputs 
of the nominal and perturbed systems is 180* . If X«0 then the perturbed 

system has trannnission zeros at ijw.. 

Q 

Alls fact is significant since Theorbos 5 and 6 can never guarantee 
stability with respect to model uncertainty when the phase of the system 
outputs is completely uncertain above some frequency or with respect 
to sensor or actuator failures in the feedback channels. 

Note that in general that condition 1 of Theorems 3 to 6 and the 
state-space description of the nominal and pertxirbed systems implicitly 
place other restrictions on L(s) (and also E(s}). These are simply that 
L(s) represent a finite dimensional linear time -invariant system that is 
possibly unstable and that L(s) has no purely imaginary poles. Similar 
conditions may be derived for each of the four forms of errors used in 
Theorans 3 to 6. 

3.6 Interpretations of Robustness Theorems 

Up to this point, it is prob^U3ly unclear to the reader what the 
significance of the various error criteria are and how they are related. 
This can be partly clarified by an understanding of how each error 
enters into the structure of the perturbed system from a bi.ek diagrmn 
perspective. This is done in Fig. 18 where a very pleasing s^sonetry 
occurs that corresponds to the four basic arithmetic operations of 
addition, subtraction, multiplication and division. As can be seen 
from Fig, 18 the absolute type of errors correspond to addition and 
subtraction whereas the relative errors correspond to multiplication 
and division. Other types of errors can be represented as combinations 


Block Diagram 
of Perturbed System 



Feedforward 

(Addition) 



Feedback 
( subtraction) 


Error Criterion 
Perturbed Systems and 
Stability Test 


E(s) - G(s) > 6(s) 
G(s) - G(s) -t- E(s) 


o . (I+G(s)) > 0 „.^(E(s)) 
min max 


E(s) - G ^(s) - g"^(s) 


G(s) = (g‘^(s) +E(s))"^ 


o . (I+G "^(s)) > a ,^(E(s)) 
min max 



(Multiplication) 


(Division) 


E(s) = G (s) (G(s)-G(s)] 
G(s) = G(s) (I+E(s)) 


O . (I+G (s)) > O (E(a)) 
min max 



E(s) = [g"^(s)-g‘^(s)]G(s) 


G(s) = G(s) (I+E(s)) 


O . (I+G(s)) > O ^(E(s)) 
nl.n n&x 


Fig. 18: Physical Representation of Perturbed Models 

Corresponding to Various Error Criteria and 
Associated Stability Test. 


C -X 
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of these basic types of errors. 

For now however, we will defer that t<^ic and discuss the inter- 
pretations of the robustness theorens that dead with relative errors 
ai»3 give sene pictorial illustrations of why these theorens ensure that 
I4G(s,e) is nonsingular amd how they aure related. 

We shall work mainly with the relative error type theoresis since 
from them we may derive gain amd phase margins for which design engimers 
have a more intuitive feel. In the theorems dealing with absolute 
type errors it is difficult to account for the effect of the compensator 
ing>licit in G(s) on the model error (i.e. , the model error depends on 
the coB^nsator used) . This does not happen with the relative error 
criteria. 

To begin with, recall that in Theorems 4 and 6 that L(s) is given 
respectively by 

L(s) >= I + E(s) (51) 

for E(s) given by (35) and 

L(s) = (I+E(s))‘^ (52) 

for E(s) given by (44). If we solve these last two equations for E(s) 

we obtain from (51) 

E(s) = L(s) - I (53) 

and from (52) 


E(s) - l"^(s) - I . 


(54) 
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Anis by making ||e(s )||2 small in (53) L(s) is kapt doss to tha idantity 
matrix whereas in (54) it is L ^(s) that is kapt close to tha idantity 
matrix by making | |e(s) 1 1^ small. Biis points out the difference in the 
types of errors since the sane L(s) may make one error quite large 
while making the other <»ily soderately large. 

The basic inequalities in Theorems 4 and 6 written in terms of 
L(s) are given respectively 1:^ 

a . II+g‘^(s) 1 > O [L(s)-Il (55) 

min max 

for Theoran 4 and by 

o . Il+G(s)l > a [l"^(s) - I) (56) 

min max 

for Theorem 6. The inequality (55) is the MIND generalization of tl» 

SISO inequality (7) of section 3.2 but written in terms of £(s) rather 

than g(s). Thus in (55) we see that o . [I+G ^(s)] is just the multi- 

min 

variable version of the distance to the critical point (0,0) and 

o [L(s)-i] is just the generalization of the distance between 8(s) 
max 

and G(s). Similar interpretations of (56) can be made. 

The SISO analogs of (55) and (56) are given by 

ll+g ^(s)| > a > |£(s) - l| (57) 

and 

|l+g(s)l > a > |£"^(s)-l| (58) 

respectively. In the form using £(s) rather g(s), the inequalities 
(57) and (58) provide a geometric insight to the relationship of g(s) 
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and i{s) . The adnisszd>le region of the ccnplex plane for -g ^(s) and 
£(s) satisfying <S7) is depicted in Fig. 19. Fig. 20 gives the 
analogous regions for g(s) and -£ ^(s) satisfying (58). Frost Fig. 19 
it is clear that 

-q'^(s) £(s) (59) 

and from Fig. 20 that 

-£"^(s) g(s) (60) 

which simply ensure that l-rg(s)£(s) = l-f^(s) # 0 or that g(s) does not 
pass through (-1,0). Recall that in Theoron 2 iK>t only must l+g(s) # 0 
but we must be able to construct a g(s,e) such that l-fg(s,6) ^ 0 for 
6 in [0,1]. However, due to the way that g(s,e) was constructed this 
merely results in the requiranent that 

|l+g"^(s)| > e!£(s)-l| (61) 

in the case of Theorem 4, auid 

jl+g(s) 1 > e |r^(s) - li (62) 

for Theorem 6. These inequalities are obviously guaranteed by (59) 
and (60) since 6 is between zero and unity. 

The main point of this discussion was to show the use of circles 
to divide the complex plane into disjoint regions, one in which il(s) 

(or £ ^(s)) lies and its complement in which -g ^(s) (or g(s)) must lie. 
The fact that the radii of the circles can be interpreted as the magni- 
tude of an error or the distance to the critical point is not crucial. 
Later an in this chapter we will use the idea of separating the complex 
pleme into tvo disjoint regions to derive additional robustness theorems 



Fig. 19; Admisseble regions for g > (s) and £{s) 
satisfying |l+g”^(s)| > a > |]l{s)-l| 


Im 



Pig. 20: Admissed>le regions for g{s) and -Jl”^(8) 

satisfying |l+g(s)l > a > |)l”^(8)-l| 


and clarify their relationships to the well-hnown snail gain theorra 


[ 12 ]. 


To continue the discussion on the relationship of Theorens 4 and 


6, we make the following observation that in (55) and (56) as 0 ^^^(I+G) 

and a , (I-k: ^) increases the bounds on L(s) and the error becones 
min 

less stringent. Therefore* to tolerate both kinds of modelling errors, 

one would like to medce both o . (I-H3) and a . (I4-G ^) as large as 

min min 

possible. However, these two quantities cure related algebraically so 
that we cannot make them both independently large. Their algebraic 
relationship can be derived trivially from the matrix identity [49] 


(I+G) ^ + (I+g"^)“^ E I 


using the triangle inequality auid the simple relationship O . (A) > 

mm 

o ^ (A ^) . There are three inequalities relating 0 , (I+G) to O . (I+G ^) 
max ^ rain min 

which are given by 


(I+G) + o"!' (I+g”^) > 1 
min min — 


(I+G) + 1 > (I+g“^) 

min — min 


a~] (l+G~^) + 1 > o"l- (l+G) 
min — min 


Two other inequalities relating 0 . (G) and o (G) to a . (I+G) eund 
^ min max min 


a . (I+G ) are given by 
min 


0 . (I+G) 

0 (G) > -SHI > a . (G) 

“ 0 . (I+g'^ " 


0 


1 


2 





Fig. 21: Shaded Area Represents Allowable 

Values of (o . [1+G] » o . [I+G"^l) 
nin min 

ordered pairs. 

Prom Fig. 21 it is clear that when o^.^(G) is large (i.e., large loop 

min 

gain in every feedback loop) that is necessarily neeur one 

and o . (I-Ki) is large. This indicates that Theorem 6 will give a 
min 

better indication of control system robustness with respect to the model 
error criterion (44) in the typically high performance ~ low frequency 
region than will Theorma 4. Likewise when o (G) is small (i.e., all 
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feedback loops are rolled off) Theorem 4 gives a better indication of 

the robustness of the system with respect to the model error criterion 

(35) since o . (I-fG is large and o . (Z-K3) is near unity. The exact 
min mm 

sense in which one theory gives a better robustness indication de- 
pending on the nature of G(s) will be made precise in the corollaries 
of Theorems 4 and 6 that specify different types of stability meurgins 
discussed in the next section. 

3.7 Multiloop Stability Margins 

J 

In this section we shall derive guarmteed minimum gain and £diase ^ 

margins for MIM3 systems as functions of both o^^[I+G(s)] and 

a . [I+G ^(s)]. He shall also introduce the notion of a crossfeed 
min 

tolerance which again is specified by o , (l+<3(s)] or a . [I+G ^(s)]. 

— ■ min min 

These stability margins are simply corollaries to Theorans 4 and 6 and 
are easily obtained by assiiming specific forms for L(s). 

3.7.1 Multi loop Gain and Phase Margins 

In contrast to the SISO case, it is not clear what gain and phase 
margins are in a multiloop system since gain or phase changes in one loop 
may affect the calculation of the gain emd phase margins in another loop. 

Therefore, to avoid this problem we shall define what we mean by multiloop 
gain and phase margins. This can be done with reference to Fig. 22 
where L(s) is cltosen to be a diagonal matrix. 

Definition; The multi loop gain margin is the pair of real numbers c^^ 
and defining the largest interval^ (c^^, c^) such that when It^(s) 

Hie could also use closed-intervals in the definition of these margins. 
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i “ in Fig. 22 are all real constants and satis^ the 

inequalities 

< 2-^ < Cj t i ■ l,2,...,m (68) 

the closed-loop system remains stable. 

Definition ! The multiloop phase margin is pair of real mnbers and 
-c. defining the largest interval (-c ,c.) such that when S.. ( jcj) « 
i “ l,2,...,m in Fig. 22 are of the form'*' e tdiere are real 

and satisfy the inequality 

-Cl < (j)i(a)) < Cl , i - l,2,...,m (69) 

and the closed- loop system remains stable. 

We will denote the multiloop gain margin of (68) by 

GM - (Ci, cj) <70) 

emd similiarly we denote the multiloop phase margin of (69) by 

PM - (-Cl* Cl) . (71) 

Note that in the SISO case (refer to Fig. 7) that these multiloop 
steibility margins reduce to the usual single stability margins but 
that in the MIMO case they differ from the stability margins obtainable 
a single loop at a time since these stability margins ai^ly in all loqps 
simultaneously . Of course, the word "simultaneously" does not mean that 
we can apply gain anA phase changes simultaneously in the same feedbacic 

_ j^TTu) 

we assume also that e has a stat^ space representation in order 

to ensure that S(s) has a state space representation. 
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lcK 9 but that only strict gain ohuigss or only strict i^iasss bk^^s n^f 
occur in ssparats fssdback diannsls sisultansously within ths prssoribsd 
li4Bits of ttM Bultiloop stabili^ margins. Ws oa^hasiss that thsas types 
of Bultiloop margins cc»isiter only a nail class of modslling «ncnrs ds- 
scribabls l/f a diagroal L(s) . With thsss prsliminariss ws arm ready to 
prasant tha following mrollarias to fhs«rna 4 and 6 rm^mctiwaly. 

Corollary It If has no CRHP saros and 

a . (l+<s“^(s)) > o (72) 

mn 

for all • d ^ than tiia multiloop gain and phase margins are bounded^ in 
the following manner 

CM D ll-o, l+ol (73) 

and 

PM D (-2sin"^ |, 2 ain" ^ | 1 . 

Proof » Fren Th^.ram 4 and (51) %#a know that no CRMP saros 

if for all s e 

(■)-!) < <’ 5 > 

and thus also if 


^rhe symbol 0 refers to sat inclusicm. Thus h D 1 moans tiiat I is 
contained in A. fhaa (73) means that tha upward gain margin is at 
least as big as l4o and that tha gain reduction margin is at least 
as nail as 1-a. Similar statements iqpply to (74) . 
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MUC “ 


( 76 ) 


for all 8 e r^. Now if L(s) is given by 


L(s) -diagtf^(s), £,^(s)I 


(77) 


then (76) implies that for all i 


|K,^(s)-l| < a 


(78) 


If f^(s) is real and is denoted by ti\en 


1-a £ £ 1+a 


j0, (w) 


(79) 


and if i^(ju)} is of Uie form e ^ with ^^(u) being real, then 


(u) 


- 1 < a 


(u)) 


( 0 )) 


- e 


< a 


4>. (tD) 

-a £ 2 sin £ a 


(80) 

(81) 

(82) 


or 


|4».(u))| < 2 sin"^ (|) 


(83) 


The bounds on the nultiloop stability margins follow from (79) and (83) 


Q.E.D 


Corollary 2 t If ''O CIQIP zeros and 


». i „ <>«<•' 1 ’ “ 


(Ml 
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for all • * ^ a £ 1 th«i tho analtiloop fain an>3 phaaa Mucfina ara 
boundad in tha folloainf aannar 

ife] 

and 

PM D 1-2 ain"^ I , 2 ain'^ |] . (86) 


Proof « Polloaing tha proof of Corollary 1 wo aiatilarly daduea that tiia 
Mrra^)Midin9 analog to ^nation (76) ia 

|l^^(a)-l|<a (87) 

and thua for raal i^(a) danotad wa aniat lutva that 


— ^ < i < — ^ 
l+a - i - 1-a 

j04(w) 

and for of tha fo» a , 6j^(b>) raal, wa hava 


( 88 ) 


-j^^(t»J) 


- i < o 


( 89 ) 


%fhich la^liaa 

|♦^(u)| < 2 ain"^(o/2) . 


( 90 ) 


Q.B.D. 


Xt Buat ba Mphaaiaad that corollariaa 1 and 2 provida worat 
caaa analyaia bounda on «ihat tha actual atability aargina ara. Ihia 
can ba illuatratad in tha SX80 caaa by Pig. 23 whara |l4g(a)| ^ a for 


all 8 € D so that Corolleury 2 is applicable 

R 


Possible Nyquist / 
diagram of qis)^ 


Fig. 23; Nyquist Diagram Illustrating Bounds of Corollary 2. 


Since the only information about the systan g(s) is contained in 
the single parameter a, the only information utilized by Corollary 2 is 
that g(s) touches the circle of radius a centered at -1 but that the 
Nyquist locus of g(s) never penetrates the interior of the circle. Thu 
to derive the worst case upward gain margin the corolleury assumes that 


H»e have used |l+g{s)| ^ ot rather than |l+g(s)| > a for convenience. The 
only modification of Corollary 2 is to m£Jce the bounds in (85) and (86) 
{^en sets. 
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g(s) passes thorough point A. Similarly , the worst ctise gain redia:ti<m 
margin wd worst case phase margin are (dstained by uswdng that g(s) 
passes through points B md C (or C) respectively, niese worst case 
margins are then useful bcnu^s <m the actual gain and phase margins. We 
refer to ttese bounds as guaranteed miniimm gain or plase margins . 


3. 7. 2 Crossfeed Tolerance 

The previous stability margins have asswed that I (s) is di^onal. 

If this is not true then there are cross couplings from one feedback 

channel to another as in the example ccmsidered in section 3.3. The 

ability to tolerate crossfeed type of perturbations is also determimd by 

the two quantities o (I+g”^{s)] and o II-»G(s)l as in the following two 

mxn min ^ 

corollaries to theorems 4 and 6 respectively. 

Corollajry 3 ; The polynosiial has no CRHP zeros and hence the 

pertvirbed feedback system is stable if the following conditions hold: 


1. condition 1 of Theorem 2 holds 


2 < 0„. (I-K3’^(s)) 

max min 


and 


,-l, 


a^^lY(s)] < a . (i-k; *(s)) 
max min 


for all s 6 and where L(s) is given by 


L(s) 


I X(s) 
Y(s) I 


(91) 
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Proof ; Iimediate from the form of L(s) in (81) and Ateoran 4. 

In this corollary L(s) of (91) represents a bilateral crossfeed 
perturbation %diere X(s) is the fraction of the control signals of the 
secmd group of feedback channels fed into the first gro\q> of feedback 
channels and Y(s) is tl^ fraction of the control signals of the first groi 9 > 
of fee<fi>ack channels fed into the sec<xid. If either X(s) or Y(s) is 
identically zero, then L(s) of (91) represents a unilateral crossfeed 
free one group of feedback channels to another. This is the particular 
form of crossfeed considered in Corollary 4. 


Corollary 4; The polyncmdal 4 >qj^{s) has no CRHP zeros and hence the 
perturbed feedback system is stable if the following conditions hold: 


1. condition 1 of theory 2 holds 


o lX(s)l < 0 . [I4C(S)1 
max min 


for all s e 0^ and where L(s) is given by 



— — 

I X(s) 


I 0 

L(s) = 


or 



0 I 


X(s) I 


(92) 


Proof ; Again immediate from Theoran 6 £md the form of L(s) in (92). 


3.8 Example of Section 3.3 Continued 

If was shown that the system of Fig. 11 under the feedback u » 
is nearly unstable if the value of b^^ very large. This neaumess to 
instability is easily detected using Theorems 3, 4, 5 or 6 because 


Magnitude (db) 


-IW)- 


a . II+g"^(s) 1 or O , [I+GCs)] becone very snail at frequencies below 
min min 

1 r«KS/sec. Fig. 24 s)»ws a plot of o^^^ [I-K»(ja)) 1 m a function of w 
with bj ^2 “ 50. 


-2.5 


-7.5 


-12.5 


-175 



-22.5 


-25' ' L. 

0.001 0.01 0.1 


Frequency (rad/sec) 


100 1000 


Fig. 24: Plot of o . (I+G(jo))) for Example 

mm 

of Section 3.3 (bj^2®50) , see 
Fig. 3.16. 


If we use Corollary 2 we obtain the following bounds on the 


multiloop gain and phase margins 
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GH 

0 

(.93, 1.08) 

(93) 

PM 

D 

(-4.1», 4.1«) 

(94) 


which are very conservative estiaates of the Bultiloop gain and phase 
margins. Nevertheless, they indicate a robustness probl» which is 
exhibited by the very small crossfeed tolerwce of Corollary 4 %ihich 
gives 

< O .^(I+G(jl)] s 0.071 = -23dB. (95) 

maX Bixn 

l^is again is a worst-case bound on the allowable amount of crossfeed 
at w=l but in this case it turns out that the magnitude of the error 
(i.e., E(s) = (2 ^(s)-G ^(s)lG(s)) induced by the crossfeed perturbation 
of Fig. 16 is -20dB, nearly the smallest necessary to destabilize the 
closed-loop system. 

3.9 Separating Functions and Additional Robustness Theorems 

At this point after having given several different robustness 
theorems, whose method of proof depended vpon the ability to ensure that 
I+G(s,6) was nonsingular on D x(0,l],we shall consider a OKsre general 
framework that allows us to generate stability tbeorens not necessarily 
derived from any particular error criterion as Theorems 4 to 6 were. 

After these additional theorems are generated we shall look for a possible 
associated natural definition of model error which if bounded in magnitude 
c£m not induce instability . 

In this section, we will define G(s,6) of Theorem 2 in terms of 


an L(s,e) giving 
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G(s,e) - G(s)L(s,e) (96) 

where L(s,e) is now continuous on D^x[0,lland such that (28) and (29) 
hold and find conditions on G(s) and L(s,6) to gueurantee closed-loop 
stability. Hec€dl from section 3.6, that the explanati«i of «diy Theoreu 
4 and 6 worlced is that they ensured (in the SISO case) that 

l+g(s,e) = l+g(s)£(s,e) 4 0 (97) 

on D x[0,l] . That is tt^y divided the congplex plane into t%fo disjoint 
R 

regions by using a circle and then ensured that -g ^(s) was in one region 
eind jl(s,€) in its camplement. The different theorems used different 
circles and thus give different allowable regions where I.(s,e) may be 
located. It thus seems natural to generate other theorems l>y choosing 
different circles to separate the values of il(s,€) «md -g ^(s). 

A simple way to specify a circle or line in the cosplex plane is 
to use a function f(-) known as a bilinear fractional transformation 
[52] given by 

f (z) = (98) 

cz+d 

where ad-bc ft 0 6tnd z is a ccxnplex variable as are a,b,c emd d. A 
circle or line can be specified by the equation 

|f(z)| * constant (99) 

where different values of a,b,c,d and the constant may give different 
lines or circles (refer to Fig. 25 for an example). The function f(>) 
has the property that it always maps circles and lines into circles and 


2z-(l-t-_1 ) 

z+1 


k . 


If(je(s,€))i<c,(s) 
region for iL(s,<) 
for c in [O, l] 
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lines. The Inverse fvnction of f(>) denoted f ^(.) also Is a bilinear 
fractional transformation and thus shares these properties. 

Now if we put -g ^(s) on the outside of a circle and ]l(s«e) on 
the Inside of a circle to separate th^ (Fig. 26) then we have a pair 
of inequalities of the form 

|f(-g"^(s))I > c^{8) > |f(£{s,e))| (100) 

on D^x[0,l] where c^(s) is a positive scalar. It may be that in (100) 

-g ^(s) is on the inside of the circle and i(s,e) on the outside de- 
pending on how f(«) is chosen but the key point is that f (•) in (1(H)) 
separates -g ^(s) and )l(s,e) and thus will be called a separating function . 

In order to develop a test that does not depend explicitly on € 
as in (100 ) , we may define JL(s,6) 

i(s,e) = f‘^[(i-e)f(i) + ef()l(s)] (loi) 

so that 


f(i(s,e)) = (i-e)f(i) +ef(Jl(s)) 


( 102 ) 


Now since (100) must hold for £(s,0) 1 and also for S.(s,l) * i(s) , 

(102) implies that 

lf(H(s,e))l < (i-e)|f(i)| +e|f(Jl(s))| (i03) 


or 

|f(£(s,e))l < max{|f(l)|, |f(i(s))l} 


( 104 ) 


and so we need only verify (100) at e«0 eind SbI. 


Now if 
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|f{>g"^(8))| > c^(8) > (105) 

«• are assured that (100) holds for all 6 in [0,1] since f(>) must be 
picked to separate -g ^(s) from 1 or g(s) from -1 and thus 

lf(-g"^(s))| > |f(l)C (106) 

Of course, the definition of l(a,e) in (101) may place restrictions on 
Ks) in order that Ha,e) be continuous on D x[0,l] that may need to be 
checked in addition to (105) . 

The preceeding discussion of the scalar case can be directly 

extended to the matrix case except that the circles becone hyperspheres 

and the cd>solute value signs in the inequalities are replaced by a . (•) 

nun 

or o (•). The objective now becomes to make sure that I+G(s)L(s,6) 
max 

is nonsingular or equivalently (assuming G ^(s) exists) that L(s,6) - 
[-G ^(s)] is nonsingular on D x[0,l] . 

Now suppose that we can find a function f(*) mapping to 

(t”*” such that f(A) - f(B) is nonsingular if and only if A-S is 
nonsingular for all A and B in for which f (A) and f (B) are defined. 

This means that the nonsinguleurity of L(s,6) - [-G ^(s)] can be checked 
in terms of the nonsingularity of f(L(s,6)) - f(-G ^(s)). A simple 
sufficient condition that guarantees the nonsingulaurity of f(L(s,6)) - 
f(-G ^(s)) is the singular value inequality (see (2.41)) given by 

o . [f(-G"^(s))l > o^_[f(L(s,e))l . (107) 

min max 

We again call f(0 a separating function since through (107) f(>) 
"separates" -G ^(s) and L(s,6) (i.e. , L(s,6) + G ^(s) is nonsingular). 
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In (107) I<(s,e} may be defined by 

L(s,e) - + 0f(L(e))] (108) 

which is amalogous to (101) . In (107) it is assmed again that 

0 , [f(-G‘^(s))] > o^^^tf(l)] (109) 

min max 

since with L(s,0) - I the nominal syst^ must satisfy (107) if f(«) is to 
be an af^ropriate separating function. Therefore (107) can t>e guaranteed 
for all e in [0,1] if 

for all s in 

In the matrix case, the separating functions f(*) may also be given 
by the matrix bilinear fractional, transformation 

f(X) * (AX+B) (CX+D)“^ (111) 

where A,B,C,D and X au:e complex matrices. To verify that they are indeed 
separating functions we present the following lemma. 

Learoa 1 ; If the matrices A, C3C+D and CT+D are nonsingular then X-Y is 
nonsingular if amd only if (AX-B) (CX+D) ^ - (AY+B) (CY+D) ^ is non- 
singular. 

Proof; Suppose X-Y is singular. Then there exists a vector z such that 
te ■ Y^ (112) 




aind thus 


-107 


{Mt+h)z - (AY+B)c (113) 

and 

(CX-H))e - (CY-»-D)£ . (114) 

Since CX4-D and CY-H) are nonsingular, let the nonxero ^ctor v be giv<ioi 
by 

V - (C3C40)* ■ (CY«»e (115) 

or 

z - (CX'«))‘^v - (CY+D)"^v . (116) 

Now substituting in (113) for z given by (116) we obtain 

(AX+B) (CX-rt»”^v - (AY+B) (CY+D)“^v (117) 

that is (AX4-B) (CX+D) ^ - (AY-fB) (CY+D) ^ is singular. To show the con- 
verse, assume that (117) holds for sme nonzero v and define a m>nzero 
z as in (116), then (113) holds and implies that 

A(X-Y)z • 0 . (118) 


Since A is nonsingular it must be that X-Y is singular. 

Q.E.D. 

One problem that occurs with the use of separating functims %d\ich 

are not defined over all of c”**”, as when CX4D is singular in (111) , is 

that the matrices X for which f(X) is not defined smst be examined for 

their effect on the continuity of L(s,8) defined by (108) as well as its 

effect on o . [f(-G ^(s))l in (110) which may alternatively defined as 
min 
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I |{f(-G~^(«))} ^ 1 ( 2 ^ (••• 2.39) %fhlch may ba well defined even thom^ 

0 . [f (-c'^Cs))] is not because (AX-t-B) in (111) may be invertible, 
nin 

Although in lensta 1 and (111) the A,B>C and D coefficients are 
matrices we will only use scalars a.b,c and d in the presentation of 
several additional theorems. The first of these theoms is a someidiat 
unusual version of the well-known small gain theorem [12] . This theorem 
is obtained by choosing f(X) ■ X, using origin centered circles. 

Theorem 7 (Small Gain Theorem) : The polynomial (s) has no CBHP xeros 

Cl* 

imd hence the perturbed feedback systm is stable if the following 
conditJ r ns hold : 


1. condition 1 Theorem 2 holds 


2. o;J^(-g"^(s))-|(g(8)||2 < 1 

for all s 8 




or equivalently. 


G(s) 




for all 8 e 




Proofs In this case Lis.€) is given by 


L(s,e) - (l-e)I * 6L(s) (119) 

and thus | |G(s)L(s,e) 1 1 ^ is sisqply botuided by use of conditions 2 and 3 as 

j |G(8)L(s,e)l I2 - lld-e)G(s) + eG(8)L(s)||2 (120) 

< (1-«)1|G(s)||2 + e||G(s)|l2l|l.(s)|| 


< (1-e) ♦ e - 1 
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ClMrly, fron (120), X4C(i)L(s,t) ■ X-H3(t,9) it noaaiiigalar and fcy 
TtiaoraB 2 $q^(*) lui* no CR8P aaros. 

Q.B.D. 

Savaral reaMrIca about thia thaorMn ara in ordar. Firat tba naaM 
8»all gain thaoran ariaae froa tha fact that cm^ititm 3 raqulraa tha 
loop gain to ba lass than tmity (saall anou^ not to dastabilisa tha 
closad-loop systan) . Furthanaora this version of tha thaoraai is rathar 
untMual in that typically tha conditions that hava iw CBHP sarcM 

(condition Ic of Thaoran 2) and that | |o(s) | I 2 < 1 are replaced by tha 
sinpla conditiOTs that G(a) and L(a) ara qpan-loop stable. Itote that 
||G(s )||2 < 1 and having ik> ClUIP aaros guarantees that G(s) is 

open-loop stable. Also L(s) need not ba stable as long as 0(s) and 
C(s) have tha sn»e nunbar of CRHP poles. Aecall fron section 3.4 the 
reason we require (*) have no CRHP zeros is that the noninal closed- 
loop systra nust be stable in order to detersiiiw its stability awrgins 
and determine if it is robustly stable, tie are not nerely detensining 
the stability of sorc arbitrary systtft with loop transfer sntrix G(s}L(s) 
where L(s) ■ I has no special significance. This the a^tin difference 
between robustness theorens and stability theorems . In robustness theory 
we are trying to determine when stability will be preserved and in stabi- 
lity theory' w..' are trying to determine \iseful conditions under which 
stability will occur without the benefit of knowing that with L(s) = X 
the feedback system (the nominal systm) is stable. Note also that 
condition 2 is sisq>ly randition 3 with L(s) ■ I and that condition 2 


I 


is the condition given in (109) . 
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Pinally wa point out tiuit G*^(a) naod not «xii»t ainM lay altar- 
noting fonnulating 1 ** ^ coi^plataly avoid tha 

^obl«n. Howavar, it is convaniant to parform 90 m formal mani^alatimts 
with thm f(-G~^(s)) in ordai' to gaiji insist on how to salact cartain 
usaful circlas and than go back and datarmina what assiaptions am 
actually nacasaary. 

It will now ba shown how all r^nistnass thaormra using tha S ■ GL 
toxm can ba understood as a mall gain tlMoram on an aquivalant faadbMk 
which is stabla only if tha original systM is stable. For this purpose 
wa introduce Fig. 27 where for TOnvmienca tha matrix L(s) ap|>aart in 
tha feedback loop instead of in series with G(s) and where %fa have 
st^ressed the depei^nce of G(s) and L(s) cm t. In Fig. 27-1 %fa have 
the original perturbed system which is transformed into Fig. 27-2 by 
use of a constant scalar multiplier a. Obviously, the systMS in Figs. 
27-1 and 27-2 are equivalent in terms of stability, that is, Mia is 
stable if and only if the other is stable. To go from tha system of 
Fig. 27-2 to that of 27-3 wa as^toy idut is known as a loop shifting 
transformation. This simply adds a pair of feedback loc^s with feed- 
back gains of tbX arouixl tha system G that cancel each other out 
because they have opposite polarity. Then cleverly, tha feedback 
loop around 4 O is moved so that it becomes a feedforward loop around 

Sk 

the «L system, again it is obvious tiiat the systesw in Figs. 27-2 and 
27-3 are equivalent in tenss of stability. Next, in order to go from 
Fig. 27-3 to 27-4 we define tha systems and shown by the dotted 
boxes in Fig. 27-3. Now we simply a^ply the sasM type of multiplier and 



Fig. 27: Loop Transformations with Multipliers Ilustrating 

the Relationship between the Small Gain Theorem 
and the use of the Bilinear Fractional Trans- 
format ion f ( • ) . 
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amd loop shifting transforaation in one step to the systea cos^sed of 
and to obtain Fig. 27-4. Now, however, we have in the loq^ shifting 
transforaaticm involving started out by adding two feedforward lo<9S 
of opposite polarity around cG^ and then aoved the tdl feedforward tera 
auraund the systea so that it becoaes a feedback loop, ^ain, we 

claim that the system in Fig. 27-4 is stability equivalent to the systan 
of Fig. 27-3 and thus stability equivalrnt to tl» original systen in 
Fig. 27-1. 

The next thing to notice is that the syst^s G 2 and L 2 defined by 
the dotted boxes in Fig. 27-4 cam be associated with a bilinear fractional 
tramsformation f(>) by the following equations 

G 2 - t-f(-G"^)l"^ - [(c-bd)G - adiHal + bG]"^ (121) 

amd 

L 2 = f(L) = [aL-bl] [adL - (c-bd)I]“^ . (122) 

Svq>pose now that we may prove the stability of the G 2 , syston of 
Fia. 27-4 by meams of the small gain theorem which has the basic in- 
equality 



This last condition, however, is equivalent to the condition 

< 124 ) 

mxn mauc 

where in both (123) amd (124) the dependence on s has been st;q>pres8ed 
amd must hold on D^. This shows that amy paurticulao: robustness test as 
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in (124) involving bilineeu: fractional transfomation may be formulated 
as a small gain test on an equivalent system. Also, thov^h it was not 
d(»ie here, the parameters a,b,c and d of the function f(.) can in 
general be stable minimun f^ase rational trsmsfer functions instead of 
constant scalars. 

^nfo final theorems which use a different separating function f(*) 
will be discussed. They are the well-known passivity theorem [12] and 
its generalization due to Barrett [47] . The passivity theorem we shall 
state has the same unorthodox assumption that closed-loop systmn is 
stcUale rather than the usual assvmption that tlte open-loop syst^ is 
stable. This again happens because we are using the theormn for deter- 
mining robustness of a nominal system under modelling errors rather 
then to ascertain the stability of an arbitreu:Y system. 

Theorem 8 (Passivity Theorem) ; The polynomial ^Q^fs) has no CRHP zeros 
and hence the pertxirbed feedback systma is stable if the following 
conditions hold: 

1. condition 1 of Theorem 2 holds 

2. G(s) + G** (s) ^ 0, s 6 

3. L(s) + l“(s) >0, s e Qjj 

Proof ; This proof is accomplished most simply without resorting to 
explicit use of separating functions and, therefore, they will not be 
used. Let I+G(s,6) by given by 


I-H3(s,e) = I+G(s) (d-€)I + 6L(s)] - l+G(s)L(s,e) 


(125) 
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and notice that L(s,€) is such that 

l”{ 8, e) + L(s,e) > 0 (126) 

on D x[0,l). Now suE^se that I-K3<s, 6) is singular for soi&e (s, 6) 
in D x[0,l]; then, there exists a nonzero vector x such that 
[l4G(s)L(s,e) ]x = 0 and hence 

x = -G(s)L(s,e)x . (127) 

Defining z = L(s,e)x, we note that z is nonzero else x in (127) is zero 
and thus 

z = -L(s,6)G(s)^ . (128) 

Condition 2 and (128) in5>ly that 

II 

z®G(s)z + z‘*g“(s)z = -z“G“(s)[L(s,e) + L(s,6) ]G(s)z > 0 

(129) 

cuid since G(s)z j* 0 a contradiction to condition 3 is obtained and thus 
I+G(s,e) is nonsingular on x[0,il. Theoron 2 again holds and the 
desired result follows. Q.E.D. 

Remark: In conditions 2 and 3 the strictness of the inequalities can 

be reversed and Theorem 8 still holds. 

Specializing to the SISO case illustrates the types of G(s) and L(s) 
that are required in Theorem 8. Conditions 2 auid 3 keeps g(s) and £(s) 
from entering the OLHP amd show (see Fig. 28) that since g(s) cannot 
encircle the -1 point, it must be open- loop stable in order to apply the 


Fig. 26: Admissable Region for g(s) and Us) 

in Theorem 8 (shelled) . 

theorem. It is fairly obvious that the phase of g(s)S-(s) is strictly 
less than 180** emd thus g(s)il(s) ^ -1. We can interpret conditions 
2 and 3 as separating -g ^(s) and S.(s) by the ju)-6ucis since -g ^(s) 
lies in the CLHP and £(s) in the ORHP. The ju)-axis can be viewed as 
a degenerate circle of infinite radius end we will use this notion to 
relate the passivity theorem to the next theorem which generalizes it. 

To derive the generalization of Theorem 8 we perform sane 
algebraic manipulations on condition 2 (and also condition 3) to relate 
these conditions to equivalent singular value conditions. First, note 
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that from conditio 2 we can deduce that 

2(G**+G) - (I+G)**{I-KS) - {1-G)“(1-G) >0 (130) 

and hence 

(I+G)“(I+G) > (1-G)“(1-G) . (131) 

This last inequality c£ui be rewritten as 

[(I+G) (I-G)"^l®t(I+G) (I-G)"h > I (132) 

or 

Oj|jinl(I+G) (I-G)'^l > 1 . (133) 

Similarly we can deduce from condition 3 that 

0 [(I-L) (I+L)"^] < 1 (134) 

max 

which when combined with (133) results in 

O . [(I4G) (I-G)‘h > 1 > O [ (I-L) (I+L) "^1 . (135) 

min — max 

In this last inequality, one wonders whether the 1 .^n the middle of (135) 
is really necessary and that if we use an inequality of the form 

0 .„l(I+G(s)) (I-G(s))"^ > 0^^^I(I-L(s)) (I+L(s))'^l (136) 

min max 

if it will guarantee closed-loop stability when (136) holds on The 

answer to this question is yes , provided that we impose sosie i4ditional 
restrictions on L(s) . The next theorem formalizes this. 

Theorem 9 (471 ; The polynomial (s) has no CRHP zeros and hence the per- 

' Cli 


turbed feedback system is stable if the following conditions hold: 
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1. condition 1 of Theorcn 2 holds 

2. L(s) has no real eigenvalvies less than or equal to -1 
for all s 6 

3. 0 j^^^I(I-k;( 8 )) (I-G(s))"^] > 0 ^,,[(I-L(s)) (H-L(s))“h 
for all s 6 

Remark ; If condition 3 is satisfied and (I+G(s) ) (I-6(s) ) 5. 1 

in condition 3 then it ceufi be easily slKuwn that condition 2 is auto- 
Buitically for all s 6 0^ for which this inequality holds. 

Proof; L(s,e) is given by 

L(s,e) « f‘^[e f(L(s))] (137) 

since f(I) =: 0 where 

f(X) » (I-X)(I+X)'^ - f”^(X) (138) 

and thus 

-1 -1 
L(s,e) = [I - e(I-L(s)) (I+L(s)) [14«(I-L(s)) (I+L(s)) 1 

(139) 

and is contimKJus on x[0,l] because of condition 2. Since f(0 is 

a separating function, condition 3 iiqplies that I-Ki(s)L(s,e) is 

nonsingular on D x(0,l) and hence 'Rieorein 2 holds, 
in 

Q.E.D. 

Note that in condition 3, the invert ibility of I-Ci(s) is not 
essential as long as G(s) ^ 1 for all s 8 since 
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a^^^Ul-*G(a)) {1-G{8))~h - II (I-G(«)){l4G(s))"^l|’^ (140) 

anA (I-fG(s)) ^ must exist because ^--(s) has no CRHP seros and thus 

wL 

1 1 (I-G(s)) (I+G(s) ) ^11^ is not zero unless G(s) ■ I. Utilizing condition 
3 of Theorem 9 we may derive 8K»ie corollaries on the stability SMrgins 
of the feedback system. 

Corollary 5 ; If ^q^(s) has no CRHP zeros and 

O , t(l4<3(8)) (I-G(s))"^l > a (141) 

fliixn 

for all s e and oi 1 then the multi loop gain and phase margins are 
bounded in the following manner 



and 

PM 3 [-2 tan”^a, 2 tan"^a] (143) 

Proof; Analogous to the proofs of corollaries 1 euid 2. 

Note that in the case ct>l we obtain t)^ boiuids on the multiloop gain 
eu)d phase margins associated with the passivity theorem which are given 
by 


GM D (0,<») 


(144) 


and 


PM D [-90», 90“] 


(145) 
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Also, as for Theoraos 4 and 6 we may derive a corollary involving the 
tolerance to crossfeed for systssui that satisfy condition 3 of Theorem 
9 but this will not be done here because it merely repeats the essential 
natures of corollauries 3 emd 4. 

The similarity between the theorems of section 3.5 and the theorems 

of this section is incomplete because it is not clear what type of 

modelling error is being bounded in Theorems 7,8 and 9. It hi^ppens that 

Theorem 7 cMnot be interpreted in terms of bounding the magnitude of any 

type of modelling error and Theorem 8 always bounds the magnitude of the 

model error by unity. This can be seen by identifying o „tf(L(s))l as 

the magnitude of the model error which is true in the case of Theorems 

4,6,8 and 9. In Theorem 7, o „[L(s)] 4 0 when L(s) ■ I and thus when 

in&x 

9(s) G(s) (no model error) the magnitude of the "error" (i.e., O [L(s)]) 

nax 

is not zero. Therefore a^^^ [L(8) ] does not correspond to the magnitude 
of a modelling error. Another manifestation of the lack of similarity 
between Theor^s 4 £md 6, and Theorems 7 and 8 is the fact that Theorems 
7 and 8 cannot be applied to all G(s) and S(s) that satisfy condition 1. 
of Theorem 2 whereas Theorems 4, 6 and 9 can. Theorems 7 and 8 place 
additional conditions on the allowed G(s) (i.e., in the SZSO case g(s) 
must lie inside of the unit disk in the complex plane for Theorem 7 
(Fig. 29) and the CRHP (Fig. 28) for Theorem 8). In Theorems 4, 6 £uid 
9 (again in the SISO case) the Nyquist diagram of g(s) may approach the 
critical point (-1,0) from any direction. This is not true for Theorems 


7 and 8. 
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Im 



Fig. 29: Adaissable Region for g(s) in Theoren 7 

To discover the underlying error criteria associated with nteorMM 
8 and 9 nudce the following identification between L(s) and B(s) given 
by 

E{s) - -f(L{s)) (146) 

since in ^eorens 4 and 6 it is f(-G ^(s)) + B(s) that is tested for 
singularity. Thus# in the case of Theorems 8 and 9 we have that 

E(8) - (L{s)+l)"^(I.{s)-l) (147) 

or since fi(s) - G(s)L(s) 

E(s) - {g“^(s)G(s)+1)"^(g"^(s)G{s)-D (148) 

and thus 

E(s) - (fi(s)4G(s)r^[fi(s)-G(s)] . (149) 

New note that we can write 2E(s) in the two following forms (drof^ing the 
s dependence) 


-1 
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(ISO) 


or 


2E 


(g+G) 

2 


IS-GI 


2B 


(- E*^ 
‘2 “l 




-1 


USD 


%^re 


and 


- g“^(8-g) 


Ej - -t5 "^ -g“^)g 


(152) 


(153) 


Fran (150) 2 E(b) can l>e interpreted as a relative error between fi and G 
vheta the base value is taJcen as the aritlumetic average of 8 and G. 
Another interpretation that is suggested by (151) is that 2E can be coe- 
pared to a resistance as can the errors E^^ and E^. Then (2E) ^ is li)ce 
a conductance that is merely the average of the c<Mvd\Ktances E^^^ and 
E^ . We note that E^^ is merely the usual relative error between G and G 
and that E^ is the negative of the relative error between S ^ and G 
In a sense the error criteria for Theorem 9 is a ^agproadse between the 
error criteria of Theoresis 4 and 6. Mote that as E^ (or E2) beccme 
small that 2E approaches (E ^) , that is 2E picks out the smaller of 
the two types of errors and uses that as a measure of the error. Fig. 

30 illustrates the nature of this error in a block diagram where by 
(147) L(s) is simply [l-E(s) ] (I-fE(s) ] 

This type of error criterion is pleasing in that it le^s to the 
symmetric (in a logarithmic scale) gain and phase margins of corollary 
5 and correlates well with classical single-loop sinultaiwous design 
requirements on gain and phase luirgin [47] . To put all the various 
theorems presented here in perspective, the Table 1 describes the 

















Fig. 30 1 Physical Representation of Perturbed Motel in nieorM 9 

sepsrating functions, error criteria and the nultiplicative L(s) factor 
corresponding to each the r^mstness test. 

In the next chapter, rteustness tests will be f emulated that 
utilize the str\x;ture of the modelling errors that were discussed in 
this chapter. This means that having an understanding of how errors 
enter into the system models will be important if any judgment about 
their structure is to be made. 

3.10 ExtMisions to Monl inear Systems 

The preceeding sections have dealt solely with tlie ^cabilit> /robustness 
properties of linear time-invariant systuts. The purpose of this sectim 
is to demonstrate that some of the theorems of the previous sections have 
corresponding nonlinear counterparts. These theorems SMy be proved by 
use of the well-known circle Theorem f'.O,!!] fonsulated by Zames and 
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later f^naraliaad by Savmwv [7, IS]. However « theee tbeoroM will not 
be proved hem doe to tiie lensthy dieoouim of oxtendM fWMtimi lene^ 
q>aces ai^ other neceaaary nathaMtical developwent reqolred. The key 
obaervation to be recognised is that the guaranteed gain nargin for 
TtMoresui 4« 6 and 9 raswin «cactly the easM when tM Multiplicative type 
of perturhatiem r^resented by the matrix L(s) is replaced by a nmlimar 
MBoryleas operator denoted as M (see Fig. 31) . This naans the gain of 
each feedback loq^ nay be changed as nml inear functi<m of tiie ou^ut 
signal of the plant provided the effective linear gain change is within 
the bounds specified by the guaranteed gain margin. This notatim of 
a gain margin for nonliiMar systems is made more precise in the Uteorsms 
of section 3.10.1. 

3.10.1 (hiaranteed Gain Margins for nonlinear Systests 

One of the first problems mcowiterjd in determining the stability 
of nonlinear systasw: is to clarify what is meant by the notion of 
stability. Various autlMrs define stability differently but the basic 
ewteept is tiut of boundedness. Thus, stability snist be defined before 
discussing the generalisations of Theorems 4, 6 and 9. For the purpose 
of this section we define stability in the following sumner. 

Definition (Stability) i h causal system with an arbitrary input u(t) 
aid correspondii^ output y^(t) is stable if there exists a nonnegatlve 
scalar k such that 


I|y<t)|l2 < k l|u(t)|l2 


(154) 
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%flMre I ! * 1 1 2 defined as 


I !*(t) 



T 

z 


(t)z(t)dt 


1/2 


(155) 


The nozn in (155) is proportional to the energy in the tine signal 
^(t) . 

Using this definition of stability^ we will examine the stability 
o£ the feedback syst^ shown in Fig. 3 where G is a linear time-in- 
variant convolutional operator representing the ncxainal loop operator 
and N(.) is a stemoryless nonlineau: operator given by 

N(x(t)) = tn (x (t)), n (x (t)) ,n (x (t))]^ (156) 

— XX ^ A lu lu 


where each n^(>) is a merooryless time -invariant nonlinearity and x^^(t) 
are the components of x (t) . 



Fig. 31: Nonlinear Syston 


In the completely linear time-invariamt case this definition of sted)ility 
requires that a stedile system have all its poles in the open-left-half- 
plcme. 
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ll) Fig. 31 %#e msubb that the nominal feedback ayst«B with ll(x(t)) > x(t) 
is stable am3 that the tremsfer function of the loop operator.G is given 
by G(s). Here M is playing the role of L(s) in the cosq^letely liiwar 
case (i.e., the perturbed loop curator 8 is given by GN) . 

A graph of the n^(x^(t)) conponents of N(x(t)) in (15C) sd^t be 
a saturation type nonlinearity shown in Fig. 32. 



In the next two simplified theorems, it is shown that by bounding the 
graph of the n^(x^) appropriately the stability of the closed- loop 
system is ensured and nonlineeu: guaranteed gain margins obtained. 

Theorem 10; The closed- loop system of Fig. 31 is stable if: 

1. it is stable with N(x) * x 

2. N(x) is memory less and time-invuiant and given by (156) 

3 for c « inf a . (I4G ^(ju))) and for all scalar x 
. «in 


(l-a)x < t\^(x) < (1-Hx)x 


for all i 
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Remark: Here it is asstsned for conveience that (ju>) ^ 0 for all (i). 

“ vis 

This theorem is the corresponding analogue to Theor«n 4 and gives 
bounds on the slope of the graph of n^(x) as riunm in Fig. 33. The 
bounds on a SISO case for the Nyquist locus of g(s) is shotm in Fig. 34. 
This is a simple application of the celebrated Circle Theorem [10,11] as are 
the next two Theorems. 


Theoran 11; The closed-loop system of Fig. 31 is stable if; 

1. it is stable with N(x) = x 

2. N(x) is a manoryless, time-invariant nonlinearity given 
by (156) 


3. for a = inf o . (I+G(i 0 ))) < 1 for all scalar x 
. - min — 

03>0 


1 , ... 1 

— X < n. (x) < ■; X 

1+CJ 1 1-a 


for all i 


Again we assume that (})„ ( ju>) 0 for all u) and observe that the gain 

OXi 

margin for Theorem 11 is the same as in the ccxnpletely linear case of 
Theorem 6. A similar nonlinear extension for Theorem 9 is averllable. 


Theorem 12; The closed-loop system of Fig. 31 is stable if; 


1. it is stable with N(x) ax 


2. N(x) is a memory less, time -invar iamt nonlinearity given by 
(156) 

3. for a * inf o . [(I-G(ju)) ^(I+G(ju)))] < 1 and for all scalar x 

U )>0 

X < n. (x) < X for all i. 

1-a 1 1-a 
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In these three theorms the quantity n^(x)/x is « 9 per and loifer 
bounded axx! can be considered as tl^ effective linear gain for the i~ 
leedback channel when the component x^(t) of the vector x(t) takes on 
the nunerical value x. Thus with this interpretation of n^(x)/x as 
an effective liiwar gain the guaranteed gain siargins for the nonlinear 
system are the same as those for the linear case. 

Notice that in the case of a saturation nonlinearity as in Fig. 

32 that Theorem 11 cannot be applied since a £ 1 in condition 3 Implies 
that n^(x) > l/2x which cannot be satisfied for a saturation nonlinearity. 

3. 11 Concluding Ranarks 

This section will attmnpt to give a perspective on the usefulness 
cind relationship of the robustness resxilts of this chapter. This chapter 
has presented a variety of robustness results and one wonders if there 
is a best robustness theorem to use in determining the largest class of 
model errors that the feedback system will tolerate. Practically > the 
newer to this question is no but theoretically Theorem 2 characterises 
the largest^ class of allowable perturbed loop transfer matrices {S(s)}, 
namely those whose multivariable Nyquist diagram is a deformed version 
of the multivariable Nyquist diagram for G(s) having the sasie number 
of encirclements of (0,0). However, the only practical way to detenairw 
if this is true is to use one of the rc^stness theoorcAS of sections 
3.5 and 3.9. These theoresis work with different types of SKxSel error 

^.argest under the restriction G(s) and G(s) have the saaw 
unstable poles. 


nusaber of 
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and one can only say that one theorem is better than another if a 

particular characterization of the model uncertainty has been selected 

to be the sense in which better is meant. 

For example, if one wanted to use the gain reduction margin as the 

criteria for the best theorem, that is, the best theorma would be the 

one that gave the smc\llest nianber for the gain reduction m 2 urgin upper 

bound. Then of Theorems 4, 6, and 9 one would say that Theory 4 is 

the best theorem to use since given any G(s), the upper bound on the 

gain reduction computed from o . ll+G ^(s)l is always less than or 

mxn 

equal to the upper bounds on the gain reduction margin conputed from 

0 . tl+G(s)l of Theorem 6 or frcan o . [(I-G(s)) ^(I+G(s))] of Theorem 9. 
min mm 

Simileurly, if one wanted the best indication of the upward gain 

margin, the lower bound computed from a . [I'fGCs)] of Theorem 6 is best. 

min 

These observations can be easily deduced via the relationships of a . [I'fGl 

m.n 

and o . [I+G of Fig. 21 emd similar relationships that may be derived 
min 

for 0 . t (I+G) (I-G) in relation to O . [I+G] or o . [I+G ^] . It 
min min min 

seems likely that in some sense that Theorem 9 should prove best but 
at present it is not clear what the particular criteria might be. 

Another way to compare Theorems 4, 6, emd 9 in the SISO case is to 
ccanpare the regions for allowable £.(s) given a r.ominal g(s). This is 
illustrated in Fig. 35 where gCjco^) * 3/4 for scnne As can be seen 
from Fig. 35, Theorem 9 places the least restriction on )l(ja)Q) in the 
sense that forbidden region for i(juiQ) in the complex pleme has the 
smallest area for any of the theor€ms. In general, these regions for 
)l(ju)Q) may overlap but may not be contained in each other, so that each 


Fig. 35: Allowable Regions for iUju^) if gCjw^) 

using Theorems 4, 6, and 9. 

Theoron 4: strictly inside circle A 

Theorem 6; strictly outside* circle 6 
Theorem 9: strictly outside* circle C 

*Except for 0 and the negative real axis 
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theorem may indicate a tolerance to a certain modelling error not 

guaranteed by either of the other two remaining theorems alone. Note 

that each of the circles in Pig. 35 pass through -4/3, the value for 

which makes l-fg(ja)^)£(jo]Q) - 0, and since cannot be on 

the circle's boundaries l+g( ju)^) £( jo)^) # 0 is ensured. 

In the MIMO case, drawing apprqpriate circles cannot easily be done 

and cong>arison of the theorems snist proceed by devising some ott^r 

appropriate criteria that is easy to check. 

The observations made so far, have been made on the basis of only 

the algebraic properties of the robustness inequalities of the theorems. 

However, using the typical frequency dependence of G(s) some additional 

typical cosq>arisons may be made. In order to obtain a good response to 

command inputs, typically of low freqtjmncy content, the loop gain in SISO 

systems is large in the frequency bamd where good following of the 

inputs is desired. The MIND generalization of the loop gain is given 

by 0 . (G(s)) and a (G(s)) where the former represents the lower 
min max 

bound on the loop gain of the "slowest" lo<^ of the feedback and the 

latter represents an upper bound on the loop gain "fastest" feedback 

loop. The crossover frequency of the SISO case becomes the frequency 

range where o . (G(s)} < 1 and a (G(s)) > 1 in the MIMO case, 
min — nmx — 

In the high performance (good command following) low frequency 

range o . (G(s)) is large and thus so is a . U+G{s)] (refer to Fig. 21) 
min min 

and in this region the tolerance to the modelling errors of Theoroas 5 
and 6 is generally good. 

In the frequency region above crossover o ^^^ (G(s) ) is small and 

thus a . [ItG ^(s)I is large emd, therefore, the tolerance to the modelling 
min 
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errora of Theorena 3 and 4 ia ganarally good. Thla la tha fraqutncy 
region where it ia inportant to have good tolerance to nodel error 
since in general nodelling error increaaea aa frequency doee (jecauae 
of unmodelled high frequency dynamics in ncminal design nodel) . This 
reccenends that Iheorem 4 always be used since Theorems 6 or 9 cannot 
be allied when the phase of the plant becoates coeipletely uncertain as 
it surely will at high enough frequency. A similar discussim of this 
nature is given in [43] . 

As mentioned previously in this chapter, the theoresmi using a analti- 

pllcative model perturbation or relative error Masure are generally 

fa\rored over the ones that use the additive model perturbatiion or 

absolute error measure, because the compensation mployed does not 

affect the nieAsure of modelling error. To make this clear, let <*p(e) 

denote the open- loop plant transfer matri.\ and G (s) the compensation 

c 

transfer matrix. Ttwn for the relative error criteria of Theoreaui 4 , 6 

and 9 with G(s) - G (s)G (s) we have that 

c p 

Theorem 4: G*^(s)(S(s) - G(s)J ■ G ^(s) [6 (s) - G (a)] (157) 

P P P 

Theorem 6: le'^(s) - g"^(s)]G(s) - ie"^{s) - c"^(a)lG (s) (158) 

P P P 

Theorem 9: (5(s) 4 G(s)r^ie{s) - G(s)) - rS^(s) 4- G^{s)r^lS.,(») 1 

P P P P 

(159) 

where 5 (a) is the perturbed open-loop plant model. Thus we see that 
P 

the compensation G^(s) does not affect the error computation. Ails is 
not true for the additive or absolute error oriteris. 
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4. Hostsmsis Muaysxs tor lxhbm sysibis with strucvurd 

NCffiSL ERROR 
4.1 Introduction 

Th* robuatnetfs tests of chapter 3 used only the magnitude of 
the model error in their formulation. It was shown there that if the 
model error magnitude is bounded by a MIMO generalisation of the 
"distance to the critical (-1|0) point" then the closed- loop stability 
of the perturbed fe^back system is guaranteed. Ho we ve r, there are 
many model errors whose magnitude is greater than the MIMO general- 
ization of "distance to the critical (-1,0) point" and yet the 
perturbed feedback system remains stable. 

In this chapter, the rc^stness tests of chapter 3 are refined 
to distinguish between those model errors which do not destabilize 
the feedback system and those that do, but both of which have 
magnitudes larger than the MIMO generalization of the "distance to 
the critical (-1,0) point". To do this it is necessary to be able 
to distinguish between model errors that increase the margin of 
stability for the feedback system and those that decrease it. litis 
cannot be done on the basis of the magnitude of the laodel error. 
Therefore, it must be done on the basis of the structure of the 
model error. 

The structure of model error, in general terms, is simply the 


ntmerical relationship of the elmnents of the error matrix E(s), 
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r^r«s«nting tin diffarwic* b«twMn th« noadMl Md tiM p«rtari»«d 
loop troAsfor aatrioos. Zti othor words, th« stnetturo of tho aodol 
is ^sclfisd by a^nltuds and phsM rslstimships bstwsMi ths 
s^^(s) slsBsnts of l(s). In this ohaptsr ths struetwrs of B(s) 
which is iapOTtsnt to dstonlM ths stsbili^ of ths p«rtorbod fsod- 
bsck systsBi is sxtrsctod using ths results of ehsptw 2 and tho 
singular value decomposition (SVD), to gei^ste an wthononsl basis 
for the eaepansion of E(s). Zt will be shown that the j^rojeotiMts of 
E(s) on only certain elMents of the basis need be known precisely 
to extract the infomation relevant for stability analysis. T1»s, 
only a partial character isatiwi of the aodelling error is iMccssary 
and its structure is constructively prodiMed by the iMthod of analysis 
used in chapter 2. Knotimx recently proposed uethod, principal gain 
and phase analysis 157] , %diich uses a somewhat different partial 
characterization of the s»del error to extend the robustness test in 
Theorsn 3.4, is discussed in chapter 6. 

In order to make a practical use of these results that utilize 
the structure of the model error, it is necessary to detezsdne if thm 
model error of minisnss iugnit\jde that will destabilize the f e ed bac k 
system can be guaranteed not to occur. This asses m s en t snist be made on 
the basis of engineering ju^esMnt about the type of model uncertainties 
that are reasonable for the nominal design model representing the phys- 
ical system. For discussimui on how to practically determine what 
constitutes a reaswiable modelling error, the reader is referred to 
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142] f«r A diacusalon of aodol orrora In an autoa»tiva MiglM 
control ayataai and 146] for a aimilar diacuaaion with ragard to 
powar ayat« aodala. 

SoBM knowladga of what la a raaaonabla aodal arror la abaelutaly 
aaaantial ainca all a»dala ara ut^artain In aao» fraquaney band. 

Modal arror al*qkya occura «dMn t}» fraquaney la aufflclantly high 
ai^ thla uncartainty auat ba aoomintad for. In Nmo control ayataaia, 
tha autxlinn croasovar fraquat^ whara tha lo^ tranafar aiatrlx# 

G(a), haa a nom of unity (l.a. tha naximai u for which 
I |G(ju) I Ij"!) Must occur in a frequency band where the aiodal atlll 
adequately rapreaanta the phyalcal ayatan If atablllty la to ba 
enaurad. It la to tha dealgner to decide how and In what my 
the nodal la \incertaln. 

Having now briefly daacribad the key role of tlM aaidal arrmr 
atructure for tha raaulta of thla chapter It la appropriate to 
outline the ramlnlng aactlonaof this chapter. In aactlon 4.2, It 
la shown exactly how tha structure of the aodellng error can ba 
used to obtained iaprovad versions of the theoroas in chapter 3. 

Theoreais 4.1 ami 4.2 show that tha necessary auignltuda of tha 
nodal error, at a particular frequency, that destabilizes tha feed- 
back ayatan, but is essentially unlike In structure to the nallast 
possible destabilizing nodal arror, nay ba such larger In aagnitude 
than the mgnltude of the snallast destabilizing a»del arror. 

This mans by differentiating tha nodal errors on the basis of their 
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•imllarity to tho Btnwtttro of tho HMllMt doatabillBliii orror* 
th« foodbock oystM eon guarantood to tolarata a pcMaibly moh 
largar Bndal arror. 

Thia is aaplainad flrat for tha 8X80 eaaa and than g«Mraliaad 
to ttw MXIIO casa. 8aetiwi 4.3 Intarj^ata tiia natura of tha aodal* 
ling arror of alniwa auignituda, that daatablllaaa tha faadbadt, 
via block dii^an nanipulatiMa. Maxt in aactlon 4.4, tha aarmpla 
of ch^tar 3 (illuatrating tha daficlanolaa of tha alngla*loop typo 
of atablllty auurgina} la oontinuad to ahow that tha analyaia of thla 
chaptar pradicta Um typa of a»dal parturbatiMi uaad to d aa on atratad 
tha naar Inatabillty of tha eloaad-lo^ ayatan. Finally, in aacticm 
4.5, tha poaaibillty of coalbinlng diffarant rohuatnaaa taata aa a 
way of aactanding Uwir uaafulnaaa ia diacuaaad. 

4.2 Robuatneaa Taata Utiliaii^ Modal Error Structnra 

In the rolmatnaaa thaorana of chaptar 3, tha key conditiona 
anauring the atabllity of the perturbed cloaed*lo<v ayatan ware 
inequalitiea of the fon 

a tE(a)l < o . Ih(C(a))l (1) 

nax Bin 

tdtare h(«) ia aone bilinear fractiMial trAnafomation 
(i.e., 1+0, I+o’^, (I-G)“^(l+<!)) and where tlj nu3c hold for all 

ia tha portion of in Fig. 3.10 
for which Ra(a)<0, Thia ctmdition aaauraa that the nodal arrwr 


Recall froB (3.33) that 
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Is sufficiently snail so that a closed-loc^ syatan designed on the 
basis of G(s) will r«iain stable tdian It is replaced by G(s). 

However, the approach used to develop these rolnistness theorems 
neglects the fact that there are perturbations or modelling errors 
for which (1) does not hold, i.e., the model error is not small, 
and yet the closed-loop system remains steUdle. These chapter 3 
theorem are conservative if one restricts the allowable type of 
model error structure because they guard against aUssolutely all types 
of structure in linear iK>del errors. 

One way to reduce this conservatism is to obtain additional 
conditions that distinguish between modelling errors that do not 
destabilize the feedback system but violate the test of (1) , and 
those that violate the test of (1) but also destabilize the feedback 
system. Or better yet, obtain some conditions that discriminate 
between modelliiig errors, that violate (1), between those, that 
increase and those that decrease t:he margin of stability of the feed- 
back system. 

The problem is illustrated in Fig. 1 for SISO systems where two 
different perturbed systems g^(s) and ^ 2 ^®^ produce exactly the same 
size of relative error on the Nyquist diagram. As cam be seen from 
Fig. 1, the difference between the perturbed systons g^^is) and g^Cs) 
cannot be detenained from the meignitude of the error alone. 

Clearly, g^(s) has a smaller margin of stability than the nominal 




Fig. 1: Tho Different Perturbed Models with the same 

Relative Error Magnitude on a SISO Nyquist 
Diagram. 

system g(s), and g^(s) has a larger margin of stability than the no- 
minal g(s). Since this is a scalar system the only additional 
information about the error needed to distinguish between g^(s) and 
^(s) is the phase of the error. Thus, in the SISO case this gives 
us a con^lete characterization of the error. 
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At point A, in Fig. 2, th« Nyquist locus of g(s) is nssrsst ths 
critical -1 point and thus the worst error siiaply noves point A to 
A* by "stretching" the Nyquist locus at that particular frequency to 
just pick up an extra encirclement of the -1 point (the point A* Is 
infinitesimally close to -1) . It is ijg)ortant to point out that this 
type of perturbation could be <g>plied to q(s) in any frequency ran«» 
but that it need happen only at one p^u^ticul^u^ frequency , (i)q near A » 
in order to induce instztbility. Thus we will speak of the worst 
error at a paurticular value of 

Notice also that there are auiy nisoaber of curves that we could 
pass through A' representing perturbations of the original Nyquist 
diagram of g(s) as depicted by g^^Cs) in Fig. 2, that induce instability 
and are identical to the worst error at the frequency of point A but 
differ at other frequencies. However, these curves will also be 
considered as worst errors since it is really their nature at a 
single frequency that is important in distinguishing them from other 
curves. 

One other point must be emphasized. The system g(s) may be 
constructed quite siJiply by finding a continuous stadale 
Us) = g(s)/g(s) that meets as closely as desired the ideal speci- 
fications given by 

(2) 

i » 
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where is the frequency corresponding to point A in Fig. 2. For 
exeonple, one continuous, stable £(s) that ap]^oxlinates A.. ,(s) in 

idMl 

(2) c«u\ be generated simply by taking £(s) to be of the fon> 


where 


Kb) - 1 - q(s) 1 + g ^(j^Q) 


q(s) " 


2 £ 


s +2P0 )qS+(u 


2 

0 



(3) 


(4) 


To approximate ^ closely, p>0 in (4) must be very small so 

that |q(s) I is as small as desired whenever |s-jo)^|>e for a given 
e. The constants a^O and c«+l in (4) are used to adjust the phase 
of q(s) without affecting |q(s) | so that 

qtjw^) = exp[j{arg(l+g ^(jw^)}] . (5) 

This selection of p, a and c in (4) makes q(s) essentially zero 
everywhere except in a suitably small frequency range near 0 )^ where 
it has the value given in (5) . Thus Ks) is as close as desired to 
the specifications in (2) but is still continuous in s and stable. 

The l{si determined by (3) , (4) ^md (5) produces a g(s) essentially 
like the one of Fig. 2. 

Returning to the MIMO case, we cem m^Jce all the euialogous state- 
ments to those concerning Fig. 2, once we have specified the worst 
error. Then similarities between the SISO Md MI>W cases cem be 
easily demonstrated using the ideas of chapter 2 developed in 


-144- 


Problems A euld B and by use of the SVD on the matrix h(G(s)) of (1) . 
Suppose that the SVD of h(G(s)) is given by 

h(G(s)) - U{s)£(8)v”(s) (6) 

where 


U(s) = lu (s) ,u (s) , . . . ,u (s)] (7) 

— i — z ~n 

V(s) = tv (s) ,v (s) , . . . ,v (s)J (8) 

—1 — z ~n 

I(s) = diagto, (s) ,a_(s) , . . . ,a (s)] (9) 

L z n 

a. (s)> a. (s)> 0 (10) 

1 — 1+1 — 


where the singular values O, (s) = a (s) and a (s) * o , (s) . 

X indx n roxn 

Recall from (2.66) that the error matrix E(s) of smallest norm that 
will make h(G(s)) + E(s) singular is given by 


E(s) = U(s) 


... 

0 

0^ 

-0 (s) 


n 


H 

V(s) 


( 11 ) 


where |!e„(s)!|< o (s) but is otherwise arbitrary. Provided the 
norm of the matrix Eq(s) is bounded by dj^ts) , its structure is 
con^letely unimportant information for the test determining the 
singularity or nonsingularity of h(G(s)) + Eq(s) • Therefore, Eq(s) 
will be taken as identically zero in the following discussion and 


X ^ 

of course it must also be such that G(s) satisfies condition 1 of 
Theorem 3.2, 
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thus, E(s) giv«n by (11) reduces to 

E(s) = -o (s)u (s) 2 ^{s) . (12) 

n — ti ~n 

The E(s) given by (12) will be called the essential structure of 

the more general form of E(s) given by (11) when E^(s)»<0. The 

H 

quantity -O (s)u (s)v (s) is the component of E(s) given by (11) 
n ~n ~n 

that alone must be exactly Icnown if it is to be ascertained whether 
or not the matrix h(G(s)) -t- E(s) is singular. Hence, the description 
of the E(s) given by (12) as the essential structure of E(s) given 
by (11) is justified. 

Remark ; The fact that E(s) in (12) is singular or that E(s) in (11) 
may be almost singular will be important in chapter 6 where a method 
that assumes that E(s) is nonsingula- or not even close to 
being singular is discussed. 

Again, as in the SISO case, the error given by (12) need only 
occur at one particular complex frequency s^ to destabilize the 
feedback system. That is, we may construct a perturbed 8(s) having 
the same number of unstable poles as the nominal G(s) that has the 
property that E(Sq) satisfies (11) arbitrarily closely and hence 
destcd3ilizes the feedback system. The MIMO error matrix 
E(s^) = -0 (s„)u (s^)v**(s„) is the generalization of the model 
errors that produce the g(s) emd gj^(s) of Fig. 2 passing through point 
A' just picking up an extra encirclement of the critical point (-1,0). 
From (12) we see that for an arbitrary error matrix E(s) that the 
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X H H 

projection, <u (s)v (s),E(s}> u (s)v (s) , of E(s) onto the 
— n “Ti -n — n 

H 

one dimensional subspace spanned by u (s)v (s) can be used to 

n — n 

determine if the component of modelling error in the roost sensitive 
direction u (s)v (s) will move the multivariable Nyguist diagram 

—n — n 

of the ncaninal system nearer or farther from the critical point 
(0,0) in the complex plane. The direction of this movement of the 
MIMO Nyguist diagram is simply ascertained by determining if 

<u (s)v (s),E(s)> is nearer or farther than a distance of 0 (s) 

—n ~n n 

from the point (-o^(s) ,0) in the complex plane. However, the quantity 

<u (s)v (s),E(s)> merely determines the effect of one component of 
— n -n 

the model error and does not ta)ce into account the effect of the 

other components of the model error (i.e. , the projections 
H H 

<u , (s) (s) ,E (s) > u.(s)v.(s)) have on the multivariedsle Nyquist 

-1 -1 - 1-3 

diagram. Therefore, some restrictions on these other model error 

components must be placed if their effect on the stadaility of the 

closed-loop system is to be easily predicted. 

SupiX)se now that we restrict the component of modelling error in 

the most sensitive or worst direction u (s)v (s) to be exactly zero 

— n — n 

(i.e., <u (s)v (s),E(s)>=0) so that it has no effect on the 
-n -n 

multivariable Nyquist diagram. Naturally, for this class of modelling 
errors, one expects that the magnitude of the error required to 
destabilize the feedback system should increase since the worst possible 


1 


The innerproduct notation <•,•> was defined in (2.40) of chapter 2 
where a discussion of projections on subspaces is also given. 


-147- 


t>’pe of error has been ruled out and indeed this la the case. The 
elimination of this type of error can only be done using ei^lneering 
judgement aJaout what type of error can occur in the physical systmn. 
The next theorem assiimes that the worst model error can be ruled out 
and extends Theormns 3.4, 3.6 and 3.9, by allowing them to deal with 
errors of larger magnitudes than previously allowable. 


Theorem 1; The polynomial 4 (s) has no CRHP zeros and hence the 

' ■ Cli 

perturbed feedback system is sted^le if the following four conditions 
hold: 


(a) (s) and (s) have the same number of 

OL 


CRHP zeros, 

(b) if 


or 


(c) 4>_-(s) has no CRHP zeros 

CLi 

h{G(s)) is of the form: 


(a) h(G(s)) = I+G(s), A(L(s) )0(-«,O] and 
E(s) = Ig'^(s)-g'^(s)]G(s) or 

E(s) = G(s)-G(s) for all 

(b) h(G(s)) = (I+G(s)) (I-G(s))“^, X(L(s))^(-<»,-l] 
and E(s) = [G (s) +G (s) ] “^ IG (s) -G (s) ] for all 


or 


(c) h(G(s)) = I+G~^(s) and E(s) = g”^ ( s) IG (s) -G (s) ] or 

E(s) - [g'^(s)-G~^{s)] and X (L(s) for all 

sen . 

R 
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3- (8)a ,(8)J^^^ 

iBax n n**l 

for all sCflL where o (s) and a , (s) are the two 
R n n-1 

siaallest singular values (assuned to be distinct) 
of h(G(s)) 

4. <u (s)v*^(s) ,E(s)>"0 

— n — n 

for all s8&^ where u (s) and v (s) are the left and 
R -n -n 

right singular vectors of h(G(s)) associated with 

a (s) . 
n 

Proof ; Conditions 1 anu 2 are the same conditions used in Theorems 
3.3 to 3.6 and 3.9 to ensure that G(s,€) is continuous in 6 on 
D x[0,l] so that Theorem 3.2 can be applied. Therefore, we need only 

K 

show that h(G(s)) + E(s) is nonsingular. This, however, is 
guaranteed by conditions 3 and 4 using the solution to Problem B in 
chapter 2 (see (2.73) to (2.76)). Q.E.D. 

Note that in Theorem 1, conditions 3 and 4 are required to hold for 
all even though they need only be used in the frequency range 

where the sufficient conditions (all given by (1) of this chapter) of 
Theorems 3.3 to 3.6 and 3.9 are violated. 

The significance of Theorem 1 is that by requiring very little 
information (condition 4) in addition to the magnitude of the model 
error, the worst type of modelling error that could destabilize the 
feedbac)c system (and whose delusion might be justified on physical 
grounds) is effectively eliminated. Hence, the "size" of the error 
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neccsiary to destabilize the systsn nay increawe significantly if 

0 ,(s)»o (s). Thus, the conservatisn of the chapter 3 theoresis 

n-i n 

for this class of nwdelling errors is reduced. The essential atructure 
of the next worst error (i.e., next snallest error) that destabilizes 
the system in this restricted class of modelling error.; is given by 
(from (2.73) with because <u (s)v (s) ,E(s)>«0) 


E(s) \/a ls)a , (s) f u (s)v** , (s)e^® 
n n-1 -n -n-1 


+u , (s)v (s)e 
~n-i ~n 


•j9(s)j 


(13) 


where (a) 9(s) is real and arbitrary and (b) the vectors 

u . (s) ,u (s) , V , (s) and v (s) are the left and right singular 
~n-l — n — n-1 — n 

vectors of h(G(s)) corresponding to o , (s) and o (s) respectively. 

n-1 n 

The spectral norm of the matrix E(s) in (13) is precisely 

Ja is) a . (s) . 

^ n n-1 

However, it must be pointed out, that it is extremely unlikely that 

condition 4 of Theorem 1 will hold exactly for a realistic modelling 

error since the model error in the i>articular direction u (s)v (s) 

— n -n 

will rarely be exactly zero. A more likely expectation is that this 

con^jonent of the error not be exactly zero but sufficiently small in 

magnitude. By requiring only that the model error in the direction 

u (s)v (s) be sufficiently small. Theorem 1 may be modified so that 
~n — n 

the essential nature of its results are still valid when the class of 
model errors considered is characterized by 

|<u (s)v^(s) ,E(s)>| < c(s)< a (s) “ o , (s). 

' -n -n ' — n min 


( 14 ) 
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The positive scalar c(s) in (14) bounds the nagnitude of the worst 

nodelling error as a fimction of frequency to be less than O . (s). 

ain 

the aininuBi nagnitude of the anallest destabilizing error required 
to destabilize the feed^c)c syst«<i. Therefore, the magnitude of the 

II 

model error in the most sensitive or vforst direction u (s)v (s) 

-n -n 

is not large enough by itself to destabilize the feedback system. 

In order to destabilize the feedback system when the model errors 
satisfy (14) , other model error components, besides the model error 
component in the worst direction, must contribute to the movement 
of the MIMO Nyquist diagram through the critical point (0,0). 

This is stated formally in the next theorem. 

Theorem 2 : The polynomial no CRHP zeros and hence the 

perturbed feedback system is stable if the following conditions 
hold: 

1. conditions 1 and 2 of Theorem 1 ho.^d 

2. a lE(s)]<lO (s)a ,(s)+c(s)IO (s)-0 ,(3)]]^^^ 

max n n-1 n n-1 

for all 

3. |<u (s)v*^(s) ,E(s)>|< c(s)< O (s) 

— n —n “ n 

for all . 

Proof ; Identical to proof of Theorem 1 except that now the general 
solution of Problem B ((2.73) to (2.76)) via conditions 2 and 3 
guarantees that h(G(s)) * E(s) is nonsingular. Q.E.D. 
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Th« ••••ntial structure of the next %foret perturbetion that 
does not violate coi^ition 3 but destabilises the feedback system 
is 9 iven by (from 2.73) 


E(s) 

where 


Ic(s)u ,(s)v‘* (s)-c(s)u {s)v**(s)4y(s)u ,(s)2r”(s)4Y*(s)u (s)v** (8)1 

Ti— 1 “Ti— i “H “n— 1 ~n n“A 

(15) 

Y(s) - [lc(s) + o^(s)nc(s)-0^_^(snj^^V*‘*^ (16) 


with ^(s) being arbitrary but real. Note that as c(s)-*' 0, in con- 
dition 3 and in (15) and (16), that we recover the results of Theorem 1. 
To make the meaning of the results of Theorem 2 clearer, the following 
example is given. 


Example 1; Suppose that we wish to determine stability robustness of 
a 2x2 control system which actually has a loop transfer function 
matrix G(s) but is represented by the nominal diagonal lo<^ transfer 
matrix G(s) given by 


- 


r 1 

•» 



s+7.5 

0 


- 



0 

K) 

M 

01 

1 


0 

1 

8+0.5 J 


(17) 


so that the nominal closed- loop system has poles at -6.5 and -1.5. 


If we use the relative error criterion 


E(s) - g'^^Cs) {G(s)-G(s)] 


’ll"’ 



(■) 

(s) 


’ll"’ 

’ll"’ 

’ll "’-’ll"’ 


(lei 
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than tha ■ultipllcAtiva uncartainty factor matrix L(s) la 9 ivan by 


L(s; - 1+E(a) - 






»n<*’ 

’ll'*' 




(19) 


First, we coi^ute CF . (T+G ^(jw)) to detemlne the aagnituda of 

XRin 

the smallest destabilizing model error E(s). This is slowly given 
by 


o^indtG (ju») 


1.5+ju 


'^(1.5)^-Hu^ > 1.5 


( 20 ) 


I 


because 


I+C'^(s) 


s+8.5 0 

0 s+1.5 


(21) 


Now suppose that the error in the loop gain of each loop of 
the feedback system is known within +50% of the nominal loop gain, 
that is 


0.5 < 


and 


0.5 < 


g^^(ju)) 


g„(jw) 




< 1.5 


( 22 ) 


“ 1^22^^“^ I 1 • 


(23) 


Next, suppose that we are more uncertain about the channel crossfaeds 
in tha sense that we can only assert that 




(24) 


and that 


|e2i(ju))|= 


g^j^(jw) 


< 2 


92i<3“> 
’ll ”“> 


< 2 , 


(25) 


It follows fr<MO (22) and (23) that we can bound jej^j^(jw) j aiid je 22 {jw)| 
by 1/2 and thus, by (24) amd (25), we can only conclude that 


1 lE(jw) 1 L = a ^lE(joj)J< 2.5 . (26) 

From (26) and (20) it is clearly possible to have 

0 lE(ja))]> 0 . Il+G'^(jw)] . (27) 

max min 


Therefore, Theorem 3.4 of chapter 3 does not apply. However, we can 


use Theorem 2 to ensure the stadiility of the perturbed feedback system. 


To see this, note that the SVD of I+G ^(ju) is given by 

0 lfljaw-8.5 ! 0 


I+G"^(jw) ^ 


je^(uj) 




0 


jw+1.5 


1 0 


0 1 


= U(jw)E(jw)V (ju) 


(28) 


where 


and 


9^((jj) = arg[j(i>+8.5] 
t\(u) = arg(jLH-1.5J . 


(29) 

(30) 
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Note that condition 3 of Theoren 2 can be satisfied with c(jo))*l/2 
since from (28) defining U 2 (j“) and from (23) bounding 

?^^22^3w) and thus thit for all b) 


l<U2(jw)v”(ju)),E(j(i))>j»|u^(j(D)E(jui)v2(ja)) l = |e22(ju) If. 1/2 • (31) 

Thus, by (31) and (20) we have 

a2(jw)> 1.5 > 1/2 > l<U2(jw)v^^(ju)),E(ja))>|. (32) 

Next, we calculate the right-hand- side of condition 2 of Theor«n 2 
and a lower bound as follows 


(jw)02 ( jb))+c(jas) lo., { ju)- 0 j^ (jo)) ]j = 1^1 ju)+8.5| I j{iH-1.5| 

+ 1/2(1 joH-l.sl-l ja>+8.5|]J^'"^ > (8.5) (1.5)+(“|) > 3. . (33) 

Therefore, using (26) we have that 


a lE(jw)]< 2.5 < 3 
mauc — 


<[c, ( ;o))o.,(ja))+c(j(i)) lo (ju))-a 



1/2 

(34) 


and so condition 2 of Theorem 2 holds. Assuming condition 1 of 
Theorem 2 holds we liave shown that the perturbed feedl3ac)( syston is 
stal-le. The next smallest de.stabilizing error can be calculated 
from (15) and (16) with ^(ja))»0 and u)=0 


r 

I 


since O'nijjfl+G "^(O)) 


L E(0) 


1/2 3 

3 -1/2 


1.5 and Is given by 


(35) 


which means that L(s) may be taken as the constant matrix L given 
by 



3 

1/2 


(36) 


Thus, we see that (refer to Figs. 3 and 4) crossfeed gain errors of 
magnitude 3 and loop gain changr-^ of jt50% are required to dest2d3ilize 
the feedback system if we insist t (22) and (23) must hold. 


Nominal Open -Loop System G(s) 
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the matrices U(a) and V(s) of the SVD of h(G(a)) axe cooplex it is 
very unlikely that L(s) determined by the E(s) given in (13) or (15) 
will even have real eigenvalues. 

We can now consider placing additional constraints on the 
modelling and further res':rict the class of allowable modelling errors 
in the manner of I’roblan C in chapter 2 and derive the next theorem. 

Theorem 3 ; The polyncxnial zeros and hence the 

perturbed feedback system is stable if the following conditions hold; 

1. Conditions 1 and 2 of Theorem 1 hold. 

2. E(s) is of the form 

H 

V(s) (37) 


where ^^(s) and e^(s) are vectors whose last 
component is identically zero and where U(s) and 
V(s) are defined in (6). 

3. 0 (E(S) )< v/o. (s)0. (3) 

mcoc ^ 

where o (s)On(s) = min a. (s)o (s) (38) 

" (l.JXW * ’ 

A 

M= { (n,n) , (n-l,n) Xn,n-1) } 


E(s) = U(s) 


Ej^(s) 




^3 


and 


(39) 
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Proof: Direct application of Problem C of chapter 2 and Theorem 3.2 

as in Theorems 1 and 2. 

Theorem 3 allows us to determine the next larger magnitude of 
the "next, next worst model error" required to produce instability 
when the smallest destabilizing model error and the next smallest 
destabilizing model error considered in Theorem 1 and given by (13) 
are con^letely eliminated from consideration. Theorem 3 eliminates 
these type of errors by requiring zero model error projections in 

II 

the worst direction u (s)v (s) and the next worst pair of directions 

— n — n 

H H 

u (s)v , (s) and u , (s)v (s) . The process of eliminating each 
-n -n-1 -n-1 -n 

"successively worst direction" could obviously be continued and 
larger magnitudes of these classes of errors would then be necessary 
to destabilize the feedback system. 

4.3 Block Diagram Interpretations of Worst Model Error 

In this section, interpretations of the smallest destabilizing 
model error will be given using block diagrams revealing the role of 
the SVD of the matrices of I+G(jo)) and I+G ^(ja)) in the input-output 
properties of the feedback system. The types of model error con- 
sidered are those of Theorems 3,4 and 3.6 involving the relative errors 

'• -1 --I 

between G{jw) and G(ju)) or G (joi) and G (jo»). 
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At 8cae particular frequency W , let the SVD of I-KS” (j« ) 


be given by 


I-K3~^(jW ) => U(jW )I(jW )V^(jW ) 


( 40 ) 


SO that the closed- loop transfer matrix at u , G (jw_) is given by 

U U 


“ (i+G~^(jUo)'^= v(jUQ)E"^{jajQ)u”{ja)Q) 


( 41 ) 


and thus 


H 




ai(jwo) 


( 42 ) 


A block diagram of a closed- loop steUcle system representing equation 
(42) is given in Fig. 5. 



Fig. 5; Block Diagram Interpretation of SVD of 

I+G‘^(jO)Q). 
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This Flgiire Illustrates that the left and right singular vectors 

u^(ju>Q) and collapse the MIMO closed- loop syston into a SISO 

system through which the signal cosw^t passes with a change of 

amplitude by a factor of l/o^ ( determined by the singular values. 

These vectors can be interpreted as input and output "directions" 

where for each different value of the index j input/output direction 

pair produces a different SISO system and represents a different route 

through the MIMO system for the signal cosw^t. Therefore, from 

Fig. 5 with j=n, we see that if 0 (ju„) =• 0 . (jw^) is near zero, then 

n 0 min u 

the system will an^lify a sinusoidal signal by a large factor of 

1/0 ( jU3 ) in tlie input/output directions of u ( jw ) and v ( • As 
n 0 — n 0 ~n u 

0 (iu)„) approaches zero the amplification factor approaches infinity 
n 0 

until at 0 ( jw „)=0 the system with a bounded input produces an 
n 0 

unbounded output, that is, the system becomes unstable. This is all 

rather obvious since if 0 (ju)„)=0, the matrix (1+G ^(jOJ.)) ^ does not 

no 0 

exist and therefore there must be closed- loop poles on the jw-axis at 

±l“o- 

As in the case of I-t-G a similar interpretation of the 

SVD of I+GCjuJ^) can be made. If the SVD of I+G(ju)Q) is given by 


I+G(ju^) = U(jW- >r (ju>-)V( ju.) making the SVD of lI+G(jio.)] 
0 0 0 0 0 


-1 


-1 -1 H n 

[l+G(jWQ)]’ = (jUQ)u(j(jjQ)= I 


i=l i 


0,(ju,jVi(3CUQ)u^(jV (43) 
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then Fig. 6 gives the block diagram analogous to Fig. 5. Figure 6 
shows that the only change from the previous case shown in Fig. 5 
is that the output is generated from the error signal e instead of 
the system output signal y. 


COS 



1 


COS 


Fig. 6: Block Diagram Interpretation of SVD 

of I+G(jW|^). 

Notice that in (^1) and (43) the vectors ^ (jw.) and v. (jw ) 
aepend on the particular frequency that is selected v^en the SVD 
is acccMjplished. Thus the input-output relationship of Figs. 5 and 6 
are only valid at the frequency Note also that the unit vectors 
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( ju.) and V. ( ju.) ara in ganaral cos^lax but may be realized by 

—3 {J Q 

passive attenuating filters that give the appropriate phase shift 
or tine delay at frequency 


In the SISO case when G(jb)) is a scalar the vectors and y. 
becone the complex scalars u and v which have unit magnitude. Since 
the input-output relationship in these figures is singly the positive 
gain of l/o^ ™ust be that at the phase of the pr<-iduct 

u( ju»Q)v* ( jw^) is simply the negative of the phase of 

-1 -1 -1 

H+g Ii+g(jwQ)] . 


Using Figs. 5 and 6 we may interpret the directional nature of 
the smallest (according to a particular error criterion) model error 
in G(jO)^) that destabilizes the closed loop system. The gain fr<xn 
input to output in Fig. 5, as mentioned before, is simply 1/0^ • 

If the input-output directions u (jw_) and v (jco^) are used and if 
OntjWo^ is small, then a small eunount of positive feedback around the 
system of Fig. 5 will destabilize the system. This is shown in Fig. 6 
where the output of the system of Fig. 5 is fedback to the input with 
a gain of a. Notice in Fig. 7 that if a=0^(ju>Q) the system becomes 


unsted)le because the system anmlifies the input by la (jii)„)-a] 

n 0 


-1 


This additional feedback could be interpreted as a perturbation to the 
system of Fig. 5. However, by block diagram manipulations it Is not 
difficult to see that this perturbed system is equivalent to those of 


Figs. 8, 9 and 10. 


Osa<<rn 


MBtabilizln, F..a»ck in H»=t Sennlniv. Direction 
u ()t>„)v!l3“o> criterion E - G (G-Gl . 
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In Pig. 10 «n sxpllclt invarsion of Il-oni (jW-)v**(j«_)v**(Jw.) ] 

"Ti 0 "T4 0 0 

has bvcn parfonMd to siagillfy tha block diagrin. 

Pron Pigs. 9 or 10 it is claar that tha stability of thasa aqal* 
valent perturbed closed- loop systsias is coi^lataly characterised by 
the behavior of the loop transfer function matrix G(jb)) which at 
is given by 

S(jU)Q) » G(jUQ)L(jWQ) - G(10 )q) Il-Q^(jWQ)vJJ(jWQ)] (44) 

Thus, in the error criterion E(s) ■ G ^(s) IG(8)-G(s)] given in (3.35) 
we have that E(ju>Q) is given by ' 

E (jw ) “ -au (jto.)v^(ju ) (45) 

u ~n u —n u 

This means that the pertiarbation matrix L(jb)^) that perturbe the 
inputs to the open- loop G(ju^) has the same effect as applyii^ 
additional positive feedback in the most sensitive direction 

II 

u ('^u.)v ( ju^) . Just as have interpreted the worst error as ad- 
-n • 0 -n 0 

II 

ditional positive feedback in the direction u (jo).)v (ju.) using the 

—n 0 ~n 0 

SVD of I+G ^(ju^), we obtain similar interpretations using the SVD of 
I-K;(ja)^). Using Pig. 6 we may again determine the smallest des- 
tabilizing feedback as given in Fig. 11 and its open- loop equivalent 


in Fig. 12. 




Osa<0'’n 


Fig. 


11; Destabilizing Feedback in (tost Sensitive 
Direction 

Criterion E “ [G ^-G ]G. 
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Here again the model error criterion of (3.44) gives the model error 
as 

E(ju)-) -= -au (ju)^)v”(jw^) (46) 

0 — n O-n O 

and the perturbation LCju^) as 

L(jw.) = (l-au (jw.)v"(jm )) . (47) 

O — n 0 — n 0 

Thus, by interpreting u.(ju_) and v.(ju ) as input-output 

—DO — D 0 

directions, the model error in a certain direction can be viewed as 

the induced additional feedback in the directions specified by the 

input vectors u.(jm) and tne output vectors v (ju„) . 

~D 0 ~j 0 

Thus to differentiate between thosu model errors that increase 
the m 2 irgxn of stability and those that decrease the margin of 
stability it is necessary to examine their contributions in certain 
input-o:tput directions . 

A model error will decrease the margin of stability of the feedback 

system if it can be interpreted as additional positive feedback in 

the input/output directions u (jw„) and v (ju„), which are equivalent to the 

— n 0 — n 0 

II 

most sensitive model error direction u (ju^)v (joj^), and the contri- 

— n 0 — n 0 

buttons of the model error in the other input/output directions are 
negligible. 

A model error will increase tlie margin of stability if it can 

be interpreted as additional negative feedback in the input/output 

directions u (jw.) and v (ju)„) and the contributions of the model 
— n 0 — n 0 

error in tne other input/output directions are negligible. 


-168- 


The contributions of the nodel error in the other input/output 
directions are negligible, where the SVD of I+G is used, if 

for all i and j not both equal to n, 

II 

ficiently close to v. ( jO)„)G„, { jw„)u. ( jw ) . If the SVD 

— 1 0 CL 0 —l 0 

of I+G(jw„) is used, the perturbed closed- loop system G is 

0 CL 0 

replaced by (I+G(jw )) ^ and G (jci) ) is replaced by [I+G(joj )1 ^ in 

0 da 0 C/ 

the previous sentence. The contributions of the model error in the 

other input/output directions must be negligible because they may 

potentially cancel out the effect or contribution of the model error 

in the input/output directions u (jOJ,^) and v (jo),^). 

— n O — n O 

4 . 4 Example of Section 3.3 Continued 

In this section, the example of section 3,3 is reconsidered 
and it is shown how the model error given in Fig. 3.16, that 
destabilizes the feedback system, can be predicted by computing the 
smallest destabilizing or worst model error by the methods of this 
chapter. Also, the class of modelling errors is restricted to 
completely exclude this type of worst or smallest destabilizing error, 
and the next worst or next smallest destabilizing error is computed. 

The size or norm of this error is given by Theorem 1 and its structure 
is given by (13) . 

These computations are displayed graphically for both the relative 

-1 ~-l -1 

and inverse relative error criteria (i.e., E«G (G-G) and E=(G -G )G). 
To nukke a comparison of the results with the different error criteria, 


Boda-like plots of the elements of L(jo))t the metr^x that makes I-tG ( ju) L ( ju) 
singular at every u, are given. These L(ju) matrices correspond to the 
minimum modelling errors of a specific criteria within an appropriately 
restricted class of modelling errors. As in the SISO case, illustrated 
in Fig. 2, the pertubed systot G(jai) needs only to correspond to 
G(ju3)L(j(t)) (for the L(jii)) that nuJees I-K:(ja)}L(ju)) slnguleu* at every (a) 
at a single frequency in order to destabilize the feedback system. The 
magnitude of these model error is given by the corresponding plots of 
the singular values of I+G(ju)) and I+G ^(jw). 

In Fig. 13 the singular values of I+G(ja>), for our example, as welJ. 
as their geometric mean are plotted in dB versus the frequency, u). 



Frequency (rad./sec.) 


Fig. 13: Singular Value Quantities of I+G(jo>) 
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These quantities determine the magnitude of the smallest 

and next smallest modelling error (when the class of errors 

is restricted by <u (s)v ( s) ,E ( s) >=0) th,»t destabilize 

— n — n 

the system t^en E ■> [G ^]G. The magnitude of the smallest or 

worst error Is given by O . emd the magnitude of the next smallest or 

min 

1/2 

next worst error is given by lo . a ] . The mlnimisn of O . occurs 

^ min max min 

near u^l rad/sec, (0 . = -23dB) ; thus the required magnitude of the 

min 

1/2 

worst error Is -23dB. However, at o)=0, la . a ] = 4dB, indicating 

min max 

that the next worst modelling error occurring only in the frequencies 

edx>ut Is necessarily of a magnitude of 4dB in order to destabilize 

1/2 

the feedback system. Since [0.0 ] approaches OdB as 

min max 

there exists a modelling error In the high frequency range of the next 

worst type that need only have a magnitude of OdB that will destabilize 

the feedback system. Note, however that 0 dB is also the magnitude of 

the worst error if the error is restricted to the high frequency raxige, 

since O . = OdB as 

min 

The nature of the worst error corresponding to O^^^ in Fig. 13, 

Is obtained by plotting elements of 

L(jO)) * (I+E(ju)))"^ (48) 

where 

E(jw) = lG“^(jw)-G"^(jaj)]G(ja3) (49) 

is such that I+G(ju)) Is singular at all u. This is displayed in 
Figs. 14 and 15. Note that in Figs, 14 and 15 that the diagonal 
elements of L(jw) neM ti^l, (where 0 . = -23dB) are essentially iinity 
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while ^ considered essentially zero (=< -80dB in nagnitude) 

and that essentially -,0708. 



Fig, 14: Magnitude Bode-like plots of elements of L(juj) 

for worst model error. 



Phase (degrees) 


Thus at b>"l, L(ji Is given approotliqately by 


-.0708 1 


and represents a crossfeed type of perturbation as in Fig. 3,16 which 
has a constant crossfeed perturbation L given by (with b, »50) 



Frequency (rod./sec.) 


Fig, 15: Phase Bode<'liKe plots of elements of L(jw) for 

worst model error 
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( 51 ) 


which gives a modelling error of -20dB. Therefore, we see that the 
essential nature of the crossfeed perturbation (51) Is detected by 
the approach presented here as evidenced ]3y L(0) In (SO). 

The above discussion points out that a control system designer 
can generate these plots and determine what type of gain changes or 
channel crossfeeds, that veie neglected In his nominal design model, 
should be examined carefully; because If these gain and crossfeed 
errors occur the feedbac)c system can became unstable. The control sys- 
tem designer does not need to worry ahead of time about all the dif- 
ferent types of model uncertainties that might occur; the nature of 
these plots vs. frequency will provide him with guidance with what 
type of modelling errors and in what frequency range he should be 
most concerned with. 

For compaurlson, the plots amalogous to Figs. 13, 14 and 15 using the 

singular values of I+G ^(ju) rather than those of I+G(jo)) and the 

error criteria E = G ^[G-G] = L-I eu:e shown in Figs. 16, 17 and 18 

respe' tively. Note that, for u>=l. Figs. 17 and 18 Indicate neaurly the 

same L(j) as in (50). However, as a . of I+G(jw) or I+G ^(ja>) both 

xnxn 

increase as u increases, that from Figs. 14, 15, 17 and 18 that L(jlOOO) 


Mognitud 



Fig. 16s Sii^ular Value Quantities of I+G ^(jw). 



Magni 



Fig. 17: Magnitude Bode-like plots of elenents of L(3W) 

for worst model error. 




■jr^' 
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asBociated with O^^Il+G(jlOOO)] Is approoclnately given by 
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Thus, aa the tolerable error beoonee larger ee 0 . lZ-K«(j(d)] and 

nxn 

^(ju)] becoBM larger (as shotm in Figs. 13 and 16 respectively), 

the type of errors that the different error criteria characterise laay 

be rather different as (52) and (53) indicate. When both 0 , U+G(jw)l 

min 

^nin^^^ (ju)] are sufficiently small the different error criteria 

guard against the sane type of nodel error or equivalently L(j(i)) as shown 

by this exanple. This neans that either test using o , Il4<3(ja))] or 

min 

^nin^^^ (jw)] will detect the near inrt-ability of a control systen. 

However « they may give rather different estlnctes of gain and phase 
margin when the feedback system will tolerate a class of sK)delling er- 
rors of larger magnitude. 

Now consider, the next vrorst u»del error for the two error 

criteria used in this chapter. Thus, the class of modelling errors 

now considered must exclude the %«orst nodel error type just discussed. 

The model errors now considered will have zero component in the most 

sensitive direction u ( jO))v**( ju>) , (i.e., u**( ju)E(ja))v (ju)-O). 

~n —n -*n •— n 

-- 1 -1 

For the IG -G ]G error criterion we may ^aln draw Dode-like 

1/2 

plots of L(ju) that corresponds to the la . CT J error magnitude 

min max 

in Fig. 13. This is shown in Figs. 19 and 20 «diere the off-diagonal 
elements of L(j ) are not plotted because their magnitudes are 


insignificant. 


Mognitude (dt) 



Fig. 19t Magnitude Bode-like plot of elenents of L(jo)) 
for next v^orst error. 


4 
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Fig. 20*. Phase Bode-llke plots of elanents of L(ju) 
for next worst error. 

Recall that the error matrix E(jw) for the next worst error is specified 
by (13) where 0(jw) is arbitrary but reel. In Pigs. 19 and 20, 0(jw) 
has been ret to zc' -> in order to calculate a single L(jw). Prom 
Figs. 19 and 20, it is clear that the next worst type of error is to 








Frequency (rod./sec.) 


Fig, 22; Phase Bode-liJce plots of elements of L(jO)) for 
next worst error. 

Note tliat in this case, L(ju>) increases the gain in both feedback chan- 
nels while changing the phase of both channels. Thus, once iDore, we 
see from Figs. 19-22 that the model error, or L(jb)}, the criteria guard 
against are essentially different. 
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4.5 Improving Robustness Tests by Combining Tests 

It has been dennnstrated in this chapter that the usefulness of 
the various robustness tests of chapter 3 can be extended by restricting 
the class of allowable model errors so that the error structure is 
exploited. In this section, it is shown that the upefulness of the 
robustness tests can be extends! by a ccxnblnation of two tests forming 
a single hybrid test. This effectively enlarges the class of allowable 
model errors and, therefore, in certain circixnstances the stability 
of the feedback system may be confirmed when either test alone would 
fail to be conclusive. 

The basic idea is to use tests that employ the same G(s,e) in 
their proof via Theorem 3.2; one uses one test for a certain subset 
of frequencies and uses the other test for the remaining frequencies. 

The reason that this procedure works is that both tests guarantee 
exactly the same thing , that I+G(s,e) is nonsingular for all 6 in 
10,1] for any s in their respective subsets. Let D , the Nyquist 

iK 

contour of Fig. 3.10, be decomposed into two subsets D and D 

XR aR 

whose union is and let TESTl and TEST2 denote any of the tests 
bounding the model error magnitude in the theorems of chapters 3 and 
4 that employ the same G(s,6) in their proof. Then the following 


theorem may be stated. 
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Theorem 4 ; The polynomial ^q^(s) no OUIP zeros and hence the 
pert'irbed feedback systmr is stable if the following conditions Itoldi 


1. Condition 1 of Theorem 3.2 holds 


2 . 


3. 


4. 


°2R 


TESTl implies I+G(s,e) 


TEST2 implies I-i-G(s,e) 


is nonsingular on D^^xIO, 1] 
is nonsinguleu: on D x[0,l] 


Proof ; Condition 2 , 3 and 4 guarantee condition 2 of Theorem 3 . 2 and 

therefore by Theorem 3.2, i (s) has no CRHP zeros. Q.E.D. 

CXt 


Theorem 4 allows us to combine tests that employ both absolute 
and relative error measures as well as tests that utilize model error 
structure and those that do not, provided they can work with the saune 
I+G(s,e). Even in the case where the tests were originally derived 
by use of different G(s,e) matrices it is sometimes possible to find 
a single G(s,e) from which versions of the original tests may be 
derived via Theorem 3.2. For example, a version of Theorem 3.6 is 
derived in 151] by use of the G(s,e) = (l-e)G(s) + eG(s) which is 
used in the derivation of Theorem 3.4. 

This version of Theorem 3.6 requires more complicated conditions 
on the allowable L(s) than the eigenvalue restrictions on L(s) in the 
present Theorem 3.6. However, these more con^Jlicated conditions are 
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autonatically satisfied provided 0 (L ^(s)-I)< 1. Therefore, 

max — 

Tinder this restriction on L(s), Theorem 3.4 may be used on and 
Theorem 3.6 on to prove the stability of the feedback system 

under variations in the system model G(s). 


Bxegnple 2 ; To show how tests may be combined reconsider the example 
of section 4.2 where G(s) is given by 


G(s) 


JL 

8+7.5 

0 


0 

_1 

s^0. 5 


(54) 


and where we use the same constraints on the model given in (22 ) , 
(23), (24) and (25). Notice from (26) that 


O (E(ju))< 2.5 (55) 

max — 

and that frcxn (20) for all W>1 

a . (I+G'^(ja») = |l.5+jul> 2.5 > 0 (E(ju)). (56) 

nun — max 

Therefore we could use Theorem 2 for all u<l and Theorem 3.4 

which employs the test in (56) for all u)>l to ensure the stability 

of the feedback system. 

Obviously, the division of the frequency axis could be csurried 
out for n different tests which divide the frequency axis into the 
n subsets whose union is the whole frequency axis. In fact, as n 
becomes large this suggests that a single test that depends contin- 
uously on s could be devised. This was mentioned before in the context 


In this section, we discuss a few key relationships between the 
nominal closed-loop system, denoted by and the quantities 

involved in confuting bounds on the allowable model error and the 
most significant error structures considered in this chapter. 

We first make the simple observation that if the SVD of a squaire 
matrix A is given by 

A = UZv” (57) 

then the SVD of A ^ (assuming it exists) is given by 

a"^ = V u“ . (58) 

This relationship between the SVD's of A and A ^ is useful vdien 

the quantities a . (I+G) , a . (I+G ^) and a . I (I+G) (I-G) , 

min min min 

used in Theorems 3,4, 3.6 and 3.9, are required. The relationship 
between these quantities can be determined from the following equations; 


^~cl = 

(59) 

(Gj,j-1/2I)“^ = -2(I+G) (I-G)”^ 

(60) 

(Gj,j^-I)'^ = -d+G) 

(61) 
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where the nominal closed-loop system Is given by 

- G(I+G)"^ • (62) 

Thus by calculating the SVD’s of G__-aX for a>0f a"l/2 and 
one can easily obtain via (59), (60) and (61) the SVD's of I+G, 

I+G ^ and (I-G) ^(I+G). This is significant since G__ need only 
be calculated once for the three SVD's; also, no explicit inversion 
of G or I-G is needed since if 

G(s) - C(ls-A)"^B (63) 

then G (s) is given by 

CL 

G_. (s) = C(Is-A+BC)~^B (64) 

CL 


Another ccxnputational saving can be realized if only approximations 

to the minimum and maximum singular values of a matrix and not its 

full SVD are required. Recall from (2.38) and (2.39) that O . (A) 

mxn 

and a (A) are given by 
max 




.-li 1-1 


(65) 


and 


= 1 A) 2 • 

max z 


( 66 ) 


Using matrix linear algebra, it can be shovm that 


-5- llAll, < l|A|l < 
/n~ 


( 67 ) 
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and 


^ llAl 




Icol 


I 2 1 


( 68 ) 


if A is an nxn matrix. Since the matrix norms | |a| and | |a| 

au:e much cheaper to confute than | |a| the singular value quantities 

a . (I+G), a , (1+g“^) and a . I(I+G) (I-G)'^] may be approximated by 
min min min 


0 . (I+G) 

= 1 |g_.-i| L 

-ill? 

min 

I ' CL "2 

' ' CL ' ' 1 

. -1 

t 1 1 1 “ 1 

II 1 1 — 1 

0 . (I+G ) 

= G 

= L 

min 

' ' cl ' '2 

' ' CL ' ' i 


1 , 1 1 

1 . 1 1 — 1 

a . ni-G)' 
min 

^(I+G)l- j 1 



(69) 

(70) 


iig^-i/2iii:^ (71) 


within the bounds given by (67) and (68) when i=l or «*>, for 
ccmtputational savings. Note that as n increases these approximations 
may be poor. However, all the robustness theorems of chapter 3 can 
be formulated using the 1 or * norms rather than the 2-norm 
(singular values) and approximations need not t>e used at all. 
Nevertheless, the results on the structure of model errors are only 
applic2d3le in the case of the spectral or 2-norm. 


4.7 Concluding Remarks 

In this chapter we have shown how the structure of the model 
error may be used to improve the theorems of chapter 3. The nature 
of the worst model error was explored through the use of block 
diagrams showing that it is equivalent to additional positive feedback 
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in the Input/output directions and The exaoyple 

of section 3.3 was used to illustrate that the nature of the snallest 
destabilizing error can be obtained from Bode-like plots of L(j(0) . 
Also in^roving robustness tests by ccnbining different robustness 
test over different frequency ranges was discussed as well as com- 
putational considerations for the efficient computation of quantities 
required by the robustness tests of this chapter and those of 
chapter 3. 

To place these results in perspective, a design/analysis pro- 
cedure is suggested. This procedure assumes that some particular 
synthesis procedure, such as the LQG methodology, is used to obtain 
specific controller designs. An outline of the procedure is the 
following sequence of steps: 

1. Obtain an initial controller design that meets basic 
performance requirements but does not produce a con- 
troller with a bandwidth larger than the upper fre- 
quency limit for which the model is valid. 

2. Obtain an estimate on the allowable model error 

magnitude ae a function of frequency and cc»^are 

with the values of 0 . of I+G or 

min 

2.1. If the model error magnitude is 
of I+G or I+G , stop. 

. If the model error magnitude is 
of I+G or I+G ^ then go to step 


I+g"^. 


less than a 


min 


larger than 0 


min 


2 . 2 . 1 . 


2.2 
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2.2.1. Calculate the worst siodel error and check to see 
if this error could possibly occur. If not go 
to step 2. 2. 1.1 otherwise go to step 2, 2.1. 2. 

2. 2. 1.1 Calculate the raagnitude of the next worst error. 

If the model error magnitude in step 2 xs less 
than this, stop. If not, Canute the magnitude 
of the "next worst error" and continue with a 
step similar to 2.2.1 etc. (this gets rather 
tedious ! ) . 

2. 2. 1.2 Try to improve model to reduce model error or 
chaise controller design. Return to step 1. 

Exactly how to change the model to reduce model error based on 
the analysis methods of this chapter is an open research question. 

Also, it is not always clear how changes in controller design may af- 
fect the quantities a . (I+G) and o . (I+G . This is also an open 

min min 

research problem. However, in spite of these difficulties (which in 
the author's opinion may eventually be adequately circumvented) the 
key to making practical use of the results of this chapter depends 
on the engineer's ability to determine whether the model error mag- 
nitude in the most sensitive direction, i.e. , the projection magnitude 

|<u (ju)v*^(jw), S(ju))>|, can be bounded. Cleeurly, engineering judgement 
is necesarry; however, this engineering judgement may not easily 


translate into this type of bound. Thus, practical experience in 


-190 


obtAlnii^ these type of bounda is neceaaeiry In ordar to further 
iaprove the type of rohuatness teata that are approprie.te for the 
kinds of knowledge about model uncertainty tha^. the designer has 
at hla disposal. 
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5. ROBUSTNESS AND LQG CONTROL SYSTEMS 
5.1 Introduction 

The previous two chapters dealt with a loop transfer matrix G(s) 
that contained compensation as well as ooen-loop plant dynasdcs. 'Shs 
robustness results of these chapters hold irxiependent of tlM MIMO design 
methodology used to determine the ccmpensation required. 

This chapter will be concerned with deriving robustness result for 
feedback control systems designed using the linear-quadratic-gaussian 
(LQG) design metl»dology [41]. This includes results for the linear- 
quadratic (LQ) state feedback regulator and some of its variations as 
well as the LQG regulator. The multivariable version of Kalmw's inequality 
and Theorem 3.6 form the basis for the derivation of these results. 

In section 5.2, the LQ and LQG control problems are stated for 
completeness and the definition of t)» loop-transfer matrices for the 
feedback systems is given. Section 5.3 continues with a discussion of 
the multivariable Kalman inequality derived fron the Riccati equation. 

The stabilit' irgins for LQ regulators is then discussed in section 5.4 
where it is sho'wn that these regulators have guaranteed minimuB stability 
margins which makes them attractive and that the control weighting R 
matrix determines the coordinate system in which these stability margins 
hold. Stability margins for the state feedback regulator, whose feedback 
gain is determined by a Lyapunov or a nonst^mdard parameterized Riccati 
equation, are also given. Section 5.5 concl\xies with a discussion of the 
stability margins for LQG regulators. It is shown that stability margins for 
LQG regulators are the same as those for LQ state feedback regulators but 
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only at a point inai^ t)M LQG coi^naator. That# Kargina cannot bo 
autonatically guarantaad at tha ^tyaical input or output of tha plant 
unlaaa tha Kalaan filtar of tha LOG coo^^naator haa an aaact s»dal of tha 
parturbad opan-loop plant, a vary raatrictiva aaaia^ion. Howavar, 
r<^uatnaaa racovary procaduraa (35, 36] ara diacuaaad that allow a proparly 
deaignad LQG ccHitrol ayatm to aaymptotically reco*/ar tha IQ atata faad- 
back atability margina at tha input or output to tha phyaical plant pro- 
vided that tha plant ia miniann phaae. The significance of having Kargins 
at various points in tha feedback loop is also discussed with rafaranca 
to modelling error characteriaation. 

Some of the results of this chapter have ai^arad previously in tha 
literature. Based upon the preliminary gain margin results in Ncmg aiwl 
Athens [27] , Safonov and Athans [25] gave tha definitive treatment of 
guaranteed minimum multiloop stability margins for the LQ state-feedback 
regulator allowing for nonlinear perturbations in the feedback loop. It 
was later shown by Doyle [33] that there are no guaranteed minimua sta- 
bility margins for LQG regulatois. Kwakernaak [36], and Doyle and Stein 
[35] have outlined procedures whereby the LQG regulators may asymptotically 
recover the LQ regulator guaranteed margins. 

The cemtributions of this chapter are mainly the simplified derivation 
and characterisation of LQ regulator stability margins. The structure 
(diagonal vs. nondiagonal) of the control weighting R matrix is shown to 
have important ing>act on the stability margins. If R is nondiagonal the 
LQ regulator may i‘ave arbitrarily small gain margins If R is not a 
multiple of the identity matrix it is shown that ability to tolerate 
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crossfced perturbation! is drastically reduced. Ibese results also apply 
to the variations of the IQ regulators discussed in Section 5.4. Also, 
in section 5.5 an inequality is derived that ens«ures that stability margins 
will apply at both the input and output of the j^ysical plant. 

5.2 The LP and LQG Regulators 

For the sake of completeness, the LQ ^utd LQG regulator pr^lems and 
their solution will be given for the linear-tiae-invariant , infinite time 
horizon case [41]. 

5.2.1 LQ Regulator Problem 

For the open -loop plant given by 

x(t) « Ax(t) + Bu(t) (1) 

find the optimal control u*(t) that minimizes the quetdratic cost functional 
J (u) given by 

QD 

“ / [x^(t)Qx{t) + u'^(t)Ru(t)]dt (2) 

0 

where Q ^ 0, R > 0, IA,B) is stabilizable and [A, is detectable. 

The c^timal control u*(t) is given by 

u*(t) ■ -r”^b\ x(t) (3) 

where K > 0 satisfies the algebraic Riccati equation 

A^K + KA + Q • 0. (4) 


The block diagram of this regulator control system is shown in Pig. 1. 


— R’’b’k — 

FEED^CKGAIN 


Pig. 1: LQ Regulator 


To calculate the loop transfer matrix G(s) at the input to the plant, 
we break the loop at point in Fig. 1 so that now and are no 
longer equal. Next we calculate the tramsfer function matrix from jij^ to 
-U 2 - This is G(s) , the lo<^ transfer matrix at the plant input and is 


given by 

G(s) = r'^’^KCIs-A)"^ (5) 

To calculate the loop tramsfer matrix at the output of the plant (at point 
of Fig. 1) we follow an analogous procedure. We break the loop at 
point Fig. 1 and calculate the transfer function matrix frcati x^^ 

to -x^ which is given by 


-X 2 (s) = (Is-A)"^B r'^B^'k Xj^(s) . 
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-1 -1 T 

Thus (Is-A) BR B K is the loop transfer matrix at the output of the plant. 

In general, if G (s) denotes an open-loc^ plant with ii^nit u(s) and 
P ~ 

output /(s) and u(s) = -G^(s)^(s) where G^(s) represents the transfer 

matrix of a compensator, we call G^(s)G^(s) the loop tramsfer matrix at 

the input and G (s)G (a) the loc^ transfer matrix at the output, 
p c 

5.2.2 LQG Regulator Problon 

Let the stabilizable, detectcihle c^>en-lo<^ plant be given by 


x(t) = Ax(t) + Bu(t) + ^(t) 


(7) 


^(t) = Cx(t) + 0(t) 


( 8 ) 


where the noises ^(t) and ^(t) are both Gaussian, white, zero mean, mutually 
independent and stationary with 


E(C(t)CNT)] = E 5(t--t) ; E ^ 0 


(9) 


E[e(t)e^{T)l = 0 6(t-T); 0 > 0 


( 10 ) 


Find the optimal feedback control u*(t) depending causally on ^(t) that 
minimizes the quadratic cost functional J (u) given by 

T 

J(u) = lim ^ r [x^(t)Qx(t) t u^(t)Ru(t) ]dt 
T^ ^ 0 

where Q ^ 0 emd R > 0. The optimal control u*(t) is given by 


( 11 ) 


u*(t) = -G x(t) 


( 12 ) 


where x(t) the state estimate is generated by the Kalman filter 


x(t) = Ax(t) + Bu*(t) + G,,[^(t) - Cx(t)] 


(13) 


or equivalently 


x(t) =■ [A - BG^ - G^Clx(t) + G^(t) 


(141 


with G^ euid G^ being given by 


-1 T 

G^ - R B K 


G^ = E c'^0'^ 


(is: 

(16: 


with [A, detectable, (A, 5^^^] stabilizable, K > 0 satisfying 


A K + KA + Q - KBR 


■Vk = 


(17! 


and E > 0 satisfying 


AE + Ea’’ + 5 - E c'^0"^C E = 0 


(is; 


From (12) and (14) we see that the transfer function matrix from ^(s) to 
u(a) is given by 


- 1 , 


u*(s) » -[G^flS-A + BG^ + G^C) G^]^(s) 


(19: 


and thus the blocli diagrcmi of the LQG regulator is given by Fig. 2. 



cds-Ar^B 


y(s) 


Gr(ls“A+B6|-+GfC) ^G^ ■■ 


Lg. 2; 
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Froin Fig. 2 , we see that the loc^ transfer matrix at the input, G(s) , is 
given by 

G(s) = G^(Is-A + BG^ + G^)~^G^C(Is - A) . (20) 

5 . 3 Multivariable Kalman Inequality 

The subject of qualitative feedback properties of LQ control syst^s 
is not a new one. An early and fundamental paper by Kalman [21] detailed 
properties shared by all LQ regulators in the single-input case. Kalman 
showed that the scalar return difference transfer function of a single-input 
LQ state feedback regulator satisfies the inequality 

ll+g(ju))| ^ X ; for all u) . (21) 

This is both a classical condition for the reduction of sensitivity at the 
feedback input to the system (see, e.g., [19]) as well as necessary and 
sufficient for a (stable) state feedback regulator to be optimal with 
re aspect to scane quadratic cost index. By inspection of the Nyquist diagram 
corresponding to (21) , (Fig. 3 with 6=1) , it is straightforward to observe 
[22, pp. 70-76] that a SISO LQ state feedback regulator has a guaranteed 
infinite upward gain margin, at least a 50% gain reduction margin eUid also 
a guaranteed minimum phase margin of + 60® . (These margins were defined 
in Section 3.7). 

Anderson [23] developed a multivariedale version of condition (21) as 
a property of LQ state- feedback regulators; a similar generalized condition 
arises in sensitivity theory^ (see, e.g., Cruz and Perkins [24]). In this 

.Sensitivity refers to the variation in system responses due to infinitesimal 
changes in the nominal system parameters. Robustness refers to tha de- 
lineation of finite regions of allowable variation in nominal system param- 
eters that preserve stability. 
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Theortan 1 (LQ Kalman Inequality) : If the matrix K satisfies the matrix 

algebraic Riccati equation 

a\ + KA + Q - KBR'^'^K = 0 (22) 

with R > 0 and Q > 0 then 

(1+G(s) )” r(I+G(s)) » R + H(s) (23) 

where 

G(s) = r“^b'^K{Is-A)“^B (24) 

H(s) = [(Is-A)“^B]”(Q + 2Re(s)K) [(Is-A)“^l . (25) 

Furthermore, If Q > 0, B has full rank and K ^ 0 then (23) implies that 

(I+G(s))”r(I+G(s)) > R, ® ® ^26) 

Alternatively, if det(jcoI-A) ^ 0, for all o), ctnd K ^ 0 then (23) implies 
that 

(I+G(s))“r(I4G(s)) > R, s e D„ . (27) 

R 



(29) 


- 200 - 


Adding R to both sides of (29) gives (23). Now Q -t- 2Re(s)K will be 
positive semidefinite for s 6 if Q > 0 and the Indentations of (2^^ 
are sufficiently small or if Re(s) ^0, s € which happens if det(jul-A) 

0 for all ( 1 ). Thus under these conditions H(s) >0 or H(s) ^0 re^ectively 
for all s e Dj^. 

Q.B.D. 

It is important to point out that this theorem uses G**(s) rather 
T 

than G (-s) as in [23] . These two qucuitities are the same vdien s » jo), 
but are different when Re(s) ^ 0. This is the case when s is evaluated 
along the Nyquist D contour cuid this contour is indented along the 
imaginary axis (Fig. 3.10). It is necessary to use G (s) in order to 
apply the theorems of Chapter 3. Note, however, that when det(ju)I-A) / 0, 
for all to , (i.e., when the open- loop system has no poles on the ju-axis) , 

that Sip is just the imagineury 2 ucis from -jR to jR. In this case (27) 
could be written as 

(I+G(j 0 )))\(I+G(jU)) ) > R; for all o) , (30) 

which is the previously mentioned multivariable generalization of condi- 
tion (21) . 

5.4 Stability Margins of LQ Regulators 

We can now ^ploy Theorem 1 in conjunction with the results of 
Theoran 3.f> to estadslish the robustness properties of multivariable liQ 
regulators. Recall from Theorem 3.6 that one of the key quantities for 
multivariable robustness analysis is the minimtmt singular value • 
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where G(s) is the loop trwsfer matrix. Unfortunately, the inequalities 
(26) and (27) of Theorem 1 do not provide a bound on o . (I-fG(s)), where 

' nuLii 

G(s) is the IQ regulator loop transfer matrix defined by (24). However, 
if we define 

G(s) -= R^G(s)r"^ (31) 

then (27) (for example) can be re%irritten in the form 

(I+d(s))”(I+e(s)) >1, s 6 . (32) 

Equation (32) provides the bound 

o . (I+G(s)) >1, s e R_ (33) 

min — R 

on the minimum singular value of I+G(s). 

To work with G(s) instead of G(s) , it is necesseiry to manipulate the 
system of Fig. 4 into the equivalent (for stability emalysis) form depicted 
in Fig. 5. Then using (26) emd (27) together with Theorem 3.6 leads 
directly to the following result. (Recall from (3.27) that (s) = 
det(sI-A+BC) where G(s) = C(Is-A) ^B, was used in Theorem 3.2 to 
determine the stability of the perturbed closed- loop system.) 

Theorem 2: The polyncwnial (s) has no CRHP zeroes provided the fol- 

■■ CL 

lowing conditions are satisfied: 

1 . (a) ({) (s) and $ (s) have the same number of CRHP zeros 

CL OL 

(b) if ~ ° ~ ° 

(c) ({)__ (s) has no CRHP zeros 

CL 
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2. G(s) is specified by (4) where K>0 satisfies (2) and [A,B] 

is stabilizable , Q^] is detectable and B has full rank. 

3. With y(s) » a tR^L"^(s)R~^-I] either of the following hold: 

IRdX 

(a) Q>0 and Y(s) <1, ® ® ^ 

(b) ^ 0 ^ Y(s) ^ * 


Proof ; It is well-known that condition 2 ensures that has no CRHP 

zeros. Defining 6(s) * R^G(s)R we see that G(s) has a state-space 
realization [A, BR , R B K] and thus its open- and closed- loop charac- 
teristic polynomials identical to those of 

(A,B,r"^B^K) . Thus any assumrtions about and (s) obviously 


•^CL' 


apply to (s) and $ (s) . Similarly, by defining L(s) ■ R^L(s)R 

Oil CL 

we may work with G(s) and L(s) instead of G(s) and L(s) . The conditions 

A 

(26) and (27) of Theorem 1 are equivalent to o^^[l+G(s)J > 1 jmd 

a . (I+G(s)) > 1 respectively. The condition 3a and Theorem 1 require 
min — 

that 


0 (l" (s)-i) < 1 < 0 . (i+G(s)), sen 

mcuc — "^in t 


(34) 


and by Theorein 3.6 we conclude that CRHP zeros. Alterna- 

tively condition 3b and Theorem 1 require that 


o (L (s)-I) < 1 < 0 , (I+G(s)) 
max - min 


( 35 ) 


which again by Thecran 3.6 meems ^ (s) has no CRHP zeros. 

CL 


Q.E.D. 
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Note that the condlticm c (R^L ^(s)R ^-l) < 1 in condition 3a 

SwX 

can be rewritten as 

RL(b) + L**(8)R - R > 0, * ® ^ 

or with s ■ j(A' 

L(ju)R"^ + R"^”(jo)) - r"^ > 0; for all u). (?7) 

The inequality (37) was used by Safonov and Athans [25] to prove the L(2 
state feedback guaranteed gain and phase margins although their method, 
of proof is quite different. Thoy unplicity assume that L(ju)) is stable, 
something which we do not require. 

Theorem 2 can new be employed to establish the guaranteed minimum 
multivariable gain and phase margins associated with LQ regulators. 

Important Remaurk ; We emphasize that these margins are guaranteed 
only if the control weighting matrix R is chosen to be a diagonal matrix; 
we will svibsequently present an example showing that ’ ''e margins can be 
m«fle arbitrarily small for an ai^ropriately chosen non-diagonal R matrix. 

Corollary 1 : The LQ regulator with loop transfer matrix GCr) satisfying 

(16) with a diagonal R>0 has simultaneously in each feedback loc^ a 
guaranteed minimum gain and phase margins given by 


GM 

Dli. ») 

(38) 

PM 

D [-60", eo’] 

(39) 


if Q>0 euid 
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GM D (i, •) 

(40) 

PM D (-60*, 60") 

(41) 


if Q ^ 0 and 0 ^ ( ju)) ^ 0 tor all o). 

Proofs From Theorms 1 and 2 wo know if Q > 0 then 

- o^_(l"^(s)-1) <1, 8 e n„ (42) 

nlaX pIaX — K 

satisfy ng condition 3a when L(s) and R are diagonal. If 0^. (ju') ^ 0 

Cm 

for al! (■' then 

0 . (l"^(s)-i) <1, 8 e Q_ (43) 

satisfying condition 3b when L(s) and R are diagonal. The remainder of 
the proof is completely analogous to Corollary 3.2. 

Q . S . D. 

Note that the margins of Corollary 1 are based on the inequality 
(42) . This inequality will not hold for all s C for "real world" 
modelling uncertainties. In the SISO case, this is clearly tononstrated 
by the f^ysical fact that the phase of g(jii 3 ) is completely uncertain at 
high enough frequencies. This means that for some Idw^) is real 
and negative — that is, there is a 180* phase difference between g(ju)Q) 
emd gCjui^) . If ^( 3 ^^) is real and negative, then |t - ij > 1 

and (42) is violated. This meems that a physical system cannot actually 
have an infinite upward gain margin because its Nyquist diagram always 
will cross t^ie negative real axis at some sufficiently large frequency. 
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Results related to Corollary 1 have been derived by various authors 
[261 - [29] ; but the definitive treatment: including the multivariable phase 
margin result ir due to Safonov and Athans [25] . The ai^roach of this 
thesis, based on relatively aimple frequency domain arguments, is new. 

If R is not diagonal then the guarantees of Corollary 1 do not 
apply. The following example illustrates that the gain margins may became 
arbitrarily small. 


Example 1 ; 

Consider the LQ regulator specified of Fig. 1 with 




(44) 


where is the 2x2 identity utatrix and R '* 0 is a nondiagonal control* 
W3ighting matrix given by 


T -2 -1 -1 

R = B IK +2K ] B , 


(45) 


where K'*0 is arbitrary. By section of R in (45) , K satisfies (2) . Now 
let the rauliplicative pertur)5ation I <s) be given by the constant matrix 
L where 



(46) 


and 6 # 0 is arbitrary. The zeros of $ (s) are the eigenvalues of the 

CL 

perturbed closed- loop system matrix A where 


A„ » A-BLb"^{2I+k‘^) 

CXi 


(47) 


or 
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(p^+2) 6 e P 2-1 

U+€)P2 


ifhere w« have let K ^ be denoted by 


p^ + e e (Pj^?) 

{!+«) (Pj+2) - 1 


( 48 ) 



( 49 ) 


For to have no C3tHP eigenvalues it is necessary for tr A^^^ < 0. How- 
ever, by inspection of (48) , if then for any Si*'' there exists a 

B that will ma)ce tr A > 0 and therefore for arbitrarily small 6, the 

CIs 

perturbed closed-loop system will be unstable. Thus an design can have 
an arbitreurily mall gain margin. 

The basic probl^ explosed this example is that the margins are really 
guaranteed at a different point in the loop than where we would like. 


TTiis is i.1 lustrated in Fig. 6 where the perturbation L(s) is inserted at 
point . When £(s) is diagonal, as wlren calculating gain and phase 
margins, and R is also diagonal then R ^ and L(s) commute and points 
amd have identical guaranteed gain and phase margins. Point is 

where it is important to have margins (1 .e , at the input to the physical 
plant) , not inside the compensator at point . If R is not diagonal 


and 


.-il 




1, a "small" perturbation at point stay look 


like a "large" perturbation at point o because of the amplifying effect 
of the nondiagonal R matrix scaling. 

Returning to Example 3.1 of the Chapter 3 once isore, an LQ feedback 


control law is given that has the same closed- loop poles as before. 





pi L-J R‘^ k 


Fig. 6. LQ regulatci.' with margin: ouarantefid at point 

for a-i R > 0 emd at both and for diagonal 


R > 0. 


but avoids the near instability associated with the negative identity 

feedback. This example shows that with R=I, o . (I+G(ju>)), the \:^per 

min 

bound on tlie allowable magnitude of modelling error of (3.44) given by 

a (L ^(jo)) - I), is autcanatically greater than or equal to unity, 
max 

Example 2; 

With b = 50 in (3.20) as in the plot of a . (I+G(jo))) in Fig. 3.24, 
12 min 

an Lg design using R=I and 


Q = 3 


gives a feedback gain of 






-209- 


rVK 


1 -50 


Lo 1 J 


and a closed- loop system matrix of 


A„ = A - ® -21 

CXj 


(51) 


(52) 


This makes I-K3(s) 


I+G(s) 


s-t-2 

s+1 

0 


0 

s+2 

s+1 


and thus 


(53) 


> 1 . (54) 

As one might expect the ability of LQ regulators to tolerate cross- 
feed perturbations defined in Section 3.7.2 is also affected by the choice 
of the control weighting matrix R. This is made precise in the following 


0 . (I+G(ja)) = 
min 


I 2 , 

1 0 ) +4 

1 0)^+1 


corolleury. 


Corollary 2 ; The LQ regulator with loop transfer matrix C(s) satisfying 
condition 2 of Theorem 2 will tolerate (i.e., (s) will have no CRH? 

CXj 

zeroes) a crossfeed perturbation of the form 



X(s)‘ 

"I 

0' 

L(s) = 


or 

(55) 

.0 

I , 

.X(s) 

I . 


where 
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o (X(s) ) < min 
max 


! X . (R) X . (R )) 
min 1 min 2 f 

X (R.) 'X (R,)( ' 

max max 1 ; 


s e Q, 


(56) 


provided condition 1 of Theorem 2 is satisfied and where R is given by 


R = 


Rj^ 0 


(57) 


and is conformably partitioned with L(s) and either Q > 0 or (jw) ^ 0 

OL 

for all (i) holds. 


Proof ; Only conditions 3a and 3b of Theorem 2 need to be verified for the 
L(s) of (55) the rest are satisifed by assumption. Note that for s 6 


a (R^l“^(s)r"^-I) = a 

max max 


0 

.0 0 J/ 


or 


max 


-R^ X(s)R~^ 


(58) 


< a (X(s)) max[0 (8^)0 (R~S . O (R^)0 (r7^)] 

— may n>ax 1 max 2 max 2 max 1 


and hence if 


a (X(s)) max 
max 


X* (R,) 
max 1 

X^. (R_) 
min 2 


X (R.) 
max 2 


X’. (R,) 
min 1 


< 1 


(59) 


then conditions (3a) auid (3b) are both satisifed. However, (59) is 


equivalent to (56) . 


Q.E.D. 


Note that 


X min(R) , 

i Xmin(R^) 

Xmin(R2) ) 

^ Xmin(R) 

X max(R) — 

jxmax(R2) ' 

Xmax(R^) 1 

— Xmax(R) 

m m 


which indicates that if the ratio of Xmin(R) /Xmax(R) is very small that the 
ability to tolerate crossfeed perturbations is drastically reduced. As 
illustrated in Fig. 6 the use of R scales the inputs euid outputs such that 
the stability margins are obtained in the scaled syst^ rather than the 
original system. This means that if our original inodel has the coordinate 
system in which we would like to guaremtee margins, that R should be 
selected as R = pi for some positive scalar p. 

The effect of the R matrix on the tolerance of the closed-loop system 
to general modelling errors of the form of (3.44) can i)e accounted for by 
using the inequality 

To guarantee stability via Theorem 2 we must have 


a (R^l”^(s)r'*-I) < 1 
max — 


(62) 


which is ensxired if 




From (63) , it is clear that the tolerance to model error may be reduced 
by a factor ^Xmin (R) /Xmax (R) from the case where R ■ pi for a positive 


scalcu: p when equality holds in (61) 


5.4.1 Variations of IQ Design s 

Since LQ designs have inherently good margins provided R is selected 
appropriately, it is natural to search for variations of this method. One 
such variation, proposed by Wong and Athans [27], is to soxve a Lyapunov 
rather than a Riccati equation to compute K used in (3). 

The Lyapunov equation with Q ^ 0 given by 

a’’k + KA + Q = 0 (64) 

guarantees that the eigenvalues of A lie in the OLHP if K ^ 0 and [A, Q^] 
is detectable. The corresponding Kalmar, type inequality for lo<^ transfer 
matrices G(s) specified by (24) where K ^ 0 satisfies (64) is given by 

RG(jco) + G*^(jto)R > 0; for all 0 ) (65) 

amd is the furdamental inequality used to derive stability margins. When 
Q > 0 the inequality (65) becomes strict. The stability margins for this 
type of feedback are given in the next theorem and its corollaries. 


Theorem 3 ; For G(s) of the form of (24) , $ (s) has no CRHP zeros if the 

CXi 

following conditions hold: 

1 . ( s) has no CRHP zeros 

C/Ii 

2. K > 0 satisfies (6%) With Q ^ 0, R > 0 and [A, detectable, 
and B has full rank 


3. either of the following holds 


(a) 

(b) 


Q > e and RL(s) + L^(s)R > 0, s e 
RL(s) + l”(s)R >0, s e Q_ . 

K 
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Proof; Conditions 1 and 2 and the Lyapunov stability criterion gucurantee 
that condition 1 of Theorem 2 is satisfied. As in the proof of Theorem 2 
we may work with G(s) « R^G(s)R ^ and £{s) « R^L(s)R ^ instead of G(s) 
emd L(s). Condition (65) is simply condition 2 of Theoron 3.8 with G(s) 
replacing G(s), and L(s) replacing L(s) in its condition 3 is simply 
condition 3b. Thus by Theorem 3.8 the theorem is proved «>hen condition 
3b holds, \then Q > 0 and condition 3a is satisfied, the strictness of 
tlte inequality of condition 3 of Theorem 3.8 may be changed to ^ and the 
^ of its condition 2 to > and Theorem 3.8 r^ains valid. Thus when condi- 
tion 3a holds the theorem is proved. 

Q.E.D. 

Corollary 3: For G(s) as in Theorem 3 with R diagonal the guareuiteed 

gain and phase margins are given by 


GM 3 [0, ®) 

(66) 

PM D (-90®, 90°] 

(67) 

condition 3a of Theorem 3 holds and 

CM 3 (0,®) 

(68) 

PM 3 (-90°, 90°) 

(69) 


if condition 3b oZ Theorem 3 holds. 

Proof ; Similar to Corollary 3.5. 

The importance of Corollary 3 is that the standard LQ guaranteed 
gain reduction margin of i can be reduced to 0 by using K satisfying the 
Lyapunov equation (43) with Q > 0 rather than the Riccati equation (2). 
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Of course, it is possible to have a zero gain reduction nuurgin only for 
open- loop stable syst^s. However, standard LQ state feedback does not 
guarantee a zero gain reduction margin even in the open-loop stable 
case, amd has been criticized on these grounds [20]. Having a zero or 
small gain reduction margin is important in situations where actuators 
may fail or saturate, respectively, and there is no opportunity to re- 
configure the control system. In fact, the motivation for the thesis 
[26] (which in turn lead to most of the robustness developments reported 
in this chapter) was a study supporting the design of the automatic depth- 
keeping controller for the Trident submarine, in which sat\iration of one 
of the two hydrodynamic control surfaces produced an unstable closed- 
loop system. 

Corollary 4 ; For G(s) as in Theorem 3 the crossfeed tolerance is given 
by 

2 r Xmin(R^) Xmin(R )”j 

^ [ Xixxwp ■■ ' W(RjTJ ' ^ ® "r '™’ 

where L(s) is given by (55) , R > 0 is given by (57) and $ (s) has no 

UXi 

CRHP zeros. 

Proof : Analogous to Corollary 2. 

The significance of Corollary 4 is that using the Lyapvinov equation 
(64) to design the state teedback the toleramce to crossfeed perturbations 
has doubled over the crossfeed tolerance of the LQ state feedback regulator 
in Corollary 2. However, (64) can only be ustd on open- loop stable systems 
and thus even though the guaranteed stability margins for the Lyapunov 
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state feedback regulator are better than for the LQ regulator its use is 
limited. What is necessary is a compromise between these design approaches 
that is applicable to unstable open- loop systents. 

Another way to modify the LQ design procedure that is a compromise 
between Theorem 2 and Theorem 3 appliced>le to unsted>le open-loop systems 
involves the use of a parameterized Riccati equation given by 

A*^K + KA + Q - 6KBR"^b\ * 0 (71) 

where B is an adjustable parameter and 0 < B < 2. The feedback law is 
still given by (3) and G(s) is still given by (24) with K > 0. Since the 
B in (71) may be Jumped together with the R matrix, (71) is just a standard 
Riccati equation and therefore has a unique solution K ^ 0 under the 
appropriate assumptions (condition 2 of Theorem 2) . The st^mdard liQ 
optimal feedback law associated with (71) is given by 

u(t) = -BR"^B'^Kx(t) . (72) 

Instead of (72) we will use a(t) = -R ^^Kx(t) as in (3). Thus depending 
on whether 6 > 1 or 6 < 1 we are merely decreasing or increasing, re- 
spectively, the optimal feedback gain by a scalar factor of 1/B. Also 
with G(s) given by (24) the standard LQ loop transfer matrix is simply 
BG(s). From Theorem 1 we know that if Q > 0 

[l+BG(s)]” I R(1+6 G(s)) > I R , s e Ur (73) 

which in the SISO case becomes 

i I + g(s)l ' - . s e (74) 

ard is illustrated in Fig. 3. To obtain bounds on L(s) to ensure stability 
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we merely work with % Ha) and &G(a) and apply Theorem 2 for the standard 

p 

m regulator problem. Doing this we obtain, in the SISO case, the in- 


equality 


l6r^(s)-i| <1, s e 


illustrated in Fig. 7. Note that from that the critical (-1, 0) point 
is no longer contained inside the circle of Fig. 3 corresponding to (74) 
if B^2 and thus there are no giiaremteed margins. If the guaranteed 
minimum margins approach those of the Lyapunov feedbac)c case given 
Corollaries 3 and 4. In general, for the multivariable case the 


Im/(s 


Fig. 7. Set of allow 2 d}le values of il(s) when 
|6i!-“^(s)-ll < 1 and 0 < 6 < 2. 
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guaranteed ainunviR margins, again if R is diagonal and Q is positive 
definite, are given by 

C3M D [B/2, ») , 0 < 3 < 2 (76) 

and 

PM D [-cos"^ I , cos"^ I ) , 0 < 3 < 2 . (77) 

These guaranteed margins (when B £ 1) can also be obtained by similar 
but distinctly different procedures reported in (30J and [31] which utilize 
standard LQ regulators with v^mishingly small control weights. Recently 
[45] it has also been shown how to ensure preselected guaranteed minimum 
gain and phase margins by using a Riccati equation with an associated 
quadratic cost index, weighting the product of tlie state and the control. 

5.5 Stability Margins of LgG Regulators 

A basic limitation associated with the guaranteed stability 
margins is that they are obtained only under the assunption of full state 
feedback. State feedback can never be exactly realized, and often it is 
impossible or too expensive to provide enough sensors to achieve even an 
approximate realization. Thus one is motivated to investigate what 
guarantoevi stability margins might be associated with LQG controllers, 
in which a Kalman filter (KF) is used to provide state estimates for 
feedback , 

Since the Kalman filter is the dual of the LQ regulator, dual ro- 
bustness results are obtainable. They ensure a nondivergent Kalman filter 
under variations in the nominal model parameters of the plant whose state 
is to be estimated (see section 5.2 for the use of the Kalman filter in 
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the liQu regulator) . To make precise the connection bet%feen the regulator 
2 md filter problems, cmsider the linear system 


x(t) - Ax(t) 

+ £(t) 

(78) 

y(t) ■ Cx(t) 

+ 9(t) 

(79) 


where ^(t) and Q(t) ^u:e zero mean white noise sources with ^>ectral in- 
tensity matrices E and 6 respectively. We wish to estisiate x(t) given 
i(T) , 1 s'Jch that the meem square error is minimize. Under the 

assumption the [A,C] is detectable, it is well-known that the state esti- 
mate is specified by 

x(t) « Ax(t) + I C^G"^v(t) (80) 

v(t) - y(t) - Cx(t) (81) 

where 

AT + Ia'^ + 5 -Ec’b'^CZ - 0, Z > 0 . (82) 

If we calculate the tremsfer matrix from v(s) ♦■o £(s) - Cx(s) , we find 
that 

y(s) » (C(Is-A)"^Zc'^0"^lV(s) -F(s)y(s). (83) 

Then, if E > 0, F(s) satisfies the dual of (26) given by 
(I + F(s))9 (I+F(s))” >0, ® ® 

which guarantees the stability of the error dynamics under a range of 
perturbations in F(s). Thus, if F(s) is perturbed to F(s) - F(s)L(s) , 
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%diere usual aaauntptions about G(s) ara applied to F(s), the Kalman filter 
will remain nondiver 9 ent if 

- 1 ) < 1 , ■ ® ^ 

or equivalently, 

0l”(b) + L(a)0 - 0 > 0 . (86) 

It is now readily apparent that F(s), the lo <9 transfer matrix of the 
error dynamics loop of the Kalman filter, is the dual of G(s) in the 
IQ regulator emd has the same guaremteed margins at its input, v(s) , for 
diagonal G. 

Safonov and Athans [32] have developed these dual results for the 
nondivergence of the extended Kalman filter. Furthermore, they have con- 
sidered the stability properties of a iwnlinaar LQG control syston formed 
by the casceide of a constant gain extended Kalm^m filter and the LQ state 
feedback gain. The IQ state feedback gain and the constant gain of the 
extended Kalmar, filter are computed from the linearized model parameters. 
However, the extended 'Caiman filter must have the true nonlinear model 
of the plant. This violates the basic premise of robustness theory, that 
is, the controller has no knowledge or at most minimal knowledge of the 
model error. Nevertheless, the result snphasizes that model mismatch 
and not control or filter gains are responsible for a reduction in the 
margin of stability. We next examine these results in the completely 
linear case where the TiQG stability margins are much easier to obtain. 

The standard IQG control system block diagram is shown in Fig. 8 . 
with various points of the lo<^ marked. To determine the robustness of 
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th« LQG control aystw we ineert perturbation* at point* ai»l (the 
input and output of the ^yaical plant) and find out how large they can be 
mad* without deetabilizing the clo*ed-lo <9 LQG eyatem. It i* therefor* 
convenient to calculate the loc^^ tranafer matrice* at point* to • 

The loop tranafer matrix at point will de noted by *nd i* 

calculated by breaking the loc^ at point (aee aection 5.2) and uaing 

it a* the input a* trell aa the output. For the four point* indicated in 
Fig. 6 we have 


Tj (s) ■ G^4>{s)B 


TjCs) - C^(0"^(s) + BG^ + G.C)'‘^G^C^(*)B 


T^{») - C<i(s)BG^((>"^(s) + BG^ + GjC)'^G^ 


T^(e) - C4»(s)G^ 


(87) 

( 88 ) 

(89) 

(90) 


where 


A -IT 

G^ ■ R B K ■ regulator gain 


G^ - lC 0 
ipis) « (ls~A) ^ 


filter gain 


(91) 

(92) 

(93) 


Note that points at have the standard I^> regulator and 

Kalman filter loop transfer matrices respectively given previously in 
(20) and by F(s) in (83). Thus at points and (inside the 

LQG controller) the LQ and KF minimisB guaranteed stability margins apply. 
The fo) lowing theorem is a much simplified version of a theorem proved 
in [32] and gives LQG stability margins at points and (the 
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input and output of the physical plant) , where we denote the <^n-loc^ 
plcmt as = C(Is-A) ^ amd G^(s) = G^(Is-A + BG^ ■¥ G^C) ^G^ as the 

compensator. We use N(s) to denote output perturbations in G^(s) at 


point i 


in Fig. 8. 


Theorem 4 ; The LQG feedback control system of Fig. 8 is asymptotically 
stable under veuriations in the open- loop plamt G^Cs) = C(Is-A) ^ if the 
following conditions hold: 

(a) the perturbed open- loop plant G (s) = C(Is-A) is such (94) 

P 

tliat the det(sl-A) and det(sI-A) have the same nximber of 
CRHP zeros eind if det(jO)Ql-A) = 0 then det(jO)Ql-A) » 0. 

(b) [A,B] is stabilizable, Q > 0, k > 0 and K ^ 0 satisfies (95) 
(22) and B haa full rank. 

(c) G (s) = G (s)L(s) = N(s)G (s) (96) 

P P P 

and either 

0 (R*l"^(s)r"^-I) < 1 (97) 

max — 

or 

a (e"V^(s)0*-I) < 1 (98) 

msx — 

hold for all s Q . 

(d) the LQG controller tramsfer matrix G^(s) from the plant 
output to the plant input is given by 

G (s) = G (Is-A + BG + G,C)“^G, (99) 

c r r r f 

where and G^ respectively satisfy (91) 


and (92). 
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Proof: Breaking the loop at point of Fig. 8 we have a loop transfer 

function matrix of 

G^CsI-A+G^C)“^[G^C(Is-A)"^B + B] - G^(sl-A)‘^ - G{s) (100) 

so that 

4>ol(s) = det(sI-A + G^C]dettsI-A] (101) 

and 

^Oj^(s) = det[sI-A + G^C] dettsI-A] (102) 

Since the Kalmar filtering error dynamics are stable given (95) and since 
(94) holds, conditions la and lb of Theorem 3.2 hold. Now by direct appli- 
cation of Theoran 2 we conclude that the system of Fig. 8 is stable if L(s) 
is inserted at point • However, this is not the location we desire to 
have the margins guaranteed. Nevertheless, by manipulation of the block 

diagram of Fig. 8 we may place L(s) at point if we change B to BL(s) 

-1 T 

inside the controller leaving - R B K fixed (see Fig. 9) . This, 

however, is equivalent to chemging (A,B,C) to (A,B,C) inside the controller 

leaving G^ and G^ fixed. This also can be interpreted as giving the 

Kalman filter the correct dynamic model of the perturbed open-loop system 

without changing either the filter or regulator gains. The same result 

follows if we start with the perturbation N(s) at point ^^where the KF 

guarantees apply and move it to point chamging C to N(s)C (see 

Fig. 10.) Again the Kalman filter has the correct model of G (s). 

P 

Q.E.D. 















Control System and perturbation 
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Notice that in (96) L(s) represents the sane perturbation in G^(s) 

at the input to the pl 2 uit as N(s) represents at the output of the plant 

and that G (s) is the same in both cases. 

P 

Thus the LQ and KF guarcmteed stability margins will apply to UQG 
controllers at the input euid output of the physical plant but under the 
restrictive assumption that the system model embedded within the Kalman 
filter is always the same as the true system (i.e., the perturbed syst^) . 
For the more realistic case in which the internal model of the Kalman 
filter remains unchanged, there are unfortunately ^ guaranteed robustness 
properties^, as Doyle has demonstrated with a simple counterexample [33]. 
This coxinterexample is extreme, but it is possible to obtain LQG con- 
trollers with inadequate stability margins that look quits reasonable 
in the time dcsnain. Fig. 11 shows the Nyquist plot of a single-input 
design reported in the literature [34]; note that the phase margin is 
less than 10°. 

5.5.1 Robustness Recovery 

Fortunately, there are two dual procedures that do not require 

the Kalman filter to have the true system model and that still recover 

the IiQ and KF guaranteed minimxmi margins. These procedures use the 

asymptotic properties of the Kalman filter and LQ regulator (see [43] 

and can be used only if the plant is minimum phase. If W is a non- 

T T 

singular arbitrary matrix, then by selecting = in (82) as p BWW B 
^md letting p ■> “> the loop transfer matrix T 2 (s) in (88) approaches 

^In other words, the robustness properties of LQG designs will depend 
on the actual values of A,B,C,Q,R,H and 0. 


Re H(j(u) 


Fig. 11: Nyquist diagram for LQG design in [34] 

(H(jw) = loop transfer function). 

T^(s) of (87) if the minimum phase assumption holds [35]. Thus the 

LQ regulator guaranteed margins will be recovered at the input to the 

plant. KwaJcemaak [36] proposeed the dual of the above procedure to 

obtain low sensitivity feedback systems. His procedure makes T^Cs) 

T T 

of (89) approach T. (s) of (90) by selecting Q in (4) as p C W WC and 
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letting emd thus the KF guaranteed minim\sn margins will be recovered 
at the output of the plant^. However, it is not always the case that an 
LQG controller needs to be robustified by these procedures since in sone 
cases the LQG control system will have better stability margins than its 
full state feedback counterpart [42] . Also, in this case one must keep 
in mind that the stochastic error performance of the robustified liQG 
controller may be better (due to lower controller bandwidth) them the 
liQ state feedback regulator. 

Even when these procedures are uced, the guaranteed stability 
margins apply at the input ot output of the physical plant but not 
necessarily at both input and output. It is desircdsle to have margins 
at both these locations since the perturbations in Gp(s) are repre- 
sented as either G (s)L(s) or N(s)G (s) and we would not like small 
P P 

perturbations in either input or output to destabilize the system. 
Margins at both input and output can be ensured if the inequalities 

a . (I+G (s)G (s)) >1 (103) 

min c p — 

cmd 

o . (I+G (s)G (s)) >1 (104) 

min p c — 

both hold. The relationship between Utese two quantities when ^^(s) 
and G^(s) are square matrices is given by the following lemma. 


^owdle [37] has adapted these procedures for 
observer based compensators and their duals. 


use with minimal order 
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Lenma 1: For arbitrary complex square matrices G_(s) and G_(s) it is 

P ® 

true that 


Z O , (I4G (s)G (s)) < o . (I+G (s)G (s)) < ko , (l+G„(s)G (s)) 
k min p c — -min c p — min p c 


(105) 


where 


mrn 


0 (G (s)) 
m^uc p 


o (G (s)) 
max c 


o . (G (s)) 
min p 


a . (G (s)) 
min c 


> 1 . 


(106) 


Proof ; Use the property of matrix norms that ||ab|| < ||a|| ||b|| on the 
equation 

[I+G (s)G (s)l~^ = g‘^(s)(I + G (s)G (S)]"^G (s) (107) 

pc p c p p 

to obtain the left inequality of (105) with k = o [G (s))/0 . [G (s)]. 

^ ^ max p min p 

The right inequality in (105) is obtained by reversing the roles of 

G (s) and G (s) in (107). 

P c 

Q • E • D • 

The quantity k is the minimvm of the condition nunbers^ of 0^(3) 

and G^(s) with respect to inversion. From (105) we conclxxle that if k 

is close to xmity then approximately the same robustness guarantees will 

apply at both input and output. Note that we have no control over 

Gp(s) so that if Gp(s) is nearly singular we must design our compensator 

so that 0 (G (s)) =0 (G (s)). On the other hand, if our plant is 
max c max c 

well-conditioned with respect to inversion, our compensator G^(s) need 


not be so severely constrained, allowing more flexibility in achieving 


In the msnerical analysis of the linear equation Ax * b, the condition 

minber of A, given by o (A) /o . (A), bounds the error in the computed 

max min 

solution X in terms of an equivalent error in b [44]. 


perfonnance objectives 


5.5.2 Characterization of Model Error 

Note in (103) and (104) that ensuring that both o , [I-fG (s)G (■)] 

min p c 

and greater than unity gives upper bounds on two 

different types of allowable modelling errors. These modelling errors 

differ in that one represents the perturbed open* loop model ^^(s) as 

G (s)L(s) , cm equivalent perturbation in the input to G (s) , and the 
P P 

other represents 0^(s) as N(s)G^(s) , an equivalent perturbation in the 

output of 0^(s) . Those are both relative modelling errors between 

G ^(s) and G ^(s) and are given by 
P P 


E,(s) - l“^(s)-I - [g'^s)-g‘^(s)]G^(s) 
1 p p p 

E„(s) - n‘^(s)-I - G (s) ( g“^(s)-g‘^(s)1 
0 P P P 


( 10 ( 

( 10 < 


where E. (s) is the model error in G (s) reflected to the input of G (s) 
1 p p 


and Eq(s) is the model error in Gp(s) reflected to the output. The re- 


lationship between E^(s) and ^^(s) is given very simply by (since they 


represent the same G^(s)) 


,-l, 


E. (s) - G (s)E.(s)G (s) 

1 pop 


(iio: 


and thus we conclude that 


1 

k 


a (E„(s)) <0 (E.(s)) < k 0 (E-(s)) 

max 0 — max i — max 0 


(111 


where 


k 


(G (s)l 
max _£_ 


0 , [G (s)l 


(1121 
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If k is very large, (111) shcMS that what may loo)c like an unreasonably 
large error at one point in the loop (either the input to G^(s) or its 
output) may look like a very reasonable size error at anotlier point in 
the loop. The stability of the closed- loop system with respect to 
t)iese errors is guar£u)teed under appropriate assumptions^ if for all 

.en„ 

o [E.(S)I < a . n*C |s)G (s)l (113) 

max 1 min c p 

or 

a [E^(s)l < 0 . [I+G (s)G (s)] . (114) 

max 0 min p c 

Now suppose that a [E. (s)], from our knowledge of the open-loop system 
max 1 

physics, seems unreasonably large and condition (213) does not hold 

but will hold for all reasonable size errors. It the perturbed model 

G (s) however is such that 0 [E. (s)l seems reasonably small, that is 

p max 0 

an error of that magnitude could be justified (again by our )cnowledge 

of the open-loop system physics) , then a sufficiently small value of 

o , [I+G (s)G (s) 1 (small enough so that (114) does not hold) indicates 
mxn p c 

that there is danger of the perturbed closed- loop becoming unstable. 

If we believe that having a reasonable size error, at the input of 
the open-loop systan, completely characterizes the class of Gp(s) 
that should be considered then we may rule out the possibility that the 
perturbed closed- loop system may become unstable as a result of the 
fact that (114) is violated for seemingly reasonably sized errors at 

^Conditions I and 2 of Theorem 3.6 must be appropriately modified to 
also deal with N(s) . 
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the output. However^ our knowlodae of t«hat class of G (s) should be 

P 

considered may depend on what se«n like reasonable errors at both input 

and output to the open-loop syst«e. That is, our knowledge of what 

constitutes a reasonable class of G (s) to consider is built up frosi 

P 

our knowledge of what errors seen reasonable whsn reflected as eqpiivalent 

perturbations at different points in the feedback loop. Zn this cane, 

one could not rule out the small error at the output (i.e. , which 

violates (114) . Thus it depends on how we decide what constitutes a 

reasonable class of G (s) that determines if ve need to c^eck to isake 

P 

sure that botli the values o , [I-t-G (s)G (s)] emd o . (Z-fG (s)G (s)] are 

min c p min p c 

sufficiently large or that only one of them is sufficiently large. 

In the author's opinion, it sems most likely that the class of 

G^(s) that should be considered is a composite class of tlws perturbed 

models that arise from reasonably sized errors reflected to both input 

cind output of the open-loop system. Therefore, it would seon wise to 

-heck the size of both a . tl+G (s)G (s)] and o . [I+G (s)G (s)]. 

min c p min p c 

5.6 Concluding Remarks 

This chapter has derived MIMO stability margins for LQ regulators 
and their variations including LQG regulators. This was accomplished 
using th* MINO version of Kalman's inequality and Theorem 3.6. The LQ 
regulator was shown to have at least a 50% gain reduction margin, an 
infinite upward gain margin and t 60° phase margin provided the control 
weighting matrix R is diagonal. If R is not diagonal it was shown that 
the LQ regulator gain margin may be arbitrarily small. If R is not a 
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multiple of the identity the crossfeed tolerance is also reduced. Ihe 
R matrix determines the coordinate systm in which the stability suirgins 
hold. The margins in the R selected coordinate system stay be much larger 
than the actual margins in the coordinate system specified by the inputs 
and outputs to the i^ysical open-loop system. Similar c^msents may be 
made for the variations of LQ state feedback using the Lyi^unov and 
modified Riccati equations. 

The guaranteed maurgins for LQ regulators do not a{^ly to LQG regu- 
lators except when the Kalman filter embedded in the LQG controller has 
a correct dynamic model of the perturbed system, a rather tmrealistic 
assvnption. However, when the open-loop plwt model is minimum phase, 
there are two procedures that recover the guauramteed IQ margins 
asymptotically. These 'Tuaranteed maurgins may be recovered at either 
the input or the output of the open-lo<^ plant but can only be guaranteed 
to be recovered at both Input and output when either the open- loop 
plant tramsfer matrix or the compensator transfer matrix has a small 
condition number near unity for all frequencies. The necessity of input 
and output stability margins is shown to be dependent cn the ability 
of the designer to characterize the set of reasonaible perturbed models 
for which the perturbed closed-loop system stability must be preserved. 

It is important to point out that the liQG methodology is inherently 
a multiloop design procedure which when coupled with the robustness 
recovery methods emd used intelligently provides a systematic design 
procedure for robust multivariable ocmipensators . This is in contrast 


to the characteristic loci (4, 5, 56] and inverse Nyquist array [1,2] 
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nethods (dlscusMd in the folloning chapter) which reduce multivariable 
controller design to a series of decoupled single loop designs. Thee# 
methods obtain good stability siargins in the coordinate systoi of the 
decoupled STSO systems but not in the coordinate system of the physical 
input and output of the open-loop multivariable syst«a. Thus they may 
not produce robust controller designs and it is nontrivial, given the 
present state of the art to change these designs to obtain better r^mst- 
ness properties. 

One word cf caution is necessary when using the LQG aj^roach to 
controller synthesis. It has been popular to reduce many more general 
control problems such as the trac)ci.ig problem to a simple regulator 
problem by state augn^ntation. State augmentation has also been used 
to provide integral controllers or provide additional rolloff. When 
using augmented versions of control probl^s solved via the regulator 
problem one must be careful to determine exactly the point in the 
feedback loop whore the guaranteed stability margins will apply. In many 

augmented regulator probl^s the point at which the guaranteed margins 
apply is not the input or output of the physical open-loop plant but a 
point associated with the £iddition of the augmented states. 
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5. BOBUSTHESS MU^YSIS WITH FBfiQUEMCY DOHMN METHODS 

Introdwtion 

Th« purpose of this chapter is to place in perspective current 
freqiKincy domain techniques for controller synthesis and evaluation and 
their implications for the robustness characterization of feedback 
control systems. He shall not present a full tutorial description of 
these methods but only briefly describe their salient features with 
regard to their importance in robustness analysis » the main thmne of 
this thesis. 

In section 6.2 the characteristic loci (CL) (5,561 *nd inverse 
Nyquist array (INA) [ 1 , 2 ] methodologies are discussed. It is shown 
that these design methodologies ensure stability margins in a coordinate 
system based on the diaganalization of the open- loop plant transfer 
matrix rather than the coordinate system specified by the physical in- 
puts and outputs of the ncaninal open- loop plejit. In some cases, the CL 
and INA design methodologies will lead tc< acceptable stability margins 
at the physical input and output of the systmn; however, in other cases, 
the stability margins at the physical input and output of the system will 
be drastically reduced. The discussion of this section is not original 
but relies heavily on the work of Doyle and Stein [43] and Stein and 
Sandell (58] . 

In section 6.2, the principal gain emd phase (PGP) analysis recently 
proposed in [57] is discussed. This method of analysis allows one to 
ensure stability of a feedback control system by tedcing into account 
the structur of the imadel in a somewhat different manner than that of 
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Chapter 4. The main shortcoming of the PGP approach is that it is not 
aiqplicable to in^rtant classes of loop transfer matrices, G(s), amd 
important classes of model errors. It is shown that the PGP method 
will fail for model error matrices E(s) that are singular or almost 
singular. Recall that in Chapter 4 (equation (4.13)) it was shown 
there exists a smallest destabilizing error matrix, E(s), that is 
singular and thus the PGP analysis is not applicable for this important 
class of model errors. 

6.2 Characteristic Loci and Inverse Nyquist Array Methods 

The CL and INA methodologies for the design of HIMO feedback control 

systems teJce advantage of the Isurge body of well-developed tools for 

".iiO control design by reducing the MIMO design problem to a sequence 

of independent SISO design problems. To make this precise, consider 

the feedback system shown in Fig. 1 where the nominal open-loop pleuit 

transfer matrix denoted by 2^(3) and tlie compensator tramsfer matrix is 

denoted tr/ giving a loop transfer matrix G(s) as either G^(s)Gp(s) 

or G (s)G (s) depending on whether the loop is broken at the input or the 
P c 

outpur respectively. 


Compensator Plont 



Fig. 1; Basic Feedback System 
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Both the CL and INA methods assume thet G (s) can be diagonalized 

P 

either exactly or approximately for all s 6 D^. Therefore, assume that 
there exist matrices ”(s) and V(s) composed of rational tremsfer functions 
such that V(s)G (s)W(s) ic bounded and either diagonal or diagonally 


dominant 


P 

for all s e eind is denoted 

S (s) » V(s)G (s)W(s) . 

P P 


by 


ii (s) , i.e. 
P 


( 1 ) 


If G (s) is diagonally dominant then the diagonal matrix 
P 


e .(s) 

pd 


given 


by 


*pd'“> ' diagKjpjjis), 


( 2 ) 


can be used as a good approximation to 6^(s) within bounds specified by 
the magnitude of the off diagonal elements of 6^(s). For the purposes 
of this chapter we will assume exact diagonalization of by V(s) 

and W(s), since all the same observations to be made will apply if 
only diagonal dominance holds. 

The form of the compensator proposed by the CL and INA methods 
is shown in Fig. 2 where K(s) is a diagonal matrix given by 


An nxn complex matrix A is diagonally dominant if for i = l,2,...,n 


a.. > Z a-- 

j.l 


j/i 


or 


n 

ja. . I > Z la. .1 

jf^i 


Compensotor 

_ GdsJ__ 


C Innt 


U 



Fig. 2. Compensator used by INA and CL Methods 


K(s) = diag[k^(s), 


,k (s)l 

n 


(3) 


If the loop transfer function calculated at point i,i=l,2,3,4,is denoted 
by (s) , then 


G,(s) = (W(s)K(s^V(s))G (s) =G (s)G (s) (4) 

1 pep 


and 


G_(s) =K(s)(V(s)G (s)W(s)) = K(s)6 (s) (5) 

2 p p 

where we note that G^Cs) is diagonal. The CL and INA methods use (5) 
to design the compensator G^(s), by selecting eachk^(s) in K(s) as 
the appropriate robust compensator for each of the SISO systems repre- 


sented by the diagonal elements of G (s) , denoted g .. (s). Thus, these 

P PiJ- 
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methods produce feedback control systems with good margins inside the 
compensator at point in Pig. 2. 

The key question is; when does this design methdology yield good 
margins at point GO , the input to the physical system ? To determine 
the answer to this question, suppose that we insert the matrix Lj^(s) 
at point to account for model uncertainty and determine the 
stability margins at the input to the physical system by placing bounds 
on the allowable Lj^(s) . If we represent the equivalent nodel uncert6uinty 
at point by then the relationship between Lj^(s) and L^Cs) 

is given by 


= W~^(s)L^(s)W(s) . 


(6) 


In chapter 3, the model error criterion (3.35) was used and is 
given in terms of both L (s) or G. (s) and G. (s) by 

ill 

E. (s) = g"^(s) [G. (s) - G. (s)] = L. (s) - I . 
1111 1 

Thus, using f6) and (7) , the tolereible model error at point ©is 
related to that at point by 


(7) 


E^Cs) = w‘^(s)E^(s)W(s) . 


(8) 


From (8) and the properties of singular values, we obtain 


0 [E.,(s)] < c[W(s)]0 (E, (s)l (9) 

max 2 — max 1 

where c[W(s)] is the condition number^ of W(s) and is given by 


^The condition nimber of a matrix is very large if the matrix is nearly 
rank deficient or almost singular. 
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a (w(s)] 

C[W(S)] - > 1 . (10) 

Suppose that the model error E^Cs) at point of miniicm magnitude 

0 (E_(s)) that destabilizes tl^ feedback system is such that eq[uality 

max <£ 

holds in (9) . Then the equivalent model error Ej(s) at point 
determined from (8) is of magnitude — . r -y - v; o [E_(s)] and also de- 

C t W ( S / J RI3X <b 

stabilises the feedback system. Therefdre, if the condition number 
c[W(s)) is very large, the margins at point , i.e., at the input to 
the physical system, may be much smaller than these at point , inside 
the compensator. 

The CL method selects W(s) to be the matrix of eigenvectors of 

G (s) and V(s) = W ^(s) when this possible choice for W(s) is rational; 

P 

otherwise, W(s) is chosen to be as rational approximation of the 

matrix of eigenvectors of G (s) . Since, the designer has no control over 

P 

G (s) , this choice of W(s) does not guarcmtee that the condition manber 
P 

c[W(s)] is near unity. Similary, the IMA method seeks to find a rational 
W(s) emd V(s) that diagonalize C^(s), (this is done niznerically to 
obtain constant matrices W(s) and V(s)) but there is no guaramtee that 
c[H(s)] is near unity. Similar conclusions may be drawn by breaking 
the feedback loop of Fig. 2 at points and and working with V(s) 
rather than W(s). 

Therefore, these methods do not autcnatically produce robust con- 
troller designs. Indeed, in some cases they can lead to nearly unsted^le 
feedback systems (i.e., small stability margins at either input or out- 


put) if the diagonalizing matrices W(s) and V(s) have large condition 
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numbers at some Irequency. They do produce good margins at a point inside 
tl% compensator but that is not the appropriate place to require good 
margins, from an engineering point of view. Doyle and Stein [431 give 
a simple example which exr>oses this deficiency. 

6.3 The Principal Gain and Phase Analysis Method 

This method [57] utilizes model error structure information 
to ensure stability of a perttirbed feedback system. It uses the notions 
of principal gains , and principal phases of an nxn complex matrix A which 
are defined via its polar decomposition. 

Definition (Polar Decomposition) ; Any nxn complex matrix A cami be de- 
composed into a product given by 

A = (11) 

or 

A = H.U (12) 

Xj 

where U is a unitary matrix and, H„ and H are positive semidefinite 

K Ij 

hemitian matrices. The representation in (11) or (12) is called a 
polar decomposition of A. 

Note that the polar decomposition in (11) ^md (12) are easily calcu- 
lated from the singular value decomposition (SVD) of A as follows 

A = (UV) (v“sv) (13) 

and 

A = (UlU*^ (UV) 

where A = UEV ^ is the SVD of A. 


(14) 
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Definition (Principal Gains and Phases) ; The principal gains of the 
matrix A in (11) (or (12)) are the eigenvalues of (or 
identical to the singular values of A. The principal phases of A in 
(11) (or (12)) are the arguments of the eigenvalues of U in (11) and (12). 

Since the eigenvalues of a unitary matrix are of the form e , 
the identification of the 6^ as phases is obvious. 

Definition (Spread of less than tt) ; If the principal phases of a matrix, 

j^i 

denoted are such that the complex numbers e can all be contained 

strictly in a half-plane in the complex plane that has the origin of the 
complex plane on the boundary, then the principal phases are said to 
have a spread of less than it. 

This is illustrated in Figs. 3 and 4. 


Hj^) which are 


Im 



Fig 3: Principal phases with a spread of more than it. 
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The model error criterion to be used is given by 

E(s) - g‘^{s)[3(8) - G(s)l (16) 

where G(s) as usual represents the perturbed- loop transfer matrix. Let 

the principal gains of G (jw) be denoted as a. (w) where 0 < a. (w) < a. , (w) 

and let the principal phases of G_. ( joi) be denoted as 0.(o)) where 

CL 1 

£ ^i+l^^^‘ Similarly let the principal gains of E(jo)) be de- 
noted as 6^((i)) where 0 _< 6^(u) £ principal phases 

of E(ju) be denoted as e^(d)) where e^(u)) < 

Nert, define the condition nunbers Cj^(u)) and c^ (uj) as 


. a (w) 

C (0)) = — > 1 
1 a (u)) - 


(17) 


wltere is the minimum singular value (or principal gain) and is 
the maximum singular value (or principal gain) of G and 

Cli 




6 (u) 
m 


5 ^ (w) - 


> 1 


(18) 


where 6, is the minimum singular value and 6 is the maximun singular 
1 m 

value of E. Also defined the quantity (o)) as 

m 


\I) (u)) * tan 
m 


/(Cj^(U)) - 1] c^(ui) \ 

\ 1 - lc^(u))-l]c2(u)) j 


( 19 ) 


With these preliminary definitions the so-called Small Phase Theorem 
(SPT) of [57] may be stated. 


Small Phase Theorem ; The perturbed closed- loop system is stable 
if: 
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1. E(s) is stedjle. 

2. {0^(u))} amd {0^(u>) + C^(«)} have a spread less 
than IT for all u) and i,j ■ l,2,...,m 

3. (Cj^(u))- IJc^CiiJ) < 1 for all o) 

4. (a) 0 (oi) + e (u)) < IT - 1 ^; (w) for all u) 

tn in n 

emd 

(b) 0,(0)) + c, ( 0 )) > ( 0 )) - 71 for all o) 

1 1 in 

This theorem basically haracterizes the tolerable model errors as 
those that do not introduce a significant amount of phase shift in the 
principal phases of G (juJ) when perturbed by model error (condition 4). 
Conditions 2 and 3 place restrictions on the type of system and the type 
of model error that can be considered by the SPT. Condition 1 is simply 
a condition that automatically guarantees that the matrices G(s) and 
5(s) have the same number of unstable poles. This theorem, as those 
presented in Chapter 4 of this thesis uses the structure of the error 
matrix E(jw). However, the SPT does this by requiring restrictions on 
its principal phases c^(u>) in conditions 2 and 4, and, therefore is 
rather different from the characterization of model error in Theorems 
4.1 and 4.2 in terms of the projection of th error matrix E(jui) onto 
various one dimensional subspaces generated by the singular vectors of 
I+G"^(ju)) . 

The main drawback of the SPT is that it cannot be applied in m£my 
cases of interest. This is illustrated by the restrictions condition 
3 places on the system and the model error in the following two single 


examples. 
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Exanple 1: The SPT does not apply when the closed-loop transfer matrix is 

given by 


Gcl< 8) - diag 


10 1 
s+10 ' 8+1 


because the condition nxmaber Cj^(u)) defined by (17) is 


(20) 


(co) 


100 (u+l) 




to^+lOO 


2 for a) > 1.1 


( 21 ) 


and th\is [Cj^(o)) - l]c 2 ((i)) > 1 for u) > 1.1 since ^ 1 for any matrix 

E(s). Hence, condition 3 of the SPT is violated. The loop transfer matrix 


G(s) corresponding to G (s) in (20) is given by 

CL 


G(s) s diagt ^ • (22) 

Example 1 demonstrates that closed-loop syst^s with input-output 
channels with widely differing bandwidths violate condition 3 of the SPT 
because Cj^(u)) > 2 for some frequency u.. In (20) the first input-output 
channel with transfer function 10/ (s+10) has a bandwidth approximately 
10 times as large as that of the second input-output channel with 
transfer function l/(s+l) . If Cj^ (u)) is to be less than 2 for all u, 
it is necessary that all input-output channels of G (s) have roughly the same 
bandwidth or equivalently the same speed of response. This same restriction 
also applies to G(s). Clearly, this restrictive condition eliminates many 
systems of interest from the point of view of robustness analysis via 
the SPT. 




Example 2 : Suppose I+G (jw) has the SVD given by 

I+G ^(jto) * U( ju))E (ju))V**(iu)) « E 0 . ( jui)u. ( ju))v!*(juj) (23) 

i ■~i ~~t 


i-1 


-247- 


where o (ju) - c . < jco) . If the taodel error E(ju) is such that for sane 
n sun 

w_, for which o (jw^) < o (jw) for all w , 
u TO u “• ro 

E(jU)o) - 

then from chapter 4 (equation 4.13) we )cnow that the closed- loop feedback 

syetem is destabilized by the model error E(s) but is stable for all tsodel 

errors E^Cju)) of magnitude o tE.(ju))] < o [E(ju> )]. That is, E(jw*) 
u max u Bkax u u 

is the smallest destabilizing model error. Hov»ever, condition 3 of SPT 
requires that 

0 . (E(jw)) “ c ( 0 ))-l 

min 1 

and thus since E^jw^) is singular » +® and condition (25) is violated 

and thus also condition 3 of the SPT. Again, the SPT does not apply. 

The results of the above example are significant because in a useful 
robustness analysis method it is very important to detect the smallest 
possible errors that may destabilize the feedback system and be ed^le to 
distinguish them frcxn model errors of equal magnitude that do not de- 
stabilize the feedback system. In many systems the smallest possible 
destabilizing errors may lead to a singular error matrix E(jo!) at 
S(xne frequency. Note that even when E(ju)) is not singular, the inequality 
(25) is still easily violated if Cj^(ui) f I (i.e., all feedback loops do 
not have the same bandwidth as illustrated by Example 1) . 

These two examples have shown the SPT is a fairly restrictive 
theorem viewed from the robustness viewpoint. However, in certain cases, 
the SPT may be an easy way to most sin^jly characterize a particular class 


of model errors. 
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6.4 Concluding R«aark» 

This chapter haa briefly exanined the characteristic loci (CL) 
and inverse Nyquist array (XNX) control design procedures and shown 
that they do not guarantee good robustness properties of the resulting 
control system. The essential deficiency in these siethods lies in the 
fact that they can ixily guarantee the desired stability margins at a 
point InsMe the compensator. The true stability margins at the inter- 
face between the physical plant and i*he compensator can be bounded 
in terms of the stability margins insid? the compensator. This bound 
depends on the condition numbers of the matrices used to diagonalize 
the open-loop plant during the controller design procedure. If the 
condition nvmibers of these matrices are large, then the stability 
margins at the interface of the physical open-loop system and the com- 
pensator may be exceedingly small. When this happens, there is no 
presently known methocl to modify or correct the CL and/or INA designs 
based on information implicit in the diagonalizing matrices that have 
l 2 urge condition msnbers in order to obtain the desired robustness 
properties. Current reseeurch efforts^ are being directed toward 

developing methods of diagonalizing G (s) with raticmal trai.sfer 

P 

function matrices that are unitary for s^jo). This would guarantee that 
the margins th&t hold inside the compensator would also hold at the 
input or output to the physical system. 

In contrast, to these aforementioned difficulties the LQG design 
procedure, utilizing tliC robustness recovery techniques discussed in 

^Private conmiunication with Professor Bernard Levy. 
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Chapter 5, docs provide a systematic procedure to synthesise controllers 
that are robust or can be made robust if the nominal desi9n model is 
sdnimxm phase. Thus, if an LQG controller is evaluated for robustness 
and found lackir.g, there is a systematic way to restore the desired 
robustness properties and not merely trial and error in the redesign of 
the controller in «m ad hoc manner. A drawback of LQG controllers is 
the large dimension of compensator states, often larger than necessary 
to meet performance and robustness objectives. However, in this era 
of microprocessor and VLSI impl^entation of control compensators the 
dimensionality of the cconpensator is not as crucial a problem, as was 
the case a few years ago. Current research [59] in the design of reduced 
order robust liQG controllers is progressing. 

The other recently proposed principal gain and phase (PGP) fre- 
quency domain euialysis discussed in this chapter ensures stability 
of a feedback system in the face of a special class of model uncertainty 
based on the model error matrix structure exhibited by its principal 
phases. It was shown however that this approach is not applicable to 
systems whose closed-loop speed of response in different input-output 
channels differ significant!-, . Furthermore, it was shown that PGP 
analysis requites a nonsingular error matrix. However, in Chapter 4 
it was shown the smallest destabilizing model error matrix may be taken 
to be singular. Thus PQr‘ analysis, cannot be applied to determine 
if a singular model error matrix is stabilizing or destabilizing. 

However, it must be said that it offers the potential in the case of 
model errors known to be of a small principal phase nature, to provide 
a relatively uncomplicated test for feedback system stability. 
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7. SUMMARY, CONCLUSIONS AND SUGGESTIONS 
FOR FUTURE RESEARCH 


7.1 Summary 

This thesis has addressed the following problaa. Given a finite- 
dimensional linear-time invariant feedback control system designed using 
an inaccurate ncxtiinal model of the open- loop pleuit, how much and what kind 
of model error can the feedback system tolerate without becoming unstable? 
Thus, this thesis deals primarily with the evaluation of the robustness 
of stabi lity of a feedback control system. This rol^stness evaluation is 
absolutely essential since all models of physical processes are only 
ai)proximations to the actual relationship between the system inputs and 
outputs. In the single- input, single-output (SISO) case, this evaluation 
is readily accomplished using frequency domain plots, (e.g. using a Bode 
diagraiTi) to display the behavior and characteristics of the feedback 
system. However, in the multiple-input, multiple-output (MltC) case, 
many generalizations of the SISO methods have proved inadec[uate because 
they have not dealt with the MIMO system as a whole but as a sequence of 
SISO systems. 

This thesis has avoided this deficiency by utilizing standard 
matrix theory concepts and methods appropriate for dealing with the MIMO 
case, ncmely the singular value decomposition (SVD) and properties of 
special types of matrices. These were discussed in Chapter 2, where the 
main problem solved was the determination of the nearest singular 
matrix. A, to a given nonsingular matrix. A, under certain constraints 
on i -A. The solution to this problem (given in Problems A, B and C) 
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is fundamental to the control system robustness results of Chapters 3 and 

4 . 

The basic foimulation of the control system robustness problem 
was considered in Chapter 3 via a multivariable version of Nyquist' s 
stability theorem. There, a fimdamental robustness theorem (Theorem 
3.2) was derived that implicitly characterized the class of perturbed 
models that would not destabilize the control system, in terms of the 
nonsingularity of the return difference matrix. Various robustness 
tests (Theorems 3.3 to 3.9), were then derived which can be used to 
test the nonsingularity of the return difference matrix for several 
types of model error criteria. These results were then related to the 
small gain theorem and some simple extensions for nonlinear feed:v>ack 
control system were presented that demonstrate that the basic rc.'bustness 
results of Chapter 3 are valid even when certain types of nonlinearitics 
are introductHi. 

Chapter 4 heavily utilizes the results of chapter 2 in determining 
what types of model error will destablize a given feedbac)c system. Model 
errors that tend to destabilize the feedback system are distinguished 
from those that tend to stabilize the feedback system by examining their 
structure as well as their magnitude. The key results, contained in 
Theorems 4.1, 4.2 and 4.3, show that the magnitude of the model errcr 
necessary to destabilize the feedback system may greatly increase if the 
class of model errors that can plausibly occur does not include model 
errors that are essentially alike in structure to the model error of 
minimum size that will destabilize the feedback system. This provides 
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an important characterisation of the model errors that are important in 
feedback desion. These types of model errors were then interpreted via 
block diagrams of the feedback system. 

In Chapter 5, the robustness properties of control system designed 
using the linear-quadratic-gaussian (LQG) methodology were presented. 
Multiloop guaranteed gain, phase and crossfeed margins were obtained using 
the robustness theory of Chapter 3. Tt.e key results of Chapter 5 are 
contained in Theorems 5.2, 5.3 and 5.4 and their corollaries. These 
theorans shewed that the guaranteed LQ robustness properties hold in a 
coordinate frame defined by the control weighting matrix (R) in the 
quadratic performance index. It was also shown that LQG control systems 
cannot automatically guarantee the same stability margins as the LQ state 
feedback regulator unless tl’ ' in' ".nal model of the system embedded in the 
compensator Kalman filter is cor ’ c. This is a very restrictive condition 
and therefore robustness recovery procedures were outlined that do not 
require the compensator to have exact knowledge of the correct dyncimic 
model of the system. These procedures allow LQG controllers to recover 
the robustness properties of LQ state -feedback controllers if the nominal 
open-loop plant is minimum phase. 

Chapter 6 contrasts the frequency domain techniques for MIMO analysis 
and design (characteristic loci (CL) , inverse-Nyquist array (INA) and 
principal gain and phase (PGP)), to the methods taken in this thesis and 
demonstrates that the CL and INA design methdologies do not ensure robust 
controller designs and that the PGP analysis is not able to determine the 
robustness of important classes of systems and model errors. 
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7.2 Conclusions and Suggestions for Future Research 

A systematic procedure for the design of robust feedback control 
systems is the ultimate goal of any analysis of tlte robustness properties 
of feedback systems. The impact of the analysis techniques developed 
in this thesis upon controller design is obviously not that of a new 
method for obtaining a controller given a model of the plant. Rather, 
this thesis has dealt primarily with the characterization of model un- 
certainty and how it affects current design methodologies. 

The theoretical impact of this analysis consists of the following 
developments : 

• the formulation of a new robustness theorem that exposes 
the fundamental character of all robustness tests which 
can be derived from it. 

• the unification of various robustness tests by the classi- 
fication of the type of model error they bound 

• the derivat n of new robustness tests based on alternate 
types of model error criteria not previously considered 
in the literature 

• the fundamentally new characterization of model error, 
which requires only a partial knowledge of the model 
error, based on its projections onto certain subspaces. 

The design impact of this analysis is not as clearly defined and 
future reseai'ch to develop design techniques is necessary. However, 
the design developments that seem possible are: 

• controller redesign based on the model error structural 
characterization of model errors critical to stability 
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• model improvement to reduce model errors critical to 
stability 

The controller redesign and model improvement go hand in hand. Once 
an initial design has been accott^lished and the model errors critical 
to system ocai. iJity are identified, these model errors can be incorporated 
into a new model and a new controller design, using this new model, may 
be performed explicitly taking into account the characteristics of the 
new model critical to system stability. 

Other directions in future research that might prove productive are 
tied to specific design methodologies. For LQG controllers, the develop- 
ment of procedures for producing robust low order compensators from robust 
high order compensators would be highly desirable. To produce a robust 
reduced order compensator, the loop transfer matrices corresponding to 
the high order and reduced order ccanpensators must be nearly alike. More 
precisely, if Gj^(s) is the loop transfer matrix resulting from the use 
of the high order ccffnpensator and is the loop transfer matrix re- 

sulting from the use of the reduced order ccmpensator then G^{s) should 

approximate G. (s) in such a way that o . (T+G„(s)) - 0 . (I+G. (s)) at or 
h min i min h 

below the crossover frequency range and o . (I+G„^(s)) ~ o . (1+G.^(s)) 

min X min h 

at or above the crossover frequency range. This is basically a model 

reduction problem with the objective of matching singular values and 

singular vectors of G. (s) and G. (s) rather than miatchina step responses 

X h 

or other typical model reduction approximation criteria. 

The basic appieal of the frequency domain design methods, namely 
the characteristic loci (CL) and inverse Nyquist array (INA) methods. 


is that they, in many cases, produce lower order and simpler dynamic 
compensators than the LQG approach. However, as was d^ionstrated in 
Chapter 6, they do not autcwnatically produce robust compensators because 
they use diagonalizing matrices, V(s) and W(s), with possibly large 
condition numbers to diagonalize the open-loop plant G^(s) . This implies 
that the gain and phase margins at the physical input and output of the 
system may be much smaller than those in the diagonal coordinate system. 
This difficulty would be eliminated if V(s) and W(s) could be guaranteed 
to have condition numbers near unity and represent stable finite di- 
mensional linear time-invaiiant systems. If V(s) and W(s) are chosen 
as the diagonalizing matrices from the SVD of (s) then they are unitary 
and have condition numbers of unity. However, with these choices for 
V(s) and W(s), neither V(s) nor W(s) represent stable finite dimensional 
linear-time invari.uit systems. Therefore, approximations to this choice 
of Vis) and W(s) must bo sought which approximately diagonalize G^(s). 

This is basically a problejn in realization theory. Further research is 
necessary to determine whether sufficiently accurate approximations of 
V(s) and W(s) (so that (s) is almost diagonal) will produce compensators 
of lower order than compensators obtained by state space design methods 
such as the methodology. If not, then these frequency domain techniques 
would seem to lose their original appeal. 

The above problems (the design of low order robust comiiensators 
via the LQG method or the frequency domain CL or INA methods) assimie 
that the robustness analysis of a control system can be carried out. 

The value of this analysis depends crucially on the control system 


designer's ability to characterize the uncertainty in the open-loop 
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nominal design model in a particular mathematical fashion. Often, the 
characterization of model error by its magnitude alone is insufficient 
to guarantee the stability of the perturbed feedback system by the robust- 
ness tests of Chapter 3. In this case, more detailed information about 
the possible structure of the error transfer matrix, E(s), is needed in 
order to guarantee stability of the feedback system. The robustness 
tests of Chapter 4 require that the magnitude of the projection of the 

error matrix, E(s), onto the siibspace spanned by u (s)v (s) also be 

“n “Ti 

II 

known- It is thus necessary to determine what magnitude of <u (s)v (s) , 

~n — n 

E(s)> constitutes a physically possible type of model error. If this 
model error E(s) gives rise to a perturbed system 5(s) that violates the 
basic physics of the underlying physical process then E(s) can be elimi- 
nated from consideration. However, the control system designer does not 
usually have information about the projections of E(s) onto subspaces and 
what magnitudes of these projections are physically feasible. This 
projection information must be somehow deduced from the information that 
he does have or is able to obtain about the nature of the model error. 

Similar model error cheuracterization problems are encountered if 
the principal gain and phase (PGP) analysis is used. Again, the structure 
of the model error contained in the minimum and maximum principal phases 
of the error matrix must be known. These quantities, like subspace 
projections, are not the kind of information about modelling error 
structure usually possessed by the control system designer. Often the 
type of information available about model error is determined indirectly 
by the knowledge of the acceptable range of parameter variations in a 
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parameterized state-space model. Alternatively, it may be determined 
indirectly from an available set of SISO Nyquist or Bode plots of measured 
transfer functions for certain input 2uid output variedales of the physical 
system under a number of different operating conditions. This type of 
information does not directly provide the model error structure informa- 
tion needed in the robustness analysis. It is therefore necessary to 
devise methods for determining the required structural information about 
model error. This information might be obtained by using a more detailed 
class of models representing truth models of the physical system and 
determining if the resulting model errors for this class of more detailed 
models can ever have projections onto certain subspaces or ever have 
principal phases outside some given range. On the other hand, frequency 
domain measurements in certain input-output directions might obtained 
from the physical system to determine what model error projections or 
principal phases are possible. More experience with practical applications 
is needed in order to detemine how particular types of model error 
structure information may be ascertained. 

Another area for future research is the robustness properties and 
stability margins for time delay systems and multi-sampling rate digital 
systems. LQG based regulators for time delay systems of the form 
x(t) = A^(t) + Ax(t-T) + Bu(t) can be determined from the solution of 
Riccati-like equations. It is not unlikely that these regulators will 
have some inherent robustness properties as does the standard LQ regulator. 
However, continuous-time delay systems are infinite dimensional amd 
mathematically complex and hence it is not clear what type of robustness 
results for these systems can be obtained. A possible approach to these 
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type of prc^lems is to first consider the optimal regulator for discrete 
time delay systems of the form ^ ^1^-n ^ systems 

are finite dimensional and can by state augmentation be formulated as 
standard discrete time systems of the form “ AXj^ + Bu^^. Thus, the 

optimal regulator for a delay system can be reformulated as standard 
optimal regulator problem. The robustness properties of this augmented 
standard optimal regulator may then be related to those of the regulator 
of the time delay system. If this approach proves useful, the analogous 
type of results for optimal regulators of continuous -time delay systems 
may be possibly developed in spite of the fact that these systems are 


infinite dimensional. 


The robustness properties of multirate discrete time systems are 
closely associated with those of discrete time-delay systems and the use 
of multiple sampling rates in discrete time control systems is occurring 
more frequently in practical applications. However, there is a lack 
of robustness theory for this type of control system and most design 
is done on a heuristic basis. Sampled data systems, give rise to 
mathematically complex continuous-discrete hybrid operators if an 
precise description of their continuous time behavior is desired. 

Simple approximations are needed to describe this behavior in order to 
make design use of the robustness results obtained in terms of these 
operators. It would be practically very useful to obtain a method for 
determining the sca;ipling rate in terms of the desired stability margins. 

The mathematical tools used in the analysis of the sampled data 
systems are largely the same as for nonlinear systems. In nonlinear 
system stability analysis, one of the basic problans is not being able 
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characterize the nonlinear systems in a sufficiently simple m<mner that 
is practically useful. There are many stability results like those 
mentioned in section 3.10 but these results use only the grossest 
characterization of the feedback system operators «\cl hence are often 
conservative. To sharpen these results, finer characterizations of the 
feedback system operators must be determined and at present only a few 
results of this type are known. These stability results involve the 
concept of invarijuit limit sets {Lyapunov stability theory) or utilize 
phase plane analysis. They are not part of input-output stability theory 
to which the results and framework of this thesis eire most similar. 

The extension of the robustness results for nonlinear systems analogous to 
the extension of the robustness results for linear systems by use of 
model error structure in Chapter 4 is not possible because there is no 
orthogonal decomposition such as the SVL available for nonlinear systems. 
Other means of determining structure in nonlinear systems must be developed 
and at present it is not clear to the author how this may be accomplished. 
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